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ABSTRACT 


Characteristics  of  gamma-ray  collimators  used  in  clinical 
applications  were  studied  with  the  specific  aims  of  defining  and  of 
simplifying  collimator  design  procedures. 

This  work  involved  both  theoretical  and  experimental  studies. 

A  precision,  mechanical  scanning  device  was  constructed  and  used  in 
field  response  characteristics  measurements  for  various  collimators. 

A  modelling  theory  was  developed  to  determine  analytically  the 
field  response  parameters  for  various  collimator-source  configurations. 
Response  contours  were  then  computer  generated.  Factors  considered  in 
this  modelling  process  were:  air  path  length  attenuation;  tissue  path 
length  attenuation;  possible  bone  absorption,  and  possible  collimator 
penetration. 

The  modelling  algorithm  developed  in  this  work  provided  good 
facsimiles  of  measured  field  isoresponse  contours.  Thus,  using  this 
algorithm,  collimator  characteristics  for  various  source-detector 
configurations  can  be  reasonably  determined  before  fabrication  and 
experimental  measurements  are  made. 
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CHAPTER  ONE 


1.1  Introduction 

Nuclear  medicine  is  an  area  of  medical  practice  where  radio¬ 
isotopes  and  radioactive  materials  are  used  for  the  purpose  of  aiding 
the  diagnosis  of  disease.  In  this  diagnostic  process  the  radioactive 
materials  introduced  into  the  body  are  localized  in  specific  organs  or 
tissues  and  subsequently  counted  or  imaged,  in  normal  practice,  using 
scintillation  radiation  detectors. 

The  static  volume  distribution  of  the  radioactive  material 
in  a  tissue  delineates  anatomic  structure  whereas  changes  in  the 
distribution  or  in  the  quantity  of  radioactive  material  with  time 
demonstrate  organ  or  tissue  function. 

A  scintillation  detector  as  used  in  nuclear  medicine  is 
composed  of  several  components,  most  of  which  are  used  in  conventional 
gamma-ray  spectrometry  systems.  These  are:  a  scintillation  crystal, 
a  photomultiplier,  amplifiers,  an  energy  analyser  and  a  scaler  or 
rate  meter.  A  critical  component  of  a  medical  radiation  detector 
system  not  usually  found  in  a  conventional  spectrometry  system  is 
the  gamma-ray  collimator  which  governs,  to  a  great  extent,  detector 
sensitivity  and  spatial  resolution.  Scintillation  detector  performance 
is  dependent  upon  these  two  factors,  and  optimising  performance 
therefore,  amounts  to  maximizing  detector  sensitivity  for  a  well 
defined  f ield-of-view.  This  raises  certain  problems,  however,  since 
for  reasons  of  patient  safety,  a  minimal  quantity  of  radiation  must  be 
used  in  diagnostic  investigations.  By  demanding  that  the  detector 
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simultaneously  yield  both  high  sensitivity  and  good  spatial  resolution, 
one  is  thus  confronted  with  the  problem  of  reconciling  two  conflicting 
parameters  in  the  design  considerations. 

It  was  the  purpose  of  the  present  work  to  study  gamma-ray 
collimator  characteristics  with  the  specific  aims  of  defining  and  of 
simplifying  collimator  design  procedures  for  various  clinical  measure¬ 
ment  situations. 

The  particular  clinical  situation  that  is  considered  in  this 

thesis  is  the  measurement  of  regional  pulmonary  function  using  an  inert, 

133 

diffusible  radioactive  gas  (  Xe) .  The  scintillation  detector  system 
used  in  these  investigations  of  regional  lung  function  was  designed 
and  constructed  at  this  university  (  1  ) .  Straight  cylindrical 

gamma-ray  collimators  were  used  in  this  instrument  but  the  system 
was  not  optimal  with  respect  to  sensitivity  or  f ield-of-view. 

The  present  work  was  thus  directed  to  studies  of  gamma-ray 
collimator  characteristics  best  suited  to  regional  lung  measurements. 
Specific  aims  of  the  optimization  process  were 

(i)  to  provide  a  well  defined,  non-overlapping  f ield-of-view 
between  detectors, 

(ii)  to  obtain  maximum  sensitivity  for  this  f ield-of-view, 

(iii)  to  produce  a  uniform,  or  iso-sensitive  counting  field 
at  a  distance  from  the  collimator  face. 

1.2  Background 

Two  principal  types  of  scintillation  detector  systems  are 
used  in  nuclear  medicine;  the  counting  system  and  the  imaging  system. 
Historically,  the  first  application  of  a  scintillation  detector  in 
medical  diagnosis  was  as  a  counting  system  (  2  ) . 
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In  terms  of  instrumentation,  counting  systems  and  imaging 
systems  have  much  in  common.  However,  the  purpose  of  each  is  different 
one  is  designed  to  efficiently  measure  the  total  quantity  of  radio¬ 
activity,  the  other  to  efficiently  measure  its  distribution. 

Counting  systems  are  readily  applicable  to  studies  of 
dynamic  organ  function.  As  an  example,  a  bolus  (delta-function)  type 
of  tracer  input  can  be  applied  to  an  organ  or  tissue,  and  measurements 
made  of  the  time-rate  of  change  of  activity  following  this.  In 
similar  fashion,  the  uptake  and  clearance  of  a  radio-labelled 
compound,  which  is  localized  by  biochemical  action,  can  be  measured 
by  such  a  system. 

Imaging  systems  on  the  other  hand  are  in  a  sense 
analogous  to  optical  camera  systems.  Whereas  the  optical  photo  or 
image  indicates  the  distribution  and  intensity  of  the  different  wave¬ 
lengths  of  light  photons  emitted  or  reflected  by  the  object,  the  gamma- 
ray  camera  system  indicates  the  distribution  and  intensity  of  gamma 
photons  emitted  by  the  object. 

Imaging  systems  have  been  developed  based  on  gamma-ray 
emission  ( 5-19 ),  gamma-ray  transmission  (20-2D  ,  x-ray  excitation 
fluorescence  ( 22-2 $  and  incoherent  holographic  (25-27 )  principles .  A 
schematic  summary  of  these  techniques  is  shown  in  Figure  1. 

In  an  imaging  system,  the  observed  distributional  variation 
in  gamma-ray  intensity  results  due  to  the  fact  that  the  radioactive 
material  is  biochemically  accumulated  in  different  concentrations  by 
abnormal  and  surrounding  healthy  tissues.  In  disease  situations  an 
image  thus  contains  "hot”  spots  (a  high  local  concentration)  or  "cold" 
spots  (a  low  local  concentration)  of  radioactivity,  depending  on  the 
type  of  tumour  or  lesion  involved. 


TO  SIGNAL  PROCESSOR 
AND  DISPLAY 
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By  a  suitable  choice  of  scintillation  crystal  size,  good 
sensitivity,  or  detection  efficiency,  can  be  obtained  for  any  given 
gamma-ray  energy.  However,  in  most  medical  applications,  spatial 
resolution  is  also  an  important  consideration  and  gamma-ray  collimator 
design  for  any  particular  application  is  based  upon  a  compromise 
between  the  conflicting  requirements  of  high  spatial  resolution  and 
high  counting  sensitivity.  This  necessary  compromise  in  collimator 
design  has  resulted  in  a  large  variety  of  collimators  in  use,  each 
being  particularly  advantageous  when  applied  to  a  particular  diagnostic 
problem. 

Gamma-ray  collimators  function  by  providing  access  to  the 
detector  for  gamma  rays  only  from  predetermined  directions.  Gamma 
rays  from  other  directions  are  absorbed  or  greatly  attenuated  by  the 
collimator.  The  most  frequently  used  collimator  material  is  lead  - 
a  consequence  of  its  high  density,  high  atomic  number,  relatively  low 
cost  and  ease  of  fabrication.  Tungsten  (  28  ) >  gold  (  29  )  and 
brass  (  30  )  have  also  been  used  as  collimator  materials,  although  not 
as  frequently  as  lead. 

Gamma-ray  collimators  are  specified  in  terms  of  their  physical 

parameters  of  spatial  resolution,  sensitivity,  depth  response,  and 

penetration  fraction  (  31  )  •  In  the  following  review,  the  various 

types  of  collimators  currently  used  in  medical  detection  systems,  and 

their  particular  advantages  and  disadvantages  will  be  briefly  discussed 
(Figure  2). 

(A)  Cylindrical  Straight  Bore  Collimator 

The  first  rectilinear  scanning,  system  built  incorporated  this 
type  of  collimator  (  32  ) •  It  is  inherently  the  most  sensitive  of 
collimators  as  it  exposes  the  greatest  crystal  area  to  a  point  source 
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of  radiation.  However,  the  price  paid  for  this  high  sensitivity  is 
poor  spatial  resolution.  This  is  because  the  solid  angle  subtended 
by  the  crystal  at  the  source  is  large;  consequently  the  points  in  the 


source  from  which  gamma  rays  originate  are  poorly  defined.  These  ever 
conflicting  requirements  of  spatial  resolution  and  sensitivity  are 
better  accommodated  with  a  multichannel  focused  collimator  which  is 
discussed  later.  However,  for  situations  where  the  principal  purpose 
is  to  determine  the  gross  distribution  of  total  body  radioactivity  and 
to  follow  the  distributional  variation  with  time,  the  cylindrical 
straight  bore  collimator  is  still  the  best  device. 

Whole  body  measurements  or  measurements  which  involve  large 
areas  of  the  body  such  as  the  thoracic  volume  employ  this  type  of 
collimator.  By  using  two  opposed  collinear  detectors  with  straight 
cylindrical  collimators,  a  fairly  uniform  region  of  isosensitive  response 
through  the  cross  section  of  a  patient  can  be  obtained. 

(B)  Focused  Collimators. 

Included  in  this  group  are  multichannel  focused  collimators 
and  single  channel  conical  or  tapered  cylindrical  collimators. 

The  focused  multichannel  collimator  generally  gives  the  best 
combination  of  both  high  sensitivity  and  high  spatial  resolution.  The 
reason  is  that  approximately  50%  of  the  crystal  area  is  exposed  through 
the  collimator  holes,  thus  increasing  the  sensitivity;  the  angle  from 
which  the  gamma  rays  can  reach  the  crystal,  however,  is  much  smaller 
than  in  the  parallel  single  hole  collimator.  These  collimators  contain 
many  holes  or  channels,  which  are  both  tapered  and  angled  so  that  their 
combined  geometrical  projections  meet  at  a  focal  point  which  is  a  fixed 
distance  from  the  collimator  face. 
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The  single  hole,  tapered  collimator  (  33  )  is  a  type  of 
"focusing"  collimator,  but  its  characteristics  of  resolution  and 
sensitivity  are  more  related  to  the  straight  cylindrical  hole  collimator. 

It  has  better  spatial  resolution  than  the  straight  hole  collimator 
because  of  its  taper  but  it  is  not  as  sensitive.  On  the  other  hand, 
it  is  more  sensitive  than  the  multihole  focusing  collimator  but  does 
not  resolve  spatial  changes  of  radiation  as  well. 

(C)  Parallel  Multihole  Collimator 

This  collimator  consists  of  a  lead  block  with  1,000,  4,000, 

15,000  or  more  parallel  cylindrical  holes  drilled  into  it.  One  of  its 
main  advantages  is  this  relative  ease  of  fabrication.  These  holes  vary 
in  diameter  from  collimator  to  collimator , making  the  septal  thickness 
suitable  for  the  energy  that  the  collimator  was  designed.  They  are  most 

frequently  used  with  large  crystal  scintillation  cameras  (34).  The  spatial 
resolution  of  a  parallel  multihole  collimator  decreases  steadily  with 

distance  from  its  face.  Sensitivity  is  also  maximal  at  the  face  of 

the  collimator,  so  this  collimator  is  well  suited  for  the  imaging 

of  shallow-lying  lesions  and  tumors.  It  is  used  primarily  for  low 

energy  gamma  photons  because  a  parallel  multihole  collimator,  if  not 

properly  designed,  is  susceptible  to  septal  penetration  from  high 

energy  radiation,  resulting  in  poorly  defined  images. 

(D)  Pinhole  Collimator 

The  pinhole  collimator  is  the  most  analogous  to  an  optical 
system  —  that  of  a  pinhole  camera.  It  was  also  one  of  the  first 
collimators  used  for  medical  applications  (35  ) .  The  image  is 
inverted,  and  depending  on  the  distance  of  the  source  from  the 
collimator,  it  may  be  enlarged  or  diminished  in  size,  or  at  a  particular 
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distance  may  be  equal  to  the  object  size.  Used  primarily  for  viewing 
small  organs,  the  pinhole  collimator  provides  sensitivity  approximately 
equal  to  that  of  a  straight  cylindrical  hole  collimator,  and  in  addition 
gives  improved  resolution  due  to  the  "optical  gain"  effect  of  image 
enlargement . 

(E)  Rectangular  Collimators 

Each  type  of  collimator  considered  to  this  point,  whether 
single  hole  or  multihole,  straight  or  tapered,  can  have  an  analogous 
rectangular  counterpart.  Generally,  the  sensitivities  and  resolutions 
of  these  rectangular  collimators  are  similar  to  the  corresponding 
cylindrical  collimators,  the  only  real  difference  being  in  the  field 
of  view.  Instead  of  a  circular  field  of  view  from  a  cylindrical 
collimator,  a  "large  rectangular"  field  of  view  is  obtained. 

1.3  The  Pulmonary  Function  Studies  System  (University  of  Mberta) 

The  Pulmonary  Function  Studies  System  (PFSS)  is  a  multi¬ 
detector  counting  system  used  to  measure  regional  lung  function  using 
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the  radioactive  gas  Xenon  (  1  ) .  Certain  basic  physiological 

parameters  of  lung  function  can  be  measured  and  related  to  lung 
disease.  (For  a  detailed  explanation  of  these  basic  lung  parameters, 
see  Appendix  A) . 

Each  lung  in  the  PFSS  is  viewed  by  6  pairs  of  collinear 
detectors,  one  detector  from  each  pair  mounted  anteriorly  and  the  other 
mounted  posteriorly  (See  Figure  3).  It  is  the  purpose  of  the  present 
work  to  study  collimator  design  with  particular  reference  to  this 
system  and  to  develop  a  collimator  which  will  simultaneously  optimize 
sensitivity  and  uniformity  of  depth  response  and  minimize  field  of 
view  overlap  between  the  detector  pairs. 


[iiiiwam 
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This  work  involves  a  combination  of  both  theory  and  experiment. 

A  modelling  theory  is  developed  and  outlined,  and  a  number  of  experimental 
studies  are  done.  The  theory  is  first  discussed  -  theoretical  aspects 
of  collimator  design,  with  reference  to  the  various  parameters 
(sensitivity,  spatial  resolution,  etc.)  are  considered  in  detail  in 
the  following  chapter. 
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CHAPTER  TWO 


THEORETICAL  CONSIDERATIONS  IN  COLLIMATOR  DESIGN 

In  any  comparison  between  collimators,  sensitivity,  resolution, 
and  field  response  need  to  be  determined  for  each.  To  a  large  extent, 
these  parameters  can  be  obtained  from  geometrical  considerations  for 
the  collimators.  It  is  important  to  realize,  however,  that  radiation 
scatter,  penetration  effects,  and  the  uncertainty  of  a  source 
distribution  can  considerably  alter  the  effectiveness  of  a  collimation 
system.  This  chapter  deals  with  these  factors  and  considers  their 
individual  effects  on  detector  performance. 

2.1  Sensitivity 

In  any  radionuclide  detection  system,  the  number  of  gamma-rays 
detected  will  always  be  less  than  the  number  emitted  by  the  source, 
which  can  be  distributed  over  any  volume  geometry.  Since  a  detection 
system  is  made  up  of  several  components,  one  has  to  look  at  the 
influence  of  each  of  these  separately.  Following  Mallard  (36  ) »  we 
can  define  the  fractional  efficiency  as: 

£  _  number  of  output  events 

i  number  of  input  events 

for  device  i.  It  is  obvious  that  the  complete  detection  system  has  an 
overall  efficiency  equal  to  the  product  of  all  the  component  efficiencies 

£  =  ne.  (2.1) 

i  1 
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In  clinical  work,  the  source  is  generally  distributed  over  a 
substantial  volume,  and  thus  the  amount  of  radiation  or  activity  that  is 
viewed  by  the  detector  is  usually  difficult  to  determine  analytically, 
if  not  impossible.  It  is  therefore  appropriate  in  sensitivity 

considerations  to  differentiate  the  source  into  infinitesimal  point 
sources  and  then  integrate  the  individual  responses  to  each  of  these 
sources  over  the  whole  of  the  volume  distribution. 

If  Pp  is  the  number  of  gamma-ray  photons  emitted  per  unit  time 
into  a  4tt  angular  area  by  a  point  source  and  the  exposed  regions  of 
the  detector  subtend  a  fractional  solid  angle  Q  at  the  source,  then 
PpSTi!  is  the  number  of  photons  incident  on  the  detector  per  unit  time. 

If  the  detection  efficiency  of  the  crystal-photomultiplier  is  2^,  then 
the  counting  rate  will  be: 

CI  =  ppq?d  (2-2> 

assuming  that  the  other  components  of  the  counting  system  are  ideal, 
and  neglecting  scatter  and  collimator  penetration.  If  these  latter  two 
effects  are  included,  then  the  practical  counting  rate  will  be  larger 
than  the  ideal  one.  If  Cq  is  defined  as  the  count  rate  due  to  scattering 
alone  and  as  that  due  to  penetration  alone,  then  the  true  count  rate, 
C-p,  is  given  as: 

Cp  =  Cp  +  Cg  +  Cp  (2 . 3) 

The  fractional  increases  in  count  rate  due  to  scatter  and  penetration 
may  be  defined  as: 


fs  "  cs/ci 


f  p  —  Cp/ Cp 


(2.4) 
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Thus  (2.3)  may  be  rewritten : 

CT  =  pp^D  (l+fs+fp)  (2*5) 

In  (2.5),  only  Pp  is  directly  dependent  on  the  source,  all 
other  terms  being  concerned  with  the  performance  of  the  system.  (Although 
fg  ,  fp>  are  obviously  dependent  on  the  energy  of  ijie  photons,  they 
are  still  governed  by  the  type,  thickness  and  material  of  the  detector 
crystal  and  collimator).  Thus  the  point  source  sensitivity  (s  )  of 
the  system  may  be  defined  as 

sp  =  fi  (l+fg+fp)  €d  (2 . 6) 

and 


C-p  -  Sp  Pp  (2.7) 

It  is  found,  however,  that  a  line  source  of  radioactivity  is 
much  more  practical  to  fabricate  and  work  with  experimentally  than  a 
point  source;  it  will  also  be  shown  that  it  yields  more  information. 

If  Cp  is  the  measured  count  rate  from  a  line  source  of 
radioactivity  which  emits  gamma-ray  photons/unit  length/unit  time, 


C  =S  P 
L  L  L 


(2.8) 


and  regarding  the  line  source  as  an  infinite  number  of  point  sources, 


SL  =  /Spdy  (2.9) 

y 

where  is  the  line  sensitivity  of  the  system.  The  graph  of  line 
sensitivity  versus  transverse  displacement  of  the  line  source  from  the 
collimator  axis  is  called  a  line  spread  function  (LSF) .  (See  Figure  4  ) • 
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In  order  to  obtain  an  LSF,  the  line  source  must  be  moved 


transversely  across  the  field  of  view;  it  thus  traces  out  a  plane.  In 
this  way  a  plane  source  sensitivity  can  be  obtained: 


S 


(2.10) 


A 


P. 


L 


(2.11) 


P 


P 


Using  (2.10),  the  plane  sensitivity  is  therefore  determined  by  the  area 
enclosed  by  the  LSF  curve  in  Figure  4. 


The  plane  sensitivity  is  valuable  in  comparing  detection 


systems  because  the  response  to  a  plane  source  is  independent  of  the 
distance  from  the  detector,  provided  the  field  of  view  is  uniformly 
covered  by  the  source.  This  is  because  the  increase  in  the  area  of  the 
plane  source  in  the  field  of  view  as  the  source  -  detector  distance  is 
increased  is  compensated  by  the  decreased  response  to  each  of  the 
infinitesimal  point  sources  on  the  plane  (  77  ) .  Practically,  however, 
scatter  and  penetration  have  an  effect  on  the  plane  sensitivity  but 
the  variations  with  distance  are  slight  compared  with  other  source 
configurations.  Attenuation  of  the  gamma  rays  through  the  medium  has 
also  been  neglected. 

The  LSF  is  also  a  valuable  tool  in  defining  the  spatial 
resolution  of  a  detection  system.  This  parameter  is  discussed  below. 

2.2  Spatial  Resolution 

Spatial  resolution  is  a  measure  of  a  system's  ability  to  detect 


spatial  changes  of  gamma— ray  intensity  from  a  radioactive  source 
distribution.  The  optical  analogue  of  resolution  expresses  the 
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minimum  separation  between  two  objects  that  can  still  be  distinguishable 
as  separate  in  the  image. 

It  has  become  common  practice  to  describe  resolution  by 
considering  the  response  to  sinusoidal  distributions  of  radioactivity. 

In  this  way,  any  object  or  source  distribution  can  be  described  by  a 
set  of  spatial  frequencies.  A  diffuse  object  with  poor  definition  is 
thus  described  by  a  spectrum  of  low  spatial  frequencies  while  an  object 
that  has  sharply  defined  edges  requires  in  addition  high  frequency 
components  to  describe  it.  A  mathematical  analogue  is  found  in 
Fourier  series  and  transforms  used  to  approximate  edge  and  step 
functions . 

Sources,  or  more  meaningfully,  source  edges,  can  be  regarded 

as  alternating  "activity  amplitudes"  with  maximum  activity  p  and 

max 

minimum  activity  p  ^  .  The  source  modulation  is  then  defined  as 

m  =  ^max  -  Pmin  (2.12) 

J  i 

P  +  p  . 

max  mm 

The  imaging  system,  in  turn,  produces  an  image  of  this  modulation, 
which  is  called  the  image  modulation: 

m  =  .  ^max  -  Cmin  (2.13) 

1  C  +  C  . 

max  mm 

where  C  and  C  are  the  minimum  and  maximum  counting  rates,  respec- 
min  max 

tively.  The  ability  of  the  detection  system  in  transferring  the  source 
modulation  to  an  image  modulation  is  represented  by  the  modulation 
transfer  function  (MTF) . 


MTF  = 


(2.14) 


' 
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The  source  may  be  regarded  as  made  up  of  numerous  line  elements  and 

therefore  the  response  of  the  detection  system  may  be  calculated  by 

the  convolution  of  each  element  with  the  LSF  of  the  system.  Thus, 

from  (2.14):  . +x ■ 

J  _x  S  (x)  cos(27TVx)dx 

MTF(V)  =  o _ (2.15) 

f+xo 

J  _x  SL(x)  dx 

Mallard  (30') 

It  is  evident  from  this  consideration  that  the  MTF  is  pro¬ 
portional  to  the  Fourier  transform  of  the  LSF.  The  variable  V  is  the 
spatial  frequency  parameter  and  has  dimensions  length  ^  (See  Figure  4). 

Although  the  MTF  is  a  valuable  index  of  spatial  resolution, 
in  a  detection  system  where  resolution  is  not  the  prime  criterion  for 
design  considerations,  the  MTF  need  not  be  taken  into  account.  In 
single  hole  collimator  counting  systems  resolution  is  degraded  to 
provide  for  increased  sensitivity.  "Visualization”  is  not  required 
and  so  the  index  of  resolution  refers  mainly  to  the  field  of  view.  In 
the  PFSS  this  determination  of  the  field  of  view  is  an  important 
consideration,  since  regional  lung  function  is  studied  and  overlap 
from  adjacent  regions  will  cause  a  "blurring"  effect  in  the  data  analysis. 
The  field  of  view  is  determined  solely  from  collimator  geometry , although 
edge  penetration  will  increase  the  calculated  values  somewhat.  Appendix 
B  derives  the  configurations  of  the  fields  of  view  from  cylindrical 
and  rectangular  collimators.  It  is  shown  that  the  field  of  view  area 
increases  as  the  square  of  the  distance  from  the  collimator  face. 

2.3  Field  Response 

The  field  response  of  a  collimator  system  is  valuable  in 
determining  both  the  spatial  resolution  and  sensitivity  of  that  system. 
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It  is  defined  as  the  response  of  a  detection  system  to  a  point  source 
of  radiation  located  anywhere  in  space.  Practically,  it  is  measured 
as  the  count  rate  of  a  system  as  a  function  of  source  point  position. 
This  parameter  was  measured  for  a  variety  of  collimators  and  the 
results  are  presented  in  Chapter  4.  Analytically,  the  field  response, 
F(x,y,z),  is  Ct  from  equation  (2.5)  measured  over  a  source  volume: 

F(x,y,z)  =  pp  f)  (x,y ,  z)  £D(x,y,z)  [1  +  fg(x,y,z)  +  f  (x,y,z)] 

(2.16) 

F(x,y,z)  is  usually  normalized  to  C-p  .  If  all  points  (x,y,z)  that 

max 

satisfy  F(x,y,z)  =  constant  are  joined  together,  an  isoresponse  or 
isocount  response  line  is  formed.  When  several  values  of  constants 
are  used,  isoresponse  contours  are  generated.  The  depth  of  the  field 
response  contours  indicates  the  sensitivity  of  the  detection  system  as 
a  function  of  distance;  the  broadness  or  width  of  the  isoresponse 
lines  indicates  spatial  resolution  capabilities  of  the  system.  By 
comparing  these  isorespanse  curves  from  different  collimators,  and 
collimator  systems,  both  sensitivity  and  resolution  can  be  optimized 
for  particular  purposes.  The  methods  used  to  measure  field  response 
are  considered  in  Chapter  3  -  Experimental  Design. 
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CHAPTER  THREE 


EXPERIMENTAL  DESIGN 

3.1  Instrumentation 

3.1.1  Raster  Scanner 

In  order  to  measure  field  response  of  the  several  collimators 
studied,  a  volume  source  distribution  was  needed.  It  was  decided  that 
the  required  source  distribution  could  be  most  easily  attained  by 
integrating  a  point  source  of  radiation  over  the  appropriate  volume. 

In  practice,  this  integration  is  achieved  by  mechanically  "scanning" 
the  point  source  over  the  volume.  Major  design  criteria  for  this 
scanner  were  the  precision  and  accuracy  required  in  source  -  detector 
positioning.  Commercial  scanning  devices  were  initially  considered 
for  these  studies,  but  none  was  available  that  could  be  modified  to 
suit  our  requirements.  The  adaptation  of  a  micro-manipulation  device, 
such  as  is  used  in  scanning  microscopy,  was  then  considered  as  a  basis 
for  our  scanning  unit.  However,  the  available  range  of  movement  with 
these  devices  is  generally  limited  to  about  5.0  cm  on  a  side,  which 
was  insufficient  for  our  proposed  measurements  of  about  25  cm  a  side 
(an  approximation  to  the  thorax  dimensions) . 

The  cost  of  these  several  commercial  devices  was  also  high, 
thus  it  was  decided  to  undertake  the  design  and  construction  of  the 
raster  scanner  for  this  particular  application.  The  mechanical 
arrangement  in  this  design  is  shown  in  Figure  5. 

Two  Hers t  synchronous  motors  drive  the  source  platform  at  a 
constant  speed  of  0.25  mm/sec  in  x  or  y  directions  (Figure  6  ). 
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FIGURE  5 
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The  maximum  range  of  the  scanner  is  25  cm  in  the  x  direction  and  26  cm 
in  the  y  direction.  The  ranges  can  be  varied  within  these  maximum 
values  by  means  of  aluminum  ’'stops”  and  microswitches. 

An  electronic  logic  device  was  designed  to  control  the  motion 
of  the  scanner  (Figure  7  ).  In  operation,  the  scanning  system  is 
initialized  with  the  source  platform  at  (x,y)  equal  to  (0,0) .  A 
pushbutton  "start”  provides  power  to  the  x  motor  which  drives  the 
platform  to  the  end  of  the  desired  x  scan  length.  A  microswitch 
at  this  end-point  provides  a  signal  to  stop  the  x-drive  motor  and  start 
the  y-drive  motor  for  a  preset  time.  This  preset  "on-time”  of  the 
y  motor  is  controlled  by  means  of  a  variable  resistor  and  ranges  between 
40  nsec  to  40  sec  in  duration.  For  the  scans  reported  in  this  thesis 
the  y  motor  was  set  to  be  driven  for  10  seconds,  during  which  time 
the  platform  and  x  motor  assembly  moved  2.5  mm.  At  the  end  of  this 
preset  time  the  y  motor  is  shut  off  and  the  x  motor  started  in  reverse, 
bringing  the  platform  back  to  the  initial  x  value.  At  this  point  it 
is  again  stopped  by  a  microswitch,  and  the  y  motor  driven  for  another 
period  of  10  seconds.  One  cycle  of  this  system  is  completed  after  this 
latter  10  second  period  when  the  y  motor  is  shut  off. 

In  this  way,  the  source  platform  traverses  the  area  desired 

in  a  raster  pattern.  The  scanning  assembly  stands  at  a  height  of 
about  15  cm  with  the  collimator-crystal  combinations  mounted  underneath 
the  motor  platform  and  the  source  suspended  from  the  platform  on  a 
rigid  rod.  The  source  was  positioned  in  the  same  plane  as  the  detector 
longitudinal  axis  so  that  the  (assumed)  symmetrical  response  of  the 
crystal  could  be  used  for  calculation  purposes. 
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3.1.2  Data  Aquisition  System 

Each  detector  used  was  coupled  to  a  pulse-shaping  preamplifier. 

The  pulse  from  the  preamplifier  was  negative,  2  ysec  in  length  and 
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nominally  10  mv  in  height  for  the  80  keV  photopeak  of  Xe.  The 
signal  was  then  routed  to  an  amplifier  (HAMMER  RC  NA-11)  and  amplified 
to  a  2-8  volt  range  (depending  on  the  energy  of  the  original  incident 
gamma  ray) .It  was  then  routed  to  a  single  channel  analyser  which  was 
set  to  accept  pulses  corresponding  to  an  energy  range  equal  to  that 
of  the  incident  gamma  rays.  If  the  input  pulse  was  in  the  correct 
energy  "window"  then  a  positive  10  volt  pulse  was  provided  to  a  scaler 
to  record  the  event.  Six  10  MHz  scalers  were  linked  in  a  "daisy  chain" 
with  a  timer-scaler,  and  magnetic  tape  controller  for  a  CIPHER  seven 
track,  write-only  tape  transport.  The  data  on  the  scalers  could  thus 
be  stored  on  magnetic  tape  after  a  preset  data  gathering  period.  After 
each  record  was  written  on  the  tape,  the  scalers  were  initialized  and 
a  new  counting  period  was  begun.  A  "record"  on  the  tape  was  composed 
of  a  42  digit  BCD  number.  Each  group  of  six  digits  comprised  the 
count  from  one  scaler,  and  therefore,  each  individual  scaler  could  be 
identified. 

At  the  end  of  each  x-scan,  by  means  of  the  microswitch  on 
the  source  platform,  a  positive  5  volt  pulse  was  sent  to  a  "y  position" 
scaler.  By  inspection  of  this  scaler,  the  position  of  the  source  could 
be  determined  for  any  time  interval,  or  for  any  particular  data  record. 

At  the  end  of  a  complete  scan,  a  microswitch  provided  a  pulse  to  the 
appropriate  scaler,  thus  "labelling"  the  termination  of  the  scan.  In 
the  studies  reported  here,  five  scalers  were  used  for  data  acquisition 
and  coordinate  position  determination.  The  scalers  in  the  "daisy  chain* 
not  in  use  were  recorded  on  tape  as  having  a  zero  count. 
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The  data  collection  period  was  chosen  to  be  10  seconds,  during 
which  time  the  source  moved  2.5  mm.  At  the  end  of  the  data  acquisition 
period,  the  contents  of  the  several  scalers  were  written  on  tape  and  the 
scaling  system  reset.  Each  record  was  thus  10  seconds  apart  in  time 
and  2.5  mm  apart  in  space. 

In  these  studies  a  normal  scan  of  20  cm  by  23  cm  required 
about  20  hours  to  complete  and  comprised  some  7500  individual  data 
records  on  tape.  Data  "blocking”  could  not  be  achieved  with  the 
magnetic  tape  unit  available  to  this  investigation,  thus  each  scan 
occupied  about  150  meters  of  tape.  These  data  tapes  were  subsequently 
blocked  and  rewritten  prior  to  processing  using  the  IBM  360/67.  Under 
normal  circumstances  about  20  individual  scans  were  stored  on  one 
180  meter  tape  reel. 

3.1.3  Data  Analysis 

The  data  was  read  from  tape  and  all  records  stored  in  core, 
each  record  being  corrected  for  radioactive  decay.  A  search  was  made 
for  end-of-x-scan  points  by  *  looking  for  two  consecutive  non-zero 
records  from  the  particular  scaler  attached  to  the  x-scan  microswitch. 

By  comparing  x-scan  lengths,  artefacts  due  to  line  noise  could  be 
rejected,  and  the  data  then  segmented  into  a  number  of  one  dimensional 
scans  as  detailed  schematically  in  Figure  8  . 

The  maximum  count  in  the  scan  matrix  was  then  found  and  all 
records  normalized  to  this  value.  This  normalization  procedure  reduced 
the  data  from  "count  rate"  to  fractions  of  the  maximum  count  rate.  A 
plot  file  was  then  prepared  by  storing  these  fractions  in  a  matrix 
array,  the  individual  matrix  elements  corresponding  one  to  one  to  an 
individual  record  in  time  and  space.  When  the  matrix  was  completed  a 


—not*— ico 


27 


* 


'D"  MICROSWITCH  TRIPS 


ONE  X-SCAN-LENGTH 


ONE  X-SCAN-LENGTH 


10  SECS  APART  IN  TIME 


J 


D"  MICROSWITCH  TRIPS 

=>  TWO  CONSECUTIVE 
RECORDS  WITH  NON-ZERO 
ENTRIES 


"D" 

MICROSWITCH 
ATTACHED  TO 
SCALER 


FIGURE  8 


28 


least-squares  curve-fitting  technique  joined  equal-valued  records  into 
an  isocount  response  curve.  This  plot  was  stored  into  a  final  plot 
file  and  an  isoresponse  contour  map  drawn  by  a  Calcomp  plotter.  The 
number  of  contours  in  this  map  and  their  values  could  be  arbitrarily 
chosen  by  the  operator  and  could  be  entered  as  parameters  into  the 
program  via  keyboard. 

The  interactive  program  (Appendix  E)  was  written  in  Fortran  IV 
and  ran  under  the  MTS  operating  system  of  the  IBM  360/67  computer. 
Similar  programs  were  used  to  analyse  scans  from  one,  two,  or  three 
detectors . 

3.2  Modelling 

In  addition  to  making  experimental  measurements  of  field 
response  for  various  source-detector  configurations  it  was  decided  to 
also  simulate  these  configurations  and  determine  field  response  using 
the  theoretical  considerations  discussed  in  Chapter  2.  The  purpose 
of  this  simulation  was  to  determine  the  validity  of  the  theoretical 
specifications  of  resolution  and  sensitivity,  with  a  direct  comparison 
to  the  experimental  results.  Also,  using  this  "modelling"  theory, 
lung  configurations  encountered  by  the  PFSS  detectors  which  could  not 
be  measured  experimentally  could  be  simulated.  For  example,  rib 
absorbers  were  introduced  into  the  modelled  program  and  corresponding 
field  responses  calculated;  this  could  not  be  done  with  the  experimental 
scanning  system.  The  approach  taken  in  devising  the  modelling  algorithm 
follows . 

As  discussed  in  Chapter  2,  the  simplest  way  to  model  a  source 
distribution  is  to  calculate  point  responses  from  a  number  of  equally 
spaced  points  in  the  source  volume  and  integrate  these  responses.  With 
the  raster  scanner,  a  plane  source  response  is  measured.  This  plane  is 
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parallel  to  the  longitudinal  axis  of  the  cylindrical  crystal.  Making  the 
approximation  that  the  crystal  response  is  circularly  symmetric,  the  plane 
response,  can  be  rotated  about  the  crystal  axis  and  a  three-dimensional 
volume  response  obtained.  This  method,  however,  assumes  the  source 
medium  is  homogeneous  in  absorption  characteristics.  In  most  clinical 
situations  this  is  not  true.  For  example  in  the  lung  (as  viewed  by  the 
PFSS)  the  location  from  which  a  gamma-ray  is  emitted  from  the  interior 
of  an  alveolus,  and  its  direction,  govern  the  amount  of  air  and  tissue 
absorber  that  it  must  pass  through  before  reaching  the  detector.  Part 
of  this  path  might  also  lie  through  a  rib. 

A  model  is  needed  that  accounts  for  these  situations  of 
absorber  inhomogeneity.  The  approach  taken  was  the  use  of  a  ray¬ 
tracing  grid  system.  The  source  plane  was  broken  up  into  an  array  of 
source  points.  The  crystal  face  was  also  replaced  by  a  grid  system 
of  points,  lying  in  the  z=0  plane  in  normal  3-space.  The  collimators, 
absorber  configurations  (such  as  ribs) ,  and  tissue  planes  were  described 
analytically.  Each  point  on  each  grid  had  a  unique  set  of  cartesian 
rectangular  coordinates  assigned  to  it.  The  detector  grid  (DG)  was 
used  as  the  reference  grid.  From  each  point  (Pp)  on  the  DG,  lines  were 
generated  to  all  points  (P<.)  on  the  source  grid  (SG)  .  The  equations 
of  the  lines  were  determined,  and  the  various  collimator,  bone,  and 
tissue  attenuation  distances,  if  any,  calculated.  The  theory  of  this 
ray  tracing  method  is  detailed  in  Appendix  C. 

Besides  considering  attenuation  of  the  gamma  rays  by  bone 
and  tissue,  and  also  the  attenuation  due  to  possible  collimator 
penetration,  the  solid  angle  factor  for  each  detector— source  point 
and  each  detector  point-source  point  arrangement  had  also  to  be  determined. 
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The  calculation  of  these  solid  angles  is  detailed  in  Appendix  D.  The 
method  used  is  that  described  by  Gardner  and  Verghese  (  38  )  and 
consists  of  replacing  the  crystal  disk  by  an  equal  area  right  polygon 
with  total  number  of  sides  n,  where  n  is  even. 

This  solid  angle,  Q  (s,R,p),  where  &  is  the  height  of 
the  source  point  above  the  crystal,  R  the  radius  of  the  crystal, 
and  p  the  distance  between  the  center  of  the  crystal  and  the  projection 
of  the  source  point  on  the  crystal  plane,  must  be  calculated  for  each 
source  grid  point.  Let 


ft.  =  ft  (S,R, p)  (3.1) 

for  source  point  i  with  coordinates  (xg,  y^,  z^) .  ^  will  be  the  solid 

angle  subtended  by  the  total  crystal  face  from  point  i.  However,  the 
crystal  will  be  broken  up  into  an  array,  and  therefore,  each  detector 
point  must  be  considered  separately.  This  is  achieved  by  weighting 
all  grid  points,  relative  to  each  other,  with  respect  to  the  solid 
angle.  Since  a  solid  angle  depends  on  the  area  exposed  and  the 
distance  from  the  subtended  point,  these  were  the  criteria  chosen  to 
determine  the  weighting  factors.  For  a  detector  point  j  and  a  source 
point  i  (Figure  9),  the  weighting  factor  becomes: 

w.  .  =  cos  0 .  .  (3.2) 

13  _ 11 

2 

r  .  . 
ij 

If  the  area  of  each  square  on  the  detector  is  identical,  then  the 
relative  area  exposed  by  each  square  will  be  directly  proportional  to 
the  angle  from  the  vertical  that  the  ray  ij  makes.  Hence  the  need  for 
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the  cosine  term  in  (3.2).  Now,  define 


wii°i 


Zw 

k 


ik 


(3.3) 


Q^j  is  the  assigned  relative  fraction  of  the  total  solid  angle  fi. 
for  the  detector  point  j .  In  summarizing,  if  the  source  point  i  emits 
gamma-rays  per  unit  time  into  a  4tt  geometry,  then  X^  gamma  photons 
per  unit  time  are  incident  on  the  crystal.  The  detector  point  j  receives 
^ijXi  these.  Of  course  the  beam  is  further  attenuated  by  tissue,  air, 
and  possibly  bone  and  collimator  material.  Thus  the  portion  of  activity 
from  point  i  that  point  j  detects  is 


A.  . 
1J 


-y.t . .  -y  a .  .  ■ 
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T)  ij  ^c 
e 


(3.4) 


where  y.  y  ,  y,  ,  y  are  the  linear  attenuation  coefficients  for 
t  a  b  c 

tissue,  air,  bone,  and  collimator  material  respectively  for  a  given 


gamma-ray  energy  and  t..,  a  .,  b .  . c..  are  the  corresponding  path 

ij  ij  ij  ij 

lengths  through  these  materials  for  the  ray  ij .  By  summing  A _  over 

either  of  its  indices,  interesting  values  are  obtained.  ZA. .  is  the 

.  ij 
l 


total  amount  of  radiation  that  reaches  detector  point  j.  ZA. .  is  the 

j  1J 

amount  of  activity  that  is  detected  by  the  total  crystal  from  a  source 
point  i,  and  is  equal  to  the  field  response  parameter  F(x,y,z)  from 
equation  (2.16).  For  a  raster  scanner  simulation,  the  value  ZA..  is 

j  1J 

stored  in  a  location  corresponding  to  source  position  i  (x^,  yc,Zg)  and 
after  all  rays  have  been  traced,  identical  values  are  joined  to  form 


isocount  response  curves. 
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Most  configurations  scanned  were  also  modelled,  as  described 
above;  the  algorithms  used  are  described  in  Appendices  C  and  D. 

The  results  of  both  the  experimental  and  modelling  studies 
are  presented  in  the  next  section. 


CHAPTER  FOUR 


RESULTS  AND  DISCUSSION 

4.1  Introduction 

Four  types  of  collimators  were  studied:  (1)  straight  cylindrical, 
(2)  tapered  cylindrical,  (3)  straight  rectangular  and  (4)  tapered 
rectangular. 

A  total  of  twenty-eight  scans  were  done.  Each  collimator 

arrangement  was  assessed  for  a  photon  source  in  both  air  and  water.  The 

22  133 

responses  for  both  Na  (551  keV)  and  Xe  (80  keV)  were  determined. 

The  experimental  source-detector  configurations  studied  are  summarized 
in  Table  1. 

The  experimental  studies  were  carried  out  in  parallel  with 
theoretical  simulations  of  the  same  source-detector  geometries,  employing 
the  modelling  theory  developed  earlier.  These  simulated  isocount  contours 
were  compared  with  the  experimental  results  in  order  to  determine  the 
range  of  applicability  of  the  theoretical  model.  In  this  way,  response 
curves  from  source-detector  geometries  that  we  were  unable  to  duplicate 
experimentally  could  be  obtained  from  the  theoretical  model. 

The  ray-tracing  algorithm  derived  was  general  and  allowed  the 
study  of  several  source-attenuator-collimator-crystal  models.  Both 
homogeneous  and  inhomogeneous  distributions  of  radioactive  source  and 
attenuating  media  were  evaluated.  Source  to  tissue  surface  distance 
and  tissue  surface  to  collimator  and  crystal  distance  could  be  varied. 

The  gamma  ray  energy  considered  was  also  varied  as  required  for  any  one 
model  but  monoenergetic  sources  were  assumed.  Collimator  type  (e.g. 
cylindrical,  conical,  rectangular  )  and  construction  material  could  also 
be  changed  as  required. 
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TABLE  1  -  Experimental  Detector  -  Source  Configurations 


COLLIMATOR 

SOURCE 

ENERGY 

(keV/ 

MEDIUM 

PAGE  NO. 

CC 

133 

Xe 

80 

AIR 

40  ,  41 

CS 

133 

Xe 

80 

AIR 

40  ,  41 

CD 

133 

Xe 

80 

AIR 

40  ,  41 

R-0° 

133 

Xe 

80 

AIR 

42  ,  43 

R-30° 

133 

Xe 

80 

AIR 

42  ,  43 

R-45° 

133y 

Xe 

80 

AIR 

42  ,  43 

R-60° 

133 

Xe 

80 

AIR 

42  ,  43 

CC 

133 

Xe 

80 

H2° 

48  ,  49 

CS 

133 

Xe 

80 

h20 

48  ,  49 

CD 

133 

Xe 

80 

h2o 

48  ,  49 

R-0° 

133 

Xe 

80 

H20 

50  ,  51 

R-30° 

133v 

Xe 

80 

V 

50  ,  51 

R-45° 

133 

Xe 

80 

h2o 

50  ,  51 

R-60° 

133 

Xe 

80 

h2° 

50  ,  51 

CC 

22 

Na 

511 

AIR 

56  ,  57 

CS 

22 

Na 

511 

AIR 

56  ,  57 

CD 

22 

Na 

511 

AIR 

56  ,  57 

R-0° 

22 

Na 

511 

AIR 

58  ,  59 

R-30° 

22 

Na 

511 

AIR 

58  ,  59 

R-45° 

22 

Na 

511 

AIR 

58  ,  59 

R-60° 

22m 

Na 

511 

AIR 

58  ,  59 

CC 

22 

Na 

511 

h2o 

64  ,  65 

CS 

22 

Na 

511 

h2o 

64  ,  65 

CD 

22 

Na 

511 

h20 

64  ,  65 
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TABLE  1  -  Continued 


COLLIMATOR 

SOURCE 

ENERGY 
(keV ) 

MEDIUM 

PAGE  NO 

R-0° 

22 

Na 

511 

H2° 

66  ,  67 

R-30° 

22.. 

Na 

511 

H2° 

66  ,  67 

R-45° 

22 

Na 

511 

h2o 

66  ,  67 

R-60° 

22 

Na 

511 

H2° 

66  ,  67 

Miniature 

133 

Xe 

80 

AIR 

72 

Miniature 

!33y 

Xe 

80 

H2° 

72 

Miniature 

133 

Xe 

80 

H20  +  SKUL 

L  72 

Coincidence 

22m 

Na 

511 

AIR 

73 

Coincidence 

22„ 

Na 

511 

h2o 

74 

Legend : 


CC 

-  Cylindrical 

CS 

-  Cylindrical 

CD 

-  Cylindrical 

R-0° 

-  Rectangular 

R-30° 

-  Rectangular 

R-45° 

-  Rectangular 

R-60° 

-  Rectangular 

Converging 

Straight 

Diverging 

-  0°  Taper  from  Normal 

-  30°  Taper  from  Normal 

-  45°  Taper  from  Normal 

-  60°  Taper  from  Normal 


( 
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4.2  Details  of  Collimators  Studied 

The  straight  cylindrical  collimator  was  made  of  lead,  6  mm  thick. 
The  crystals  used  were  25.4  mm  diameter,  12.7  mm  thick  NaI(T£)  (Harshaw 
Chemical  Company),  and  were  placed  7.6  cm  back  from  the  face  of  the 
collimators . 

The  tapered  cylindrical  collimator  was  machined  from  a  stock 
brass  cylinder,  6.35  cm  in  diameter.  The  inside  diameters  of  this 
collimator  were  3.175  cm  and  1.905  cm  respectively;  the  total  length  of 
the  collimator  was  7.6  cm  and  taper  angle  4.76°.  This  collimator  was 
used  in  both  converging  and  diverging  directions. 

The  rectangular  collimator  was  fabricated  from  0.635  cm  cold- 
rolled  steel  plate,  and  its  dimensions  are  shown  in  Figure  10  . 

One  side  was  designed  to  swivel  around  its  end,  permitting 
measurements  to  be  made  of  a  tapered  rectangular  collimator,  with 
varying  degrees  of  taper. 

A  miniature  detector  system,  comprising  a  0.6  cm  diameter,  1.27  cm 
thick  NaI(T&)  crystal  was  also  studied.  The  crystal  was  recessed  9  cm 
from  the  face  of  a  straight,  cylindrical,  stainless  steel  and  lead 
collimator  1.27  cm  internal  diameter.  In  these  particular  studies,  a 
Rhesus  monkey  skull  was  placed  against  the  collimator  face  to  simulate 
gamma-ray  scattering  effects  in  bone. 

In  addition  to  the  single  photon  measurements,  coincidence  scans 

22 

in  both  air  and  water  were  performed,  without  collimation,  using  a  Na 
source. 

In  the  modelling  program,  the  crystal  was  always  assumed  to  be 
a  right  cylinder,  although  two  different  sizes  were  used.  One  crystal 
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was  taken  to  be  2.54  cm  in  diameter,  1.27  cm  thick,  of  the  type  used  in 
the  PFSS;  the  other  was  0.6  cm  diameter,  1.27  cm  thick  as  used  in  the 
miniature  system  for  cerebral  blood  flow  studies.  The  2.54  cm  diameter 
crystal  was  replaced  by  a  square  grid  with  points  spaced  at  2  mm  for  a 
total  of  129  points.  The  source  grid  points  were  set  at  a  spacing  of  5  mm 
but  the  source  array  was  made  much  larger,  with  40  x  40  grid  points. 

4.3  Experimental  Scans  and  Models 

The  experimental  and  simulated  scans  are  presented  in  Figures 
11-47  .  For  a  particular  combination  of  source  and  attenuating  medium, 

field  responses  for  the  cylindrical  straight,  converging  and  diverging 
collimators  are  presented  in  the  same  figure.  Field  responses  for  the 
rectangular  collimators  (straight  and  tapered)  for  the  same  source- 
attenuator  combinations  are  presented  immediately  following  the 
cylindrical  results  in  order  to  facilitate  comparisons.  The  simulated 
scans  for  particular  experimental  arrangements  are  presented  following 
the  experimental  data.  The  jaggedness  of  the  isocount  contours  in  the 
experimental  scans  is  due  to  random  fluctuations  in  source  activity. 

Scan  numbers  prefixed  with  an  M  are  those  produced  from  the  computer. 
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FIGURE  11 
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FIGURE  13 
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RECTANGULAR  0°  -  SUM  80  KEV  AIR  24.1cm  X  25cm 
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FIGURE  16 
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MODEL  RECTANGULAR  30°  80  KEY  AIR  24.1cm  X  25cm 
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FIGURE  19 
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FIGURE  20 
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FIGURE  22 
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FIGURE  24 


FIGURE  25 
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FIGURE  30 
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FIGURE  32 
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MODEL  RECTANGULAR  0°  511  KEY  AIR  19.6cm  X  25cm 
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RECTANGULAR  30°  511  KEY  AIR  19.6cm  X  25cm 


FIGURE  33 


63 


TTT 


5. 3 


MODEL 


MODEL 


FIGURE  34 


6 1\ 
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FIGURE  36 
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RECTANGULAR  45°  511  KEY  HO  18.5cm  X  25cm 
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FIGURE  40 
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MODEL  RECTANGULAR  30°  511  KEY  H20  18.5cm  X  25cm 
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FIGURE  43 
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FIGURE  44 
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MODEL  CYLINDRICAL  STRAIGHT  80  KEV  TISSUE+1  RIB  19.6cm  X  25cm 

TRANSVERSE  SECTION 


MODEL  CYLINDRICAL  STRAIGHT-SUM  80  KEV  TISSUE+1  RIB  19.6cm  X  25cm 

TRANSVERSE  SECTION 


FIGURE  46 
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MODEL  CYLINDRICAL  STRAIGHT  80  KEV  TISSUE+2  RIB  19.6cm  X  25cm 

SAGITTAL  SECTION 


SAGITTAL  SECTION 


FIGURE  47 
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4.4  Comparison  of  Collimators 

In  comparing  collimators,  it  is  apparent  that  the  cylindrical 
straight  collimator  has  response  lines  which  are  fairly  flat ,i.e. uniform 
with  depth.  The  enhancement  of  this  parameter  is  desired  in  the  PFSS,  but 
a  flatter  field  response  can  be  obtained  with  diverging  or  rectangular 
collimators . 

The  cylindrical  converging  collimator  has  the  highest  spatial  re¬ 
solution  of  all  the  large  collimators  tested,  as  can  be  seen  by  the  narrow¬ 
ness  of  the  response  lines.  However,  a  high  spatial  resolution  is  not 
required  in  the  PFSS  and  consequently  this  particular  collimator  is  not 
considered  with  respect  to  lung  function  studies.  Possible  applications  of 
this  type  of  collimator  could  be  in  regional  studies  of  smaller  organs  such 
as  the  kidney  or  thyroid.  Spatial  resolution  could  be  made  equal  to  that  of 
the  miniature  detection  system,  shown  in  Figure  43  but  the  sensitivity  would 
be  larger  for  the  converging  collimator,  due  to  the  larger  crystal. 

The  cylindrical  diverging  collimator  yielded  fairly  uniform  iso¬ 
response  lines  for  all  source-medium-detector  configurations.  While  this 
is  a  desired  characteristic  for  the  PFSS,  a  major  drawback  in  this  type  of 
collimator  is  the  circular  f ield-of-view.  With  this  f ield-of-view,  certain 
regions  of  the  lungs  are  not  "seen"  by  the  detectors;  consequently,  analysis 
is  complete. 

The  rectangular  straight  collimator  yielded  the  flattest  isoresponse 
pf  the  height  of  the  collimator  is  increased  from  the  height  of  iche 
measured  collimator (7 . 6  cm),  the  response  lines  would  become  flatter  still. 
The  f ield-of-view  of  this  collimator  is  given  by  equation  (B.15)  in 
Appendix  B  and  is  rectangular,  with  the  corners  slightly  rounded  off.  It 
is  this  type  of  collimator  that  can  view  the  complete  lung,  with  minimal 

overlap . 

In  the  higher  energy  scans,  using  Na  (Figures  27  and  29),  the 


‘i  • 


78 


f lareout  at  the  ends  is  due  to  collimator  penetration  due  to  improper 
shielding  of  the  crystal  at  the  back  of  the  collimator  and  should  be 
disregarded . 

For  tapered  rectangular  collimators,  only  diverging  collimators 
were  tested  as  a  large  f ield-of-view  in  the  PFSS  is  desired  and  this  cannot 
be  achieved  readily  with  converging  collimators.  Three  taper  angles  were 
measured:  30  ,45  ,  and  60  to  the  normal.  Although  the  diverging  rect¬ 
angular  collimators  give  slightly  better  flat  field  response,  particularly 
for  the  30  diverging  collimator,  the  field  lines  eventually  curve  at  the 
ends.  This  end  curvature  arises  since  the  collimator  no  longer  limits 
the  crystal,  and  a  distance  squared  dependence  results.  This  effect  be¬ 
comes  more  prominent  with  larger  taper  angle.  The  "f lareout"  at  the  ends, 
(evident  in  some  of  the  scans) ,  is  again  due  to  improper  crystal  shielding 
and  should  be  disregarded. 

With  variation  in  depth  the  "tightest"  and  most  uniform  iso¬ 
response  lines  were  obtained  for  the  positron  coincidence  scans,  which  in¬ 
volved  no  physical  collimation.  The  reasons  for  this  excellent  field-of- 
view  are  two-fold.  First,  positron  annihilation  results  in  two  gamma 
rays  travelling  in  directions  180°  apart.  In  order  for  a  coincidence 
event  to  be  detected,  the  annihilation  event  is  constrained  to  occur  in 
the  geometrical  space  of  the  cylinder  formed  between  the  two  crystals, 
which  therefore  makes  the  f ield-of-view  distance-independent.  Second, 
since  the  annihilation  radiation  is  of  a  reasonably  high  energy  (551  keV ) , 
the  attenuation  in  tissue  is  not  as  depth-dependent  as  for  80  keV  gamma 
rays.  However,  the  major  problem  in  using  positron  emitters  for  clinical 
work  is  that  most  have  short  half-lives  and  must  be  accelerator  produced. 
Also,  the  radiation  dose  to  the  patient  is  usually  larger  than  for 
’single  photon’  radioisotopes  since  large  amounts  of  activity  must  be 
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used  in  order  to  obtain  significant  numbers  of  coincidence  events. 

For  overall  comparisons,  in  comparing  source  energy,  it  is  apparent 
that  the  higher  energy  scans  (551  keV )  have  deeper  equi-valued  isocount 
response  lines  than  the  80  keV  scans  for  the  same  geometries.  This  of 
course  is  due  to  the  greater  attenuation  of  the  lower  energy  gamma  rays 
by  the  medium  they  traverse. 

In  comparing  source  media, a  higher  density  medium(water  as  compared 
to  air)  moved  the  isocount  curves  closer  to  the  crystal  face,  making  the 
depth  response  shallower.  This  again  is  due  to  the  higher  linear  attenuation 
coefficient  for  the  denser  medium  which  emphasizes  attenuation  as  a  function 
of  path  length. 

Water  or  tissue  also  scatters  gamma  rays  more  than  air  and  this  is 
evident  from  the  scans  made  in  water  since  the  isocount  curves  are  much 
broader  than  the  corresponding  scans  made  in  air.  The  broadening  results 
from  scattering  which  allows  the  crystal  to  detect  gamma  rays  whose  path 
direction  would  not  normally  have  been  through  the  crystal. 

The  high  energy  scans  are  also  slightly  broader  than  their  low 
energy  counterparts.  This  effect  results  from  edge  penetration  of  the  col¬ 
limator.  Edge  penetration  occurs  for  gamma  rays  from  just  outside  the  geo¬ 
metrical  f ield-of-view  of  the  collimator,  and  is  more  prominent  in  a  brass 
collimator  than  in  a  lead  collimator. 

All  the  collimators  yielded  roughly  the  same  depth  response  for 
any  particular  energy, with  the  50%  isoresponse  line  lying  between  8  and  10  cm 
from  the  collimator  face  for  80  keV  gamma  rays  in  air  and  5  and  7  cm  in  water. 

The  major  differences  between  the  several  scans  is  the  field— of 
view,  which  is  reflected  in  the  broadness  of  the  isocount  curves.  The  field— 
of— view,  which  was  discussed  in  Chapter  2, is  dependent  mainly  on  collimator 
geometry.  Bearing  in  mind  that  neither  the  smallest  f ield-of-view  nor  the 
largest  is  necessarily  the  best  f ield-of-view  for  any  particular  application, 
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the  cylindrical  converging  collimator  gave  the  narrowest  f ield-of-view 
among  the  collimators,  with  the  cylindrical  diverging  and  rectangular 
diverging  collimators  having  the  broadest  values.  The  isocount  contours 
for  any  one  scan  also  reflect  the  widening  of  the  f ield-of-view  with 
increasing  distance  from  the  collimator  face.  This  can  be  seen  by  com¬ 
paring  a  50%  response  contour  with  a  75%  contour  and  noting  that  the  50% 
contour,  which  lies  deeper,  is  also  wider  than  the  75%  response  line. 

The  modelling  scans,  prefixed  with  an  M  in  the  diagrams,  simulate 

the  experimental  results  quite  well.  The  major  differences  seen  arise 

from  gamma  ray  scattering  in  the  experimental  studies  which  was  not 

accounted  for  in  the  simulations.  However,  these  effects  were  not  large 

enough  to  create  significant  alterations  in  the  isocount  profiles.  Of 

particular  interest  are  scans  MR1  and  MR2  which  simulate  lung  geometry  that 

could  not  be  measured  experimentally.  ScanMRlhas  one  rib,  1.25  cm 

diameter,  across  the  middle  of  the  detector  field  of  view;  while  scan  MR2 

assumes  two  ribs  lie  in  the  field  of  view,  both  1.25  cm  in  diameter. 

The  scans  show  the  amount  of  radioactivity  that  is  incident  on  the 

crystal  as  a  function  of  crystal  position  (vector  defined  previously) . 

This  vector,  I  A...  which  is  the  total  amount  of  activity  that  reaches 
i  iJ 

detector  point  j ,  is  not  a  useful  parameter  in  a  single  collimator 

counting  system,  since  the  photomultiplier  attached  to  the  crystal  cannot 

differentiate  spatial  events  inside  the  crystal.  However,  it  is  useful 

in  an  imaging  system  employing  a  multihole  collimator.  %  A.,  would  then 

i 

be  the  amount  of  activity  that  is  imaged  in  position  j  on  the  crystal 
face.  By  comparing  simulated  images  from  known  source  distributions  to 
experimental  images,  reasonable  estimations  of  the  unknown  source 
distributions  can  be  obtained. 
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CHAPTER  FIVE 


SUMMARY  AND  CONCLUSION 

It  was  the  purpose  of  this  study  to  investigate  gamma-ray 

collimator  characteristics  with  specific  aims  of  defining  and  of 

simplifying  collimator  design  procedures  for  various  clinical  measurement 

situations.  The  particular  clinical  situation  considered  in  this  study 

is  the  measurement  of  regional  pulmonary  function  using  an  inert,  diffusible 

133 

radioactive  gas  (  Xe) .  Specific  goals  of  the  optimization  process  are: 

(i)  to  provide  a  well  defined,  non-overlapping  f ield-of-view 
between  detectors, 

(ii)  to  obtain  maximum  sensitivity  for  this  field-of  view  ,  and 

(iii)  to  produce  a  uniform,  or  iso-sensitive  counting  field 
at  a  distance  from  the  collimator  face. 

As  discussed  in  Chapter  4,  all  collimators  tested  yield  roughly 

the  same  sensitivity  with  depth  in  any  one  medium.  In  comparing  the 

response  of  the  rectangular  straight  collimator  to  the  responses  of  the 

133 

cylindrical  straight  and  diverging  collimators  for  Xe  in  water,  it 
can  be  seen  that  the  straight  rectangular  collimator  provides  the 
broadest  field  response.  Also,  since  the  f ield-of-view  of  a  rectangular 
collimator  is  rectangular,  this  type  of  collimator  on  the  PFSS  will  view 
the  complete  lung  and  not  just  cylindrical  sections  as  in  the  current 
system.  Field-of-view  overlap  can  also  be  more  easily  minimized  using 
rectangular  collimators. 

As  an  optimum  condition  in  the  PFSS,  a  field-of-view  width  of 
8  cm  is  desired  at  a  distance  13  cm  from  the  collimator  face.  The 
f ield— o f — view  is  also  required  to  include  the  chest  from  the  sternum  to 
the  side  of  the  rib  cage,  a  distance  of  approximately  18  cm.  In 
considering  the  results  of  this  study,  the  ideal  collimator  for 
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optimizing  these  several  parameters  would  be  a  tapered  rectangular 
collimator.  The  crystal  end  dimensions  should  be  square,  having  sides 
equal  to  the  diameter  of  the  crystal.  The  face  end  dimensions  should  be 
a  width  equal  to  the  crystal  diameter;  the  face  end  length  7.6  cm.  The 
overall  length  of  this  collimator  should  be  12.7  cm  and  only  the  side 
distal  to  the  sternum  should  be  tapered.  This  collimator  will  have 
good  sensitivity  and  collimators  of  this  type  used  in  the  PFSS  will  view 
the  complete  lung  with  minimal  overlap  between  detectors. 

As  an  additional  concluding  remark,  the  modelling  programs 
used  in  this  work  may  be  utilized  in  testing  characteristics  of  potential 
collimators.  Using  the  programs,  since  collimator  specifications  can  be 
varied  easily  and  quickly,  corresponding  field  responses  for  different 
configurations  can  be  determined  analytically  without  the  need  of 
fabrication  and  testing. 
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APPENDIX  A 


The  Physics  of  the  Lung 

A  continuous  and  adequate  supply  of  oxygen  is  a  fundamental 
requirement  for  the  metabolism  of  all  mammalian  cells.  In  addition, 
a  means  is  needed  to  remove  unwanted  metabolic  by-products  of  these 
cells  -  principally  carbon  dioxide.  The  organ  which  meets  these 
requirements  is  the  lung. 

Until  recently  the  lung  could  only  be  studied  "in  vivo" 
by  making  measurements  at  the  mouth,  or  by  bronchospirometry  in  which 
the  patient  must  be  anesthetized  in  order  to  study  individual  lung 
lobes.  Nuclear  medicine  has  been  important  in  furthering  our 
knowledge  of  regional  lung  function  at  the  lower  anatomical  scale 
and  this  Appendix  will  discuss  some  of  the  physical  processes  involved. 

A.l  Mechanics  of  Respiration 
A. 1.1  Expansion  and  Contraction 

The  lung  itself  is  passive  and  has  no  musculature  except 
for  some  smooth  muscle  in  the  airways  and  blood  vessels.  Its  main 
purpose  is  to  bring  blood  and  air  very  close  together  to  enable  gas 
exchange.  This  it  does  via  the  alveolar  membrane,  which  is  of  the 
order  of  one  micron  in  thickness  and  has  a  total  surface  area  for 
both  lungs  in  the  normal  male  adult  of  about  70  square  meters. 

The  act  of  breathing  depends  upon  the  fact  that  the 
thoracic  cavity  is  a  closed  compartment  and  the  lung  contained 
in  it  is  open  to  the  exterior  via  the  trachea.  (The  basic  anatomy  of 
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the  lung  is  shown  in  Figure  48).  Thus,  when  the  total  volume  of  the 
cavity  increases,  negative  pressure  (with  respect  to  atmospheric) 
in  the  cavity  causes  air  to  be  drawn  inward  through  the  trachea; 
when  the  volume  decreases,  the  pleural  pressure  goes  toward  atmospheric 
and  lung  volume  decreases  (Figure  49)  .  The  lung  is  therefore  able  to 
breathe  by  periodic  expansions  and  contractions  of  the  thoracic 
cavity  brought  about  by  intermittent  contractions  of  the  respiratory 
muscles  and  passive  recoil  of  the  elastic  lungs. 

The  enlargement  of  the  thorax  during  quiet  breathing  is 
brought  about  primarily  by  diaphragmatic  contraction  and  by  movements 
of  the  ribs,  sternum,  and  vertebrae.  As  the  intercostal  muscles  contract 
the  ribs  move  upwards  bringing  the  plane  of  the  ribs  nearer  to  a  perpen¬ 
dicular  with  the  vertebral  column,  thus  increasing  the  thoracic  volume. 

The  most  important  inspiratory  muscle  is  the  diaphragm,  which 
serves  as  a  "piston"  in  a  thoracic  "cylinder".  It  is  dome  shaped, 
due  partly  to  negative  intrathoracic  pressure  and  partly  to  positive 
intra-abdominal  pressure  which  produce  approximately  10  Newtons  of 
force  acting  over  its  surface.  Its  direction  of  movement  is  vertical, 
with  a  surface  area  of  approximately  250  square  centimeters.  A  descent 
therefore  of  1  centimeter  increases  thoracic  volume  by  about  250  cc. 

As  this  is  an  increase  in  volume  in  the  cavity  at  a  constant  temperature, 
Boyle’s  law  dictates  that  the  pressure  must  decrease,  which  forces  the 
lungs  to  expand  thus  drawing  in  air  from  the  airways.  This  process  is 
reversed  as  the  diaphragm  makes  its  regulated  ascent.  The  excursion 
of  the  dome  averages  1.25  cm  during  normal  breathing  (  39  )  . 
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(lb)  DESCENT  OF  DIAPHRAGM  AT  THE  SAME 
TIME  FURTHER  PROMOTES  CAVITY 
EXPANSION 


FIGURE  49 
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A. 1.2  Respiratory  Pressure  -  The  Lung  Pump 

As  the  ribs  expand  at  the  start  of  an  inspiration,  the 
intrathoracic  pressure  is  reduced  to  about  -5mmHg  subatmospheric . 

This  negative  pressure  in  turn  affects  the  intrapulmonary  pressure 
which  is  the  pressure  existing  in  the  air  passages  and  the  alveoli 
within  the  lungs.  During  normal  breathing,  intrapulmonary  pressure 
within  the  alveoli  ranges  from  about  -3mmHg  subatmospheric  during 
inspiration  to  atmospheric  at  the  end  of  inspiration,  to  +3mmHg  at 
the  beginning  of  expiration  and  finally  back  to  atmospheric  pressure 
at  the  end  of  expiration,  thus  completing  one  cycle  (Figure  50) • 

Normal  intrapulmonic  pressure  fluctuations  can  be  greatly 
exceeded  by  vigorous  contraction  of  the  inspiratory  and  expiratory 
muscles  wThile  breath  holding  with  the  glottis  closed.  Values  as  high 
as  -90mmHg  and  +110mmHg  can  be  produced, but  to  prevent  acute 
circulatory  depression  the  safe  range  of  maxima  are  about  -20  to 
+30mmHg .  The  reason  for  this  circulatory  depression  is  that  the 
veins  are  thin— walled  collapsible  tubes,  and  if  the  outside  pressure 
is  greater  than  the  inside  driving  pressure,  the  tube  will  collapse 
and  flow  will  stop  (  40  ) . 

A. 1.3  Compliance  and  Work 

Ventilation  of  the  alveoli  requires  that  the  respiratory  muscles 
overcome  several  types  of  resistance.  The  work  done  by  these  muscles  can  be 
estimated  by  multiplying  the  mean  force  of  muscular  contraction  by  the  mean 
distance  of  muscle  movement.  These  measurements  can  be  made  using  miniature 
isometric  strain  gauges  and  microelectrodes  implanted  in  the  muscles. 
However,  for  the  lung,  a  much  easier  and  accepted  way  is  to  express  the  work 
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INTRAPULMONIC  PRESSURE 


INSPIRATION  |  EXPIRATION 

(FROM  PERKINS  1961) 

(51) 


INTRA-ALVEOLAR  PRESSURE  (cm  HO) 


COMPLIANCE  DIAGRAM  OF  THE  LUNGS  ALONE. 

THE  SHADED  AREA  REPRESENTS  JP  dv,  THE 
WORK  REQUIRED  TO  OVERCOME  NONELASTIC 
TISSUE  RESISTANCE.  THE  DOTTED  LINE  REP¬ 
RESENTS  THE  COMPLIANCE. 

(FROM  GUYTON  1971) 


THE  OUTER  LOOP  IS  THE  PRESSURE- 
VOLUME  RECORDED  CONTINUOUSLY  DUR¬ 
ING  A  CYCLE  OF  RESPIRATION  WHILE 
THE  INNER  LOOP  IS  THE  ABOVE  STATIC 
COMPLIANCE  DIAGRAM.  THE  SHADED 
AREA  REPRESENTS  THE  INCREASED  WORK 
NEEDED  TO  OVERCOME  AIRWAYS  RESIS- 
TANCE 

*  (FROM  GUYTON  1971) 


(52) 


FIGURE  50 
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done  in  terms  of  pressure  change  (AP)  times  volume  change  (AV) ,  both 
of  which  can  be  measured  using  either  an  intraesophageal  balloon  or  a 
less  traumatic  body  plethysmograph  (  41  ) . 

A  graph  of  lung  volume  change  verses  intrapleural  pressure 
change  gives  a  compliance  diagram  (Figure  50) .  Compliance  is  the 
ability  of  the  lungs  to  be  expanded.  It  is  defined  as  the  ratio  of 
volume  change  to  pressure  change  and  can  be  either  a  static  or 
dynamic  phenomenon.  The  lungs  alone,  when  removed  from  the  chest, are 
roughly  twice  as  distensible  as  the  lungs  and  thorax  together  because 
the  thoracic  cage  itself  must  also  be  stretched  in  normal  breathing, 
which  is  an  energy  drain  on  the  respiratory  muscles. 

The  overall  mechanical  energy  expended  by  the  respiratory 
muscles  is  used  up  in  three  main  ways:  1)  Elastic  work  is  the  work 
required  to  stretch  the  elastic  tissues  of  the  lungs  and  thorax  and 
can  be  calculated  from  compliance  diagrams,  2)  inertial  work  is  related 
to  the  force  needed  to  set  tissues  such  as  the  liver,  diaphragm  and 
chest  wall  in  motion.  (However,  this  is  significant  only  at  high  rates 
of  ventilation  and  is  negligible  in  normal  breathing) , and  3)  work  done 
to  overcome  airway  resistance  is  that  energy  expended  in  moving  the 
respiratory  gases  through  the  airways.  As  can  be  seen  from  Figure  50 
airway  resistance  accounts  for  a  major  portion  of  the  available  energy . 
It  will  be  discussed  in  greater  detail  in  the  following  section. 

The  total  work,  therefore,  can  be  determined  and  has  been 
found  to  range  in  normal  subjects  between  0. 2-0.4  Joules/breath. 

The  mechanical  efficiency  of  this  breathing  is  defined  as  the 
respiratory  work  done  per  unit  of  oxygen  required  for  the  work,  or 

work  per  liter 
O2  cost  per  liter. 
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This  efficiency  is  about  8%  in  normals  and  depends  upon  the  frequency 
and  depth  of  breathing. 

A. 1.4  Surface  Tension  and  Elastic  Recoil 

The  elasticity  of  the  lung  is  not  due  solely  to  its  elastic 
fibers.  Of  equal  importance  for  elastic  recoil  is  the  surfactant 
that  covers  all  lung  alveoli.  In  1929  von  Neergaard  inflated  the  lungs 
of  anaesthetized  cats  and  measured  deflation  pressures  (  42  )  .  He 
used  both  gas  and  physiologic  liquids  to  distend  the  lungs  and  found 
that  the  pressure  required  to  maintain  any  given  lung  volume  with 
fluid,  under  static  conditions,  was  less  than  half  that  required  to 
maintain  the  same  volume  of  an  air  filled  lung.  Von  Neergaard  deduced 
that  there  must  be  a  liquid-air  interface  on  the  alveoli  which  tends 
to  retract  the  alveoli,  thus  contributing  to  the  lung  recoil.  Filling 
the  lung  with  liquid  would  eliminate  the  fluid-gas  interface  and  its 
recoil  force. 

Using  Laplace’s  Law, 


and  approximating  the  alveolus  as  a  sphere,  the  pressure  (P)  is 
related  to  the  surface  tension  (y)  and  the  radius  (r)  of  the  bubble  . 
Assuming  that  the  alveolar  film  exerts  a  surface  tension  of  50  dynes/cm 
which  is  equal  to  that  of  plasma,  it  has  been  calculated  that  the  recoil 
pressure  of  the  lung  is  roughly  2  x  104  dynes/cm2  or  about  20  cm  H^O  for 
a  fully  inflated  lung. 

This  is  verified  experimentally  but  it  is  almost  an  order  of 
magnitude  too  large  at  normal  or  low  lung  volumes.  The  implication  is 
that  the  surface  tension  of  the  alveoli  must  be  5  -  10  dynes/cm  instead 
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of  50  dynes/cm.  This,  however,  could  only  occur  if  the  alveolar 
membrane  was  lined  with  something  that  reduced  the  surface  tension. 
This  indeed  is  the  case  as  the  walls  of  the  alveoli  secrete  a  material 


with  a  high  dipalmitoyl  lecithin  concentration  which  reduces  the  surface 
tension  by  a  factor  as  great  as  15.  The  surfactant’s  ability  to  reduce 
surface  tension  depends  upon  its  concentration.  Therefore,  when  the 
alveolus  is  deflated,  the  concentration  (per  unit  surface  area)  is  high 
and  the  surface  tension  is  very  low,  so  that  the  alveolus  is  expanded 
without  difficulty.  However,  as  it  expands,  the  concentration  of  the 
surfactant  decreases  and  the  surface  tension  increases,  until  a  point 
of  equilibrium  is  reached  at  maximum  expansion.  Upon  exhaling  the 
increased  surface  tension  helps  to  collapse  the  alveolus  and  expel  the 
air  from  it. 

Consider  for  example  that  during  inhalation  the  radius  of  an 

2> 

alveolus  expands  from  about  0.5  x  10  2  to  1.0  x  10  cm  (  43  ) .  Assuming 
a  surface  tension  of  50  dynes/cm,  the  pressure  differences  required  to 
inflate  an  alveolus  would  be: 


(A. 2) 

=  2  x  104  dynes/cm2  =  15  mmlig 

This  means  that  the  gauge  pressure  PQ  outside  the  alveolus 

would  have  to  be  15  mmHg  less  than  the  pressure  =  -3  mmHg  inside  it. 
Thus  PQ  would  be  -18  mmHg.  However,  the  intrapleural  pressure  is  only 
about  -4  mmHg,  so  that  the  pressure  difference  is  only  about  1  mmHg,  or 
15  times  less  than  would  be  required.  Pulmonary  surfactant  overcomes 


this  difficulty. 
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Besides  reducing  muscular  effort  in  inflating  the  lungs, 
pulmonary  surfactant  helps  to  stabilize  the  alveolar  gas  bubbles  and 
keep  them  from  collapsing.  Alveoli  vary  in  size  throughout  the  lungs, 
some  having  radii  3-4  times  that  of  others.  Laplace's  Law  shows  that 
if  the  surface  tension  was  equal  for  bubbles  of  differing  radii,  the 
smaller  bubbles  would  collapse  and  the  larger  ones  would  overexpand. 

The  lung  therefore,  would  consist  of  one  hyperinf lated  alveolus,  with 
all  other  alveoli  being  collapsed.  That  this  is  not  the  case  is  due 
solely  to  the  surfactant.  A  collapsed  alveolus  without  a  coating  of 
surfactant  would  behave  like  a  bubble  being  formed  on  the  end  of  a 
straw  (Figure  51) .  The  alveolar  pressure  would  have  to  increase  to 
relatively  high  levels  before  the  alveolus  would  "pop"  open.  Thereafter, 
less  pressure  would  be  required  to  expand  the  alveolus  to  a  higher 
volume . 

Table  A-l  shows  the  elastic  properties  of  the  chest  wall  and 
lungs.  It  is  apparent  from  this  table  that  the  lungs  seek  a  lower 
volume  even  at  residual  volume  (20%  of  lung  capacity)  while  the  chest 
wall  has  an  expanding  tendency  until  the  lungs  are  inflated  to 
approximately  50%  vital  capacity.  Lung  volume  greater  than  50%  vital 
capacity  results  in  the  chest  wall  recoiling  inward  rather  than  outward 
as  it  does  below  50%  vital  capacity. 

A. 2  Flow  Dynamics 

A. 2.1  Volumes  and  Capacities 

The  basic  nomenclature  of  pulmonary  function  is  given  in 
Figure  52  .  Physiologically,  the  main  bronchi,  bronchioles,  and  terminal 
bronchioles,  which  do  not  play  a  part  in  actual  gas  exchange,  can  be 
classified  as  airways.  The  respiratory  bronchioles  and  the  alveoli 
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PRESSURE-VOLUME  CURVE  OF  A  BUBBLE  BLOWN 
AT  THE  END  OF  A  GLASS  TUBE. 

(FROM  COMROE  1974)  (56) 

FIGURE  51 


TABLE  A-l  (FROM  KNOWLES,  HONG,  AND  RAHN;  1959) 

(53) 


%  of  Vital 
Capacity 

Recoil  (Relaxation) 
Pressure  of  Lungs  + 
Thoracic  Cage 
(Torr) 

Recoil  Pressure 
of 

Lungs  Alone 
(Torr) 

Recoil  Pressure 

of 

Thoracic  Cage  Alone 
(Torr) 

20 

-8 

1.8 

-9.8 

40 

4 

5.0 

-1.0 

60 

14 

8.9 

4.1 

80 

21 

13.3 

7.7 

100 

30 

21.0 

9.0 

LUNG  VOLUME  (L) 
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TIME 


FIGURE  52 
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on  the  other  hand,  may  be  considered  the  basic  lung  units  involved 
with  oxygen  transfer  to  the  blood  cells.  There  are  therefore  two  main 
areas  of  the  lung,  one  involved  with  gas  transport,  and  the  other  with 
gas  exchange.  Making  the  assumption  that  all  units  are  basically  the 
same,  the  whole  organ  can  be  depicted  as  shown  in  Figure  53.  Typical 
values  are  included. 

A. 2. 2  Rates  and  Flows 

As  outlined  previously,  the  total  force  exerted  by  respiratory 
muscles  must  overcome  several  types  of  resistance  in  ventilating  the 
pulmonary  alveoli.  The  main  one  is  airways  resistance. 

In  fluid  flow. 


(A.  3) 


Driving  pressure 

Flow 


R  = 


where  R  is  the  resistance  to  flow  in  the  tube  or  tubes.  In  the  lung 
however,  there  are  two  types  of  flow,  laminar  and  turbulent,  and  each 
requires  a  different  driving  pressure.  The  first  step  therefore,  is 
to  determine  what  type  of  flow  exists  at  various  levels  in  the  system. 
This  can  be  done  using  the  Reynolds’  number. 


(A.  4) 


The  flow  of  a  fluid  of  density  P  and  viscosity  n  at  an  average  speed  of 
v  through  a  tube  of  radius  r  is  laminar  as  long  as  the  Reynolds'  number 
is  less  than  1000.  Note  that  both  the  viscosity  and  the  density  of  the 


gas  in  question  are  important.  Table  A-2  lists  both  of  these  parameters 
for  gases  involved  in  respiration.  Table  A- 3  gives  the  calculated  Reynolds' 
number  at  various  sites  in  the  respiratory  tract.  The  velocity  of  air 


flow  in  the  alveoli  is  gene 


rally  thought  to  be  small  enough  to  be  negligible. 


f 


* 
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ANATOMIC  DEAD  SPACE 


0.15  L 


BLOOD  FLOW 


70  mL 


5.0  L/MIN 

(0) 


NUMBERS  ARE  TYPICAL  VALUES  ONLY  AND  THERE 
IS  CONSIDERABLE  VARIATION.  NOTE  THAT  NOR¬ 
MAL  RATIO  OF  VENTILATION  TO  BLOOD  FLOW  IS 


FIGURE  53 
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TABLE  A- 2  (TAKEN  FROM  KAO,  F.F.  1972)  (45) 


GAS 

n  6 
(POISE  x  10  ) 

P  -3 

(g/ML  x  10  ) 

at  Atmospheric 

OXYGEN 

214 

126 

NITROGEN 

184 

110 

C02 

156 

173 

WATER  VAPOR 

103 

70.8 

HELIUM 

203 

15.6 

RESPIRED  AIR 

180 

112 

TABLE  A- 3  REYNOLD’S  NUMBER  (ADAPTED  FROM  MEAD,  1961)  (54) 


LOCATION 

RADIUS 

(MM) 

V  =  6  L/MIN 

V  =  60  L/MIN 

V  =  200  L/MIN 

NASAL  CANAL 

2.5 

200 

2000 

6000 

PHARNYX 

6 

400 

4000 

12000 

GLOTTIS 

4 

800 

8000 

24000 

TRACHEA 

10.5 

625 

6250 

18500 

BRONCHI 

8.5 

455 

4550 

13650 

4.5 

350 

3500 

10500 

3 

285 

2850 

8550 

2 

95 

950 

2850 

0.5 

17.5 

175 

525 
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Where  the  gas  flow  is  laminar,  Poiseuille's  Law  states  that: 


AP  =  V 


(A.  5) 


which  relates  the  volume  flow  per  unit  time  (V)  to  the  pressure  gradient 
AP  in  a  tube  of  length  £  and  radius  r,  r)  being  the  viscosity.  It  is  seen 


Trr 


where  is  the  resistance  to  laminar  flow. 
Where  the  flow  is  turbulent, 


AP  =  V2 

Tr“r 


R_ 


(A.  6) 


(A.  7) 


p  being  the  density  of  the  gas  and  k  a  constant  relating  to  frictional 
effects  between  the  gas  and  the  wall.  Here, 


k££ 

kt  -  2  2 

tt  r 


(A.  8) 


where  R^  is  the  resistance  to  turbulent  flow.  Note  that  R^  depends 
mainly  on  the  gas  viscosity,  while  R^  depends  primarily  on  gas  density. 

In  smooth,  straight  tubes,  turbulent  flow  occurs  only  at  high 
velocities  which  are  present  in  the  main  bronchi  and  trachea.  The  flow 
rates  in  the  smaller  anatomical  levels  are  very  low  because  the  total 
airflow  is  divided  among  hundreds  of  thousands  of  tubes  and  the  total 
cross-sectional  area  of  these  tubes  is  very  large  leading  to  very  low 
flow  velocity.  However,  eddy  formation  may  occur  at  each  branching  of 
the  tracheobronchial  tree  and  the  pressure  required  for  eddy  flow  is 
approximately  the  same  as  for  turbulent  flow.  Turbulence  or  eddy 
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formation  at  low  rates  is  particularly  apt  to  occur  when  there  are 
irregularities  in  the  tubes.  In  the  respiratory  system  airflow  is 
both  laminar  and  turbulent  with  eddy  formation.  Rohrer’s  equation 
combines  these  influences  and  the  flow  pressure  relationship  for  this 
"impure"  physiologic  state  is 

AP  =  V  +  R^,  V2  (A. 9) 

Figure  54  shows  an  experimental  determination  of  this  equation.  The  slope 
of  the  curve  at  any  point  gives  the  resistance  to  flow  at  that  particular 
pressure  and  flow  rate. 

A. 3  The  Lung  and  Gas  Exchange 
A. 3.1  The  Pressures 

Only  a  few  of  the  cells  of  human  tissues  could  survive  extended 
direct  exposure  to  the  partial  pressures  of  gases  present  in  the  earth's 
atmosphere  (  44  )  .  Oxygen  at  sea  level  pressures  is  distinctly  toxic, 
and  so  the  pulmonary  structure  evolved  in  man  has  had  to  provide  a 
portable  gaseous  environment.  This  is  the  environment  found  in  the 
alveoli  and  is  considerably  different  from  that  of  the  ambient 
atmosphere . 

In  respiratory  gas  exchange,  the  separate  behaviours  of  the 
various  gases  must  be  considered,  since  their  physiologic  roles  are 
also  different  from  each  other.  This  state  goes  well  with  Dalton’s  law, 
which  enables  the  partial  pressures  of  each  gas  to  be  treated  distinctly. 
The  partial  pressure  of  a  gas  is  of  extreme  importance  both  in  the 
physical  dynamics  of  pulmonary  -  blood  -  tissue  gas  exchange  and  in  the 
participation  of  the  gets  in  the  chemical  processes  of  metabolism. 


PRESSURE  GRADIENT  (cm  HO) 
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THE  FLOW-PRESSURE  CURVE  OF  THE  RESPIRA¬ 
TORY  SYSTEM.  THE  EXPERIMENTAL  DATA  POINTS 
LIE  CLOSEST  TO  THE  CURVE  WITH  EQUATION: 

P (cm  H  O)  =  0.0417  V  +  0.0009  V2 

(FROM  MEAD  1960)  (55) 


FIGURE  54 
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The  two  main  gases  involved  in  respiration  are  oxygen  and  carbon 
dioxide.  Oxygen  supplies  energy  needed  by  the  body's  cells  for  their 
metabolic  processes  and  carbon  dioxide  is  the  by-product  of  tissue 
metabolism.  In  addition,  there  is  the  "inert"  fraction  of  the  atmos¬ 
phere,  consisting  primarily  of  nitrogen,  but  also  including  trace 
amounts  of  argon,  xenon,  krypton,  hydrogen  and  helium  (Table  A-4)  . 

Movement  of  the  individual  gas  molecules  by  the  time  they 
reach  the  alveoli  is  largely  Brownian,  but  because  of  pressure  gradients, 
mass-movement  occurs  by  diffusion  from  the  higher  partial  pressure  of 
gas  to  the  lower  partial  pressure.  At  this  level,  each  gas  behaves 
to  a  good  approximation  according  to  the  ideal  gas  law  and  can  be 
treated  completely  by  application  of  the  kinetic  hypothesis.  It  is  in 
this  way  that  the  lung  functions ,  bringing  blood  and  air  very  close 
together,  separated  only  by  the  alveolar  membrane.  With  an  alveolar 
P02  of  100  mmHg  and  a  venous  blood  P02  of  40  mmilg,  a  driving  pressure 
of  60  mmHg  forces  oxygen  into  the  blood  where  most  of  it  combines  with 
hemoglobin  (Table  A-5 ,  Figure  55).  Figure  55  shows  that  very  little  time 
is  required  for  equilibration  of  alveolar  with  pulmonary  capillary  P02 . 
This  equilibration  occurs  even  during  exercise  when  the  pulmonary 
capillary  transit  time  is  decreased.  Figure  55 shows  how  this  rapid 
equilibration  can  be  upset  by  alteration  in  the  diffusion  properties 
of  the  alveolo capillary  membrane.  Increased  thickness  of  this  membrane 
increases  the  time  required  for  equilibration.  If  the  defect  is  severe 
enough  the  systemic  arterial  blood  will  be  incompletely  oxygenated. 
Conversely,  a  venous  blood  PCC>2  of  roughly  46  mmHg  and  alveolar  PC02 
of  40  mmHg  results  in  a  PC02  gradient  favorable  for  elimination  of 

carbon  dioxide  from  the  body . 
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TABLE  A-4  GASEOUS  COMPOSITION  OF  THE  ATMOSPHERE 
(FROM  MOUNTCASTLE,  1968)  (44) 


ELEMENT 

%  IN  DRY 

AIR 

NITROGEN 

78.09 

OXYGEN 

20.94 

ARGON 

0.93 

CARBON  DIOXIDE 

0.03 

NEON 

0.002 

HELIUM 

0 . 000!) 

KRYPTON 

0.00001 

HYDROGEN 

0.00005 

XENON 

0.000008 

TABLE  A-5  PARTIAL  PRESSURES  OF  RESPIRATORY  GASES  (mmHg) 

(FROM  MOUNTCASTLE,  1968)  (44) 


GAS 

AMBIENT 

AIR 

TRACHEAL 

AIR 

EXPIRED 

AIR 

ALVEOLAR 

AIR 

ARTERIAL 

BLOOD 

MIXED 

VENOUS 

BLOOD 

TISSUES 

0„ 

158.0 

149.0 

116 

100 

95 

40 

<40 

2 

C0o 

0.3 

0.3 

32 

40 

40 

46 

>46 

2 

h2o 

5.7 

47.0 

47 

47 

47 

47 

47 

N2 

596.0 

563.7 

565 

573 

573 

573 

573 

totals 

760 

760 

760 

760 

755 

706 

<706 
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ALVEOLOCAPILLARY  BLOCK  AND  CURVE  C  A  SEVERE  CASE  OF  DIFFUSION  IMPAIRMENT . EX¬ 
ERCISE  SHORTENS  THE  TIME  BY  AS  MUCH  AS  2/3  BUT  NOTE  THAT  AT  0.25  SEC  THE  pO 
OF  A  IS  ALREADY  AS  HIGH  AS  90  mrnHp,.  (ADAPTED  FFOM  COMROE  1974)  (56)2 
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FIGURE  55 
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A. 3.2  Diffusion  and  the  Alveolar  Membrane 

Diffusion  is  a  very  important:  physiologic  process  in  the  body 
because,  in  the  end,  it  alone  governs  the  transport  of  oxygen  from  the 
lungs  and  thereby  to  the  mitochondria  of  cells.  The  diffusion  of  0^ 
and  CO^  that  is  essential  to  pulmonary  gas  exchange  takes  place  between 
gas  and  tissue  -  that  is  between  gaseous  and  liquid  phases.  Thus  the 
solubility  of  the  gas  in  the  liquid  now  becomes  important.  Henry's  Law 
states  that  the  volume  of  a  fairly  insoluble  gas  that  dissolves  in  a 
liquid  at  a  given  temperature  is  almost  directly  proportional  to  the 
partial  pressure  of  that  gas.  Since  the  diffusion  rate  of  a  gas  within 
a  liquid  depends  on  its  pressure  gradient,  the  greater  the  solubility 
of  a  gas,  the  more  rapid  will  be  its  rate  of  diffusion.  The  relative 
solubilities  of  CO^  and  0^  in  water  are  nearly  25  to  1  which  explains 
why  CO^  diffuses  far  more  rapidly  between  alveolii  and  capillary  blood, 
even  though  the  pressure  gradient  for  0^  is  much  higher  than  CO^. 
However,  the  ability  of  hemoglobin  to  combine  chemically  with  0 ^  and 
release  CO^  is  a  factor  of  greater  importance  in  the  dif fusibility  of 
these  gases ,  but  detailed  consideration  of  these  phenomena  is  not 
within  the  scope  of  this  work. 

Alveolar  walls  can  be  as  thin  as  0.2  micron,  and  calculations 
have  shown  that  diffusion  of  gas  in  a  normal  alveolus  is  802  complete 
in  0.002  sec,  if  the  diffusion  distance  is  0.5  mm.  However,  there  are 
several  diseases  which  can  cause  the  alveolocapillary  membranes  to 
thicken  dramatically  (alveolocapillary  block).  These  give  good  evidence 
to  the  importance  of  the  length  of  the  diffusion  path  as  discussed  above 

It  is  also  obvious  that  the  area  of  the  diffusing  surface  is 
important,  for  the  body  must  not  only  maintain  a  highly  inflexible  rate 
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of  oxygen  transport,  but  a  stringently  governed  volume  of  oxygen  uptake 
as  well.  Certain  diseases  can  destroy  alveolar  walls,  decreasing  the 
area  from  its  normal  average  of  70  square  meters. 

In  addition  to  these  problems,  blood  -  gas  mixing  is  also 
important.  There  are  capillaries  in  the  lung  that  do  not  perfuse 
ventilated  alveoli,  and  the  deoxygenated  blood  in  them  is  sent  through 
the  body  a  second  time.  This  lowers  the  percentage  of  oxygenated 
hemoglobin,  which,  in  turn,  decreases  the  oxygen  supply  to  the 
metabolizing  tissues. 

In  summary,  mechanics,  flow  dynamics,  diffusion  and  physiologic 
and  anatomic  dead  space  all  contribute  to  the  broader  category  of 
pulmonary  ventilation.  The  matching  of  ventilation  with  blood  flow 
determines,  to  a  large  extent,  the  degree  to  which  the  peripheral 
tissues  will  be  oxygenated.  By  using  radioactive  materials  and  appro¬ 
priate  methodologies,  scintillation  detectors  positioned  over  the 
chest  can  provide  information  about  the  distributions  of  ventilation. 

The  detector  is  therefore  very  important  in  the  assessment  of  general 
lung  performance,  and  any  improvements  in  it  may  allow  more  accurate 
diagnosis  of  lung  impairments. 


APPENDIX  B 


Collimator  Field-of-View  Derivations 

For  single  hole  collimators,  the  index  of  spatial  resolution 
refers  to  the  diameter  or  length  of  the  f ield-of-view.  The  field-of- 
view  may  be  defined  as  the  area  exposed  to  the  crystal,  geometrically 
limited  by  the  collimator,  at  any  distance  from  the  collimator  face. 
This  appendix  details  the  calculations  for  f ield-of-view  perimeters 
and  areas,  for  both  cylindrical  and  rectangular  collimators  of 
variable  dimensions. 

B.l  Cylindrical  Collimators 

For  cylindrical  collimators  (Figure  56  )  the  diameter  of 
the  f ield-of-view,  at  a  distance  s  from  the  collimator  face  is  given 
by: 

d  =  2r  +  2s_  (r  +  r  )  (B.l) 

Zh 

where  r^  and  are  the  radii  of  the  collimator  ends  at  the  crystal 
and  at  the  face  respectively,  and  h  is  the  collimator  length.  Note 
that  (B.l)  is  valid  for  all  r^  and  r2  whether  the  collimator  is 

converging,  diverging,  or  straight. 

The  area  seen  by  this  collimator  at  a  distance  s  is: 

A„  -  ITT?  (  3(rl  +  +  l)2  (B.2) 

G  Z  hr„ 


no 


Ill 


did 


RECTANGULAR  COLLIMATORS 
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A  directly  relates  to  the  plane  sensitivity  of  the  detector, 
as  can  be  deduced  from  equation  (2.11). 

B.2  Rectangular  Collimators 

For  rectangular  collimators  (Figure  56)  the  length  of  the 
f ield-of-view  ( £’ )  and  the  width  (w')  are  given  by: 


V 


^  f  ( K  + 


(B .  3) 


w 


T  _ 


W2  +  h  (W2  +  wl> 


(B.4) 


where  and  %  are  the  lengths  of  the  collimator  ends  at  the  crystal 
and  at  the  face  respectively,  and  w^  and  w^  are  the  corresponding  widths 
Although  not  evident  from  Figure  56  ,  the  perimeter  ABCD  includes  the 
rays  from  diagonally  far  corners  of  the  collimator. 

The  area  seen  by  this  collimator  at  a  distance  s  is: 

S^2  SW2(J?+i?) 

=  V2  +  IT  (w2  +  V  +  IT  2  1 

„  (B.5) 


+  -2  U2  +  V  (w2  +  wp 

h 


In  most  cases,  single  detector  crystals  are  cylindrical, 
with  a  circular  cross-section.  Equation  (B.5),  however,  assumes 
the  crystal  to  have  a  rectangular  cross-section.  If  a  circular 
crystal  is  placed  in  a  square  collimator  of  cross-sectional  side 
length  equal  to  2R,  where  R  is  the  radius  of  the  inscribed  crystal, 
then  the  f ield-of-view  area  can  no  longer  be  described  by  equation 
(B.5).  Instead,  a  more  complicated  expression  is  derived  for  the 
locus  of  points  on  a  plane  a  distance  s  away  from  the  face  of  the 


collimator  that  describes  the  perimeter  of  the  f ield-of-view.  From 


1.1  3 


Figure  57  ,  if  the  crystal  face  lies  in  the  z  =  0  plane  in  normal  3- 
space  centered  on  x  =  x  ,  y  =  y  ,  and  if  any  point  on  the  perimeter 
of  the  crystal  face  is  denoted  by  (x^,  y  ,  z  )  then  all  such  points 
P^  satisfy: 


(XD  -  XC)2  +  (yD 


z 


D 


0 


(B .  6) 


The  line  joining  a  point  P  to  any  point  P  with  coordinates  x  y  , 

U  D  O  5 

z  is  described  parametrically  by: 

u 


x  =  ^  +  (xs  -  V" 

y  =  yD  +  (yg  -  yD)t  (b. 7) 


Now,  if  the  collimator  length  is  h,  (Figure  57  )  and  Pg  lies  on  edge  A, 
its  coordinates  satisfy  x^  <  xg <  ,  y^  =  y2>  zg  =  h*  Using  equation 

(B . 7) ,  the  line  P  P  is  therefore  described  by: 

D  JJ 

x  =  [(1  -  a)x1  +  ax  ]t  +  (1  -  Ox^  0  <  a  <  1 
y  =  y^t  +  (1  -  t)yD  (B.8) 


z  =  ht 


In  the  plane  z  =  s  +  h,  t  - 


s  +  h 


,  and  equations  (B.8)  become: 


x 


=  [(1  -  a)x1+ax2](^)  "Iv0-06-1 


/S  +  hN  s  „ 

y  =  y0(— z — )  ”  u  y 


(B .  9) 


h 


h  D 


z 


s  +  h 
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The  point  (x^,  y^,  z^)  obviously  satisfies  equation  (B.6), 
as  it  lies  on  the  circle.  Thus  x^  ,  y^,  and  z^  may  be  eliminated 
using  (B.6)  and  (B.9),  to  yield  an  expression  relating  x,y  to  x  ,  y  . 
This  expression  is: 


{x  -  [xc  +  R(2a  -  l)(S-i-!i)]}2  +  (y 


[yc  +  R^])2  =  ^ 

h 


0  <  a  <  1 


(B. 10) 


z  =  s  +  h 

Equation  (B.10)  therefore  defines  the  perimeter  of  the  f ield-of-view  of 

the  crystax  area  over  one  edge  of  the  collimator,  in  a  plane  a  distance 
s  from  the  collimator  face.  The  other  three  edges  may  be  treated 
similarly,  and  corresponding  f ields-of-view  defined.  In  all  cases, 

g 

the  expressions  define  circles  with  radius  R— ,  but  with  differing 
centers : 


Edge  A;  Centres: 


x  =  x  R(2a  -  0  <  a  <  1 

C  h 


y  =  yc  + 


(B.ll) 


Edge  B;  Centres: 


„  ,s  +  hN 
x  =  x„  -  R( — r  ) 


y  =  y  +  R(2a  -  1)  (24_!l)  ’  0  £  01  -  1 


(B. 12) 
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Edge  C;  Centres: 


x  =  x  +  R(2a  -  1)  (— *  h)  ,  0  <  a  <  1 
c  n  ~  — 


(B . 13) 


Edge  D;  Centres: 


y  =  y„  +  R(2a  -  0  <  a  <  l 

C  h 


(B.14) 


Figure  58  shows  the  f ield-of-view  at  various  distances  s  from  the 
collimator  face  as  calculated  from  these  equations.  As  can  he  seen, 
imposing  the  constraint  that  the  crystal  is  circular  rounds  the 
corners  off  in  the  field-of-view  response  described  by  equation  (B.5), 
and  decreases  the  area  viewed.  The  area  is  given  by: 


y  -  \  -  *)r2-2 


h 


(B .  15) 


where  A  is  defined  in  (B.5)  and  £  -  £ 

R  -L  Z 


w 
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FIGURE  58 


APPENDIX  C 


Mathematics  of  Ray-Tracing 


Assuming  any  two  points. 


PD  =  (V  yD’  V 


(c.l) 


ps  =  (xs’  ys>  zs} 


(C.  2) 


on  the  detector  grid  and  source  grid  respectively,  the  distance  between 
these  points  is 

DDS  =  [(XS  -  V2  +  (yS  '  yD)2  +  (zS  -  ZD)2]  /2 
The  equation  of  the  line  joining  P  and  P  is  given  parametrically  by: 

U  o 


X  =  XD  +  (XS  '  V1 


y  =  yD  +  (yS  -  yD)l: 


(C .  3) 


The  point 


Z  =  ZD  +  (ZS  '  V1 


PA  =  (XA>  yA>  ZA) 


on  the  skin  surface  through  which  the  line  passes  must  be  found, 
Assuming  that  the  skin  surface  is  parallel  to  the  detector  plane  and 
lies  in  the  plane  z  =  z  ,  the  coordinates  of  PA  become: 
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yA  yD  +  ^yS  yD^ZA  ZD^ 


(C.  A) 


The  distance  D  between  P  and  P 

o  A  S 


is  the  length  of  tissue  absorber  through  which  the  gamma  ray  must  pass. 


Neglecting  collimation,  the  remaining  length,  D  -  D  is  the  length 

D  o  OlA. 


of  air  absorber  which  the  gamma-ray  must  also  traverse. 

The  crystal  will  be  surrounded  by  a  collimator,  and  edge 
penetration  effects  must  be  taken  into  account.  In  cross  section, 
the  collimator  will  be  a  height  h  above  the  crystal  face.  Whether  it 
is  a  tapered  or  straight  cylindrical,  or  tapered  or  straight  rectangular 
collimator,  its  geometry  may  be  defined  analytically.  Then,  in tersection 
points  of  the  line  P  P  with  the  collimator  can  be  found  and  the  distance 

J  J  O 

through  collimator  material  determined  and  hence  the  attenuation 
calculated. 

C.l  Cylindrical  Straight  Collimator 

Figure  59  shows  a  section  of  the  collimator  of  length  h 

above  the  detector  grid,  whose  axis  passes  through  the  point  (x^, 

y  z  )  and  whose  inner  and  outer  radii  are  R  and  R  respectively. 

-yl  D  1  o 

The  surface  of  this  collimator  can  thus  be  described  by  the  equations. 


(x  -  +  (y  “  ~  R-l 


(C.5) 
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z 
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for  the  inside  and 


(x  -  x  )2  +  (y  -  y  )2  =  R  2 
i  1  o 


(C.  6) 


for  the  outside.  In  both  cases,  the  z-coordinate  is  restricted  to 
the  length  of  the  collimator,  i.e. 


0  <  z  <  h 


(C.7) 


Rq  -  R^  is  the  collimator  thickness  and  h  is  its  length.  Consider 
the  two  points  P  and  P  given  by  (C.l)  and  (C.2).  The  line  that  joins 

JJ  l) 

these  points  is  characterized  by  the  set  of  equations  (C.3).  By 
solving  (C.3)  and  (C.5)  simultaneously,  a  quadratic  equation  for  t  is 
derived: 


t(ys  -  yD)2  +  (xs  -  V2]t2  +  2[(yD  -  yp  <yg  -  V  +  -  xi>(xs  '  V]t 


(C .  8) 


+ 


By  substituting  for  the  differences  that  appear  in  (C.8),  it  may  be 
simplified.  Let 


a  =  x, 


s  - 


^  yS  yD 


(C.  9) 


and 


A  =  Xp  ”  xi 


(C.10) 


C  =  z. 


D 


z 
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Then  (C.8)  is  solved  for  t,  and  gives: 

tE  =  -(Aa  +  B3)  +  [R.2(a2  +  32)  -  (A3  -  Ba)2]1/2 

—  2  T~ 

a  +3  a  +  3 

Using  this  value  of  t,  equations  (C.3)  become  uniquely  defined, 
t  must  be  chosen  so  that  the  z  intercept  is  greater  than  zero.  Since 
will  always  be  inside  the  cylinder  and  on  z  =  z^,  there  will  always 
be  only  one  positive  z  intercept  (assuming  z^  will  normally  be  zero) . 

The  point  so  defined  (E) ,  will  be  the  point  of  entry  of  the  line  P^P.^ 
into  the  collimator.  The  exit  point  (F)  may  be  determined  in  an 
exactly  similar  fashion,  and  is  characterized  by  t  : 

t  =  -(Aa  +  B3)  +  [RQ2(a2  +  32)  -  (A3  -  Ba)2]  /z  (C.ll) 

a2  +  32  a2  +  32 

From  Figure  60  ,  there  are  three  possible  ways  in  which  PDPS  will 
intersect  the  collimator. 


Case  (1)  If  E (x  ,  y_,  z  )  has  z  >  h  then  there  will  be  no  penetration 

E  E  E  E 

through  the  collimator. 

Case  (2)  If  z£  <  h  and  the  exit  point  F(xE>  yE»  Zg)  which  is  uniquely 
defined  by  substituting  (Coll)  into  (C.3)  has  z^  >  h,  then 
the  ray  P  P  passes  through  the  top  face  of  the  collimator 
and  the  new  coordinates  of  F  will  be  given  from  (C.4)  as: 


-  *D  +  (XS  -  V  (h  ~  V 


(zs  '  V 


yF  -  yD  +  <ys  -  yD} (h  '  zd) 


(zs  -  V 


(C. 12) 
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CASE  I 


Zr>  h 


CASE  2 


CASE  3 


FIGURE  60 


124 


z 


F 


h 


Case  (3)  If  and  z^,  are  both  less  than  h,  the  ray  passes  through 

the  outer  side  of  the  collimator,  and  equations  (C.3)  and 
(C.ll)  are  used  to  determine  F. 

For  either  one  of  the  last  two  cases,  the  absorber  thickness  will  be 

T  =  t(xF  -  X^2  +  (yF  -  yE)2  +  (zF  -  ze)2J  7  2  (C.13) 


C.2  Cylindrical  Tapered  Collimators 
C.2.1  Cylindrical  Diverging  Collimator 

Knowing  the  equation  of  a  conic  in  3-space,  the  inside  surface 
of  the  cylindrical  tapered  diverging  collimator  in  Figure  61  can  be 
described  by: 


(x  -  x^)2  +  (y  -  y^) 2  =  (z  +  D  -  h)2 


(C. 14) 


R 


R 


D 


where 


D  =  R  tan  (j). 


The  line  P  P  will  again  intersect  the  conic  at  one  point 
13  S 

only,  using  the  same  constraints  as  before.  Since  the  outside  of 
the  collimator  is  a  cylinder,  the  outer  surface  is  given  by  (C.6). 
The  finite  length  of  the  collimator  (h)  constrains  z  as  in  (C. 7) . 
Solving  (C.3)  and  (C.14)  for  t,  and  using  the  substitutions  (C.9) 
and  (C.10)  where  z^  =  h  -  D,  a  quadratic  equation  in  t  is  formed. 

[a2  +  g2  _  *  y2]t2  +  2 [Aa  +  B3  -  R2  CY]t  +  [A2  +  B2  -  R2  C  ]  =  0 
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FIGURE  61 
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The  solution  for  t  becomes : 


tE  “ 


(Aa  +  B3  - 


R  Cy)  +  [R  (Ay  -  ca)2  +  R2(By  -  C3)' 

2  2  2 
D  D“  D 


2  1  /  2 
(A3  -  Ba)  ]  72 


a2  +  B2 


R2y2 

D2 


a2  +  32  -  R2y2 
D2 


t  must  again  be  chosen  so  that  z  from  (C.3)  is  positive.  The  value  for  t 

Ej 

for  the  exit  point  F  will  be  given  by  (C.12)  for  Case  (2).  For  Case  (3), 

F  will  be  given  by  (C.ll).  The  absorber  thickness  T  is  given  by  (C.13) 
for  both  cases. 


C.2.2  Cylindrical  Converging  Collimator 

For  a  cylindrical  tapered  converging  collimator,  from 
Figure  61,  the  equation  of  the  conic  describing  the  inside  face  will  be: 

(x  -  xx)2  +  (y  -  y1)2  =  (z  -  D)2 


The  solution  to  the  intersection  points  will  be  the  same  as  for  the 
diverging  collimator  except  that  instead  of  z^  =  h  -  D,  one  needs  to 
substitute  z^  =  D. 

C.3  Straight  Rectangular  Collimator 

For  a  straight  rectangular  collimator,  the  sides  are  described 

by  equations  of  planes,  all  parallel  to  the  z  axis  (Figure  62  ) .  If  all 

the  planes  are  considered  to  be  infinite,  the  analysis  becomes  simpler. 

However,  any  line,  P  P  ,  can  then  intersect  all  the  planes.  In 

1) 

order  to  avoid  this,  the  proper  side  of  the  collimator  involved  in 
penetration  can  first  be  determined  by  certain  selection  rules.  In 
Figure  62,  looking  down  on  the  x-y  plane,  the  collimator  can  be 
divided  into  four  quadrants.  If  the  point  Pg  given  by  (C.2)  is 
projected  onto  the  Z  =  zD  plane,  (which  will  normally  be  the  z  =  0  plane) 


I 
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PLANE  m  n  p  k 

G.  OlOy. 

1  2 

G  0  1  0  y 

o  '  1 

1  0  0  x3 

H  1  0  0  x, 

o  4 

I.  0  1  0  y3 

I  0  1  0  y, 

o  J  4 


FIGURE  62 
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the  z  coordinates  of  the  points  P  and  P  are  equal  and  can  be  ignored. 

U  u 

The  equations  for  x  and  y  in  (C.3)  can  be  solved  simultaneously, 
eliminating  t,  to  give  x  in  terms  of  y,  or  vice  versa: 

(y  -  yD)(xj,  -  V  -  (x  -  V(yD  '  ys> 

If  x^  4  Xg  and  y^  4  y^>  four  possibilities  exist: 


Case  (A)  ^  >  xg ,  yD  >  yg 


Possible 

sides  involved: 

V 

J 

o 

,  G  , 

l 

G 

o 

Side  J.: 

l 

(yD  '  y2) 

(xd 

■  V 

> 

(XD  * 

V 

(yD 

~  ys} 

Side  G. : 
i 

(yD  -  y2) 

S 

-  V 

< 

(XD  ' 

V 

(yD 

“  ys} 

Side  J  : 

o 

(yD  -  yp 

(XD 

-  XS} 

> 

(XD  - 

Xl} 

(yD 

-  V 

Side  G  : 
o 

(yD  -  yp 

(XD 

-  V 

< 

(XD- 

Xl} 

(yD 

-  ys} 

Case  (B) 

*D  >  V  5 

r  < 

D 

ys 

Possible 

sides  involved: 

V 

I 

0 

’  Ji’ 

J 

o 

Side  I . : 

l 

(y3  -  yD) 

(XD 

■  V 

< 

(XD  - 

V 

(ys 

”  V 

Side  J . : 

l 

(y3  -  yD) 

■  V 

> 

(XD  - 

V 

(ys 

'  V 

Side  I  : 
o 

(y4  -  V 

(XD 

-  xs} 

< 

(XD  - 

Xl} 

(ys 

"  V 

Side  J  : 
o 

<y4  -  V 

(XD 

-  xs} 

> 

(XB- 

Xl} 

(ys 

-  V 

Case  (C) 

*d  <  V  yD  > 

ys 

Possible 

sides  involved: 

:  G., 

G  ,  H  , 
o  1 

H 

o 

Side  G. : 

i 

(yD  '  y2) 

(xs 

-  V 

< 

(x3  ‘ 

V 

(yD 

"  ys} 

Side  H . : 

l 

(yD  -  y2) 

(xs 

-  XD) 

> 

(x3  - 

V 

(yD 

*  ys} 

Side  G  : 
o 

(yD  '  yP 

(xs 

-  V 

< 

(x4  ' 

V 

(yD 

-  ys} 

Side  H  : 
o 

(yD  -  yp 

(xs 

-  V 

> 

(x4  ‘ 

V 

(yD 

-  ys} 
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Case  (D)  <  xg,  yD  <  yg 


Possible 

sides 

in vo lved: 

H., 

l 

H 

o 

,  I.,  I 

1  o 

Side 

H  : 

l 

(y3 

-  V 

(xs  ' 

y 

> 

(X3  '  XD) 

(ys 

-  V 

Side 

I.: 

i 

(y3 

-  V 

(xs  - 

V 

< 

(x3  -  y 

(ys 

"  V 

Side 

H  : 
o 

(V 

V 

(xs  - 

y 

> 

(X4  -  XD) 

(ys 

-  V 

Side 

I  : 
o 

V 

V 

(xs  - 

y 

< 

(x4  -  y 

(ys 

-  V 

In  every  case,  only  one  outside  plane  and  one  inside  plane  will 
satisfy  the  conditions.  In  all  four  cases,  if  the  left  hand  side  equals 
the  right  hand  side,  then  the  projection  intersects  the  corner  of 
the  rectangle. 

Other  cases: 


*D 

= 

xs’ 

yD 

> 

ys 

=>  sides 

G. 

l 

and 

G 

o 

XD 

:= 

xs> 

yD 

< 

ys 

=>  sides 

I. 

i 

and 

I 

o 

*0 

> 

V 

yD 

= 

ys 

=>  sides 

J. 

l 

and 

J 

o 

< 

xs’ 

yD 

a 

ys 

=>  sides 

H. 

l 

and 

H 

o 

XD 

= 

V 

yD 

= 

ys 

=?>  no  penetration  factor  possible 

Once  the  penetration  sides  have  been  determined,  since  any 
plane  can  be  described  by 

mx  +  ny  +pz=k  (C.15) 

where  the  values  of  m,  n,  p,  k  are  given  in  Figure  62  for  each  of  the 
planes,  then  the  value  of  t  from  (C.3)  will  be 

k-mxD  -  nyD-pzD _ 

m(xs-xD)  +  n^S_yD^  +  P^ZS_ZD^ 


t 


(C. 16) 
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The  entry  and  exit  points  E  and  F  will  then  be  given  by  (C.3)  for 
Case  (3),  where  an  inside  wall  will  be  used  for  E  and  an  outside  wall 
for  F.  For  Case  (2),  if  <  h  and  z  >  h>  the  point  F  will  have 
coordinates  given  by  (C.12).  The  absorber  thickness  T  is  given  by 
(C . 13) 

C.4  Tapered  Rectangular  Collimators 

The  last  case  to  be  considered  is  that  of  a  tapered  diverging 
rectangular  collimator.  Normally,  only  two  opposing  sides  will  be 
angled,  the  other  two  opposing  sides  being  left  parallel  to  the  z 
axis  (Figure  63^*  The  problem  of  choosing  the  proper  intersection 
planes  still  exists,  but  the  method  of  selection  is  different  since 
the  equations  of  the  tapered  planes  are  more  complicated.  t  can  be 
determined  from  (C.16)  for  both  H  and  JR,  and  the  corresponding  values 
of  (x^j  y^,,  Zg)  calculated  from  (C.3).  If,  for  example,  the  zp  value 
of  plane  1L  is  negative,  assuming  z^  =  0,  the  plane  JR  must  be  used. 

On  the  other  hand,  if  the  z  value  of  plane  is  negative,  the  plane 
1L  must  be  used.  One  of  the  intersection  points  must  be  positive, 
while  the  other  is  constrained  to  be  negative.  After  one  of  the  two 
planes  has  been  chosen,  its  corresponding  y  value  is  checked  to  see 
that  the  condition 


y2  -  yE  1  ^3  (C.17) 

holds.  If  this  is  not  the  case,  then  no  penetration  of  the  tapered 

planes  is  possible.  If  (C.17)  is  true,  a  final  check  is  made  to  see 

whether  zn  <  h.  If  this  constraint  also  holds,  then  the  entry  point 
E 

E  has  been  found.  If  >  h,  no  penetration  of  the  collimator  occurs 
(Case  (1)).  Assuming  an  entry  point  has  been  found  on  H±  or  y, 
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y 


2 

1 

3 

4 


z 
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then  t  is  calculated  for  or  respectively,  and  (x^,  y^,  z^) 
determined.  The  condition 


yx  5  yF  <  y4  (c.is) 

is  checked  for  validity.  If  it  holds,  then  z  for  the  outer  plane  is 

r 

compared  with  h.  For  Case  (2),  F  is  given  by  (C.12),  and  for  Case  (3), 
the  value  of  t  for  the  outer  plane  is  inserted  into  (C.3)  to  yield 
(xp,  yF,  zp). 

When  (C.18)  does  not  hold  for  y  but  an  entry  point  E  exists 

r 

on  1L  or  J\  then  the  ray  must  pass  through  Gq  or  I  respectively.  A 

simple  test  y^  >  y^  or  y^  <  y^,  determines  if  the  plane  is  Gq  or  Iq 

respectively.  Once  this  is  done,  t  for  G  or  I  can  be  found  from  (C.16) 

o  o 

and.  the  exit  point  F  determined  from  (C.3).  Again  if  z  >  h,  F  is  given 

r 

by  (C.12). 

If  (C.17)  is  not  true,  then  the  ray  must  pass  through  either 
G.  or  I..  Again,  y  >  y  or  y  <  yQ  determines  if  the  plane  is  G  or 
I  ,  respectively.  The  entry  point  E  on  the  plane  is  then  found.  If 
z_  >  h,  no  penetration  occurs.  If  z  <  h,  the  entry  point  E  on  the 
inside  plane  is  found.  t  is  calculated  for  Hq  or  Jq,  whichever  applies, 
and  (C.18)  is  tested.  The  same  procedure  as  before  determines  if  one 
of  these  outer  tapered  planes  is  involved.  If  not,  then  G^  or  IQ, 
whichever  applies,  is  used  and  t  calculated  for  one  of  these  and  F  found. 
F  is  determined  from  either  Case  (2)  or  Case  (3).  Finally,  T,  the 
attenuation  distance,  is  given  by  (C.13). 
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APPENDIX  D 

Solid  Angle  Derivations 

The  solid  angle  subtended  by  a  circular  disk  at  an  off-axis 
point  may  be  determined  by  replacing  the  disk  by  an  equal  area,  right 
polygon  with  total  number  of  sides  n,  where  n  is  even  (38  ), 

(Figure  64  )  . 

It  has  been  shown  (38  )  that  the  total  solid  angle  for  the 
n-sided  polygon  for  any  s,  R,  and  p,  is 

(n-2) 

2 

Q(sy  R,  p)  =  2  I  (E  -  F  )  (D. 1) 

r  r 

r=l 


where 


E 

r 


arctan 


2  2  2  2  2  2  2 
s  (x_,  1+y  -2m  x  y  +m  x  +s  +m  x  .  ,+2m  x 


r+1 


r  r  r  r  r 


9  2  \  1  /  2 

«  ~  ,ny  “2m  x  x  ) 

r  r+1  r  r+1  r  r  r+1  r 


m  (s  +x  x  )  — x  y 
r  r+1  r  r+1  r 


(D.2) 


F 

r 


,  2  ,  2  2*/z 
„  s  (x  +  y  +  s  ) 
arctan  r  r 


m  (x  4 
r  r 


)-x  y 
r  r 


(D.  3) 


where 

m^  =  tan  [  (n  -  4r)  ] ,  1  <  r  <  ^  2  ' 

i,  p  sin  [  (2  (r-1)  tt]  ,  1  <  r  <  y 

yr  =  _ n_ 

sin  (7r/n) 

=  2R  [  (-)  tan  (-)  ] 
n  n 


(D.  4) 


(D.  5) 


P 


(D.  6) 
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Xr  =  P  ~  dr’  1  -  r  1  |  tt>.7) 

p  =  (D2  -  s2)  /z  (D.8) 


for 


n 

—  even: 


— 1 / 2  P  + 


.TT 

Psm  - 


TTr.  .  -TTr  ,  TT  TT. 

— )  sin  ( —  +  7  ~  -) 

n  n  4  n 


sm 


(Tr/n) 


1  <  r  < 


(n-4) 

4 

(f>.9) 


d  = 
r 

l  /  n 

/  2  P  >  r  =  - 
4 

(D. 10) 

d  = 
r 

.n+4.  .  n 

_dr  (  4  )  -  r  i  2 

(D. 11) 

where 

l  = 

.n+2. 

(  2  )  r 

(D. 12) 

for  odd: 

d  = 

,TT  TT  TTr.  .  .TT  TT 

p  sm  (y  —  —  +  — )  sm  (-  +  - - 

4  2n  n  4  2n 

‘  T>.  1  <  r  <  (n-2> 

r 

sin  (TT/n) 

4 

(D. 13) 

d  = 
r 

0  r  =  (n+2) 

(lb  14) 

d  = 
r 

-d  ,  (2±«  <  r  <  3- 

v  4  -  -  2 

(D. 15) 

where 

£  = 

(n+2) 

2  -  r 

(D. 16) 

It  was  found  that  a  36  sided  polygon  was  necessary  to  give 
0.1%  accuracy  whereas  for  a  1.0%  accuracy,  it  was  sufficient  to  use  a 
20  sided  polygon.  This  however  was  only  required  for s  / D  <  0.1  and 
0.1  <  R/D  <  2.  For  other  cases  the  number  of  sides  could  be  greatly 


reduced . 


- 
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APPENDIX  E 


COMPUTER  PROGRAMS 


RASTER  SCANNER  DATA  ANALYSIS  PROGRAM 


REAL  RDAT1 ( 10000 ) 

I  NT  EGER  ENDPT(1009),DEFT(200),l(20Q),NDEFT(200) 
REAL  A1 (  100  0  0  ), ARC  100  00  ), AS(  10  000  )/CVAL( 20) 

INTEGER  EOS, EOS 1(20), ROSS, FI NAI  ,  PARTI , PART 2 , I  0^  C  5 ) 
DIMENSION  FORM (10) 

REWIND  2 


DATA  INPUT  FROM  MAGNETIC  TAPE 


WR  I TE ( 6 , 1 2 4 ) 

READ (5,100)  FORM 
DO  1  0=1,100 

READ( 2, FORM , END  =  3 , PPR  =  2  0  0  )  A, RDAT1 ( 0 ) , AB 
IF(A.NE.O.O)  GOTO  2 
GOTO  1 

RDAT1 ( 1 ) =RDAT1 ( .H 
END  PT ( 1 ) = J 

WRITE  (6, 115)  RDATKl  ),ENPPT(1) 

GOTO  3 
CONTINUE 
KK=  0 

WRITE (6, 127) 

READ ( 5,115)  XI 
IF(XL.NE.O.O)  XI  =2 . 1 0  6  2  F  -  6 
K  =  2 
KK1  =  1 
KK2  =  0 
I  KJ  =  0 


CHECKS  FOP  AND  STORES  END-OF-SCAM  POIMTS 


17 


5 


DO  4  « J  =  K ,  10  0  0  0 
READ( 2 , FORM, ENP=5 , FRR=201 ) 

I  FCRDATH  J)  .GT.10.0*PPAT1( 

RDAT 1 ( 0 ) =RDAT1 ( U ) * E X P ( X  I  *10 . 0*( 

GOTO  8 


A , RDAT1 ( J ) , AB 


-1  )  ) 


P  P  A  T 1  ( 
T-]  )  ) 


)=pn  ATI  (,l)/l  0  0.0 


K=J 

K  K  =  K  K  + 1 

I FCKK.GE.7)  GOTO 6 
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201 
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K  =  J 

GOTO  7 

S  IF(AB)  10,10,11 

11  I K  J  = I KJ  + 1 
EOS1  ( I KJ ) =  J 

10  I F ( A )  4,4,0 

9  KK1=KK1+1 

ENDPT ( KK1 ) = J 

I F( ENDPT( KK1 ) - ENDPT ( KK1-1 ) . NE . 1 )  GOTO  4 
KK2  =  KK2  + 1 

DEFT( KK2 ) =ENDPT( KK1-1 ) 

W R I T E ( 6 , 1 1 0 )  DEFTCKK2) 

4  CONTINUE 
7  CONTINUE 
GOTO  17 
C 

c 
c 

C  ALL  RECORDS  HAVE  BFEN  READ  lM,  DECAY  COPPPOTFD,  AND  THE  SUM 
C  DET.  1  &  DET.  2  HAVE  A!  SO  BEEN  RECORDED  IN  RP  ATS .  K  IS  THE 
C  ADDRESS  OF  THE  T n T A L  NO.  OF  RECORDS. 

C 

C 

c 

C  FINDS  I  NIT! A'  SCAN  I  FN0TM5 

C 
C 
C 

f  KD=KK2- 1 

DO  12  J  =  1 , K D 
L( J)=DEFT( J+1)-DEFT( J) 

I  F  ( I  (  D.GE.55)  GOTO  12 
L  ( .1 )  =  L  (  J  -  1 ) 

12  CONTINUE 
C 

C 

c 

c  FINDS  MINIMUM  SCAN  LENGTH 

C 

c 

c 

1  K  =  1 

63  I  S  C  A  M  =  L (  IK) 

I K2  = I K  +  1 

M  2  =  0 

DO  13  J  =  I K2 , KD 
I  F  (  I  SCAN-  L  (  J  )  .  F.O  .  0)  M2=M2  +  1 

13  CONTINUE 

I F ( M 2 . G T . 1 0 )  GOTO  64 
I K= I K  +  1 

I F ( I K . EO . ( KD- 1 ) )  GOTO  65 
GOTO  63 

65  WR 1 TE ( G , 1 2  0 ) 

READ (5,117)  M2 

64  LM I N= I  SCAN 


* 


WRITE (6,104)  LMIN 
L  ENGTH  =  LM I f- 2 
WR I  TEC  6, 105 )  LENGTH 
WR I TE ( 6 , 1 22 ) 

HO  14  J  =  l,  |  Kvl 

WR  I TE ( 6 , 1 1 0 )  EOSl(J) 

14  CONTI  NO E 

WRITE(6,123) 

R  EAD (5,110)  EOS 
C 

c 

c 

C  RDAT1  IS  NORMA l  IZrD  WPT  ITS  I  ARSEST  MEMP,  Fp 
0  I.E.  MAX  I  Ml  'M  =  1.0 

C 

c 

c 

XM1=0 . 0 
DO  32  J  =  1 , K 

XM1=AMAX1(XM1/RDAT1(  ))) 

32  CONTINUE 

WRITE (6, 126)  X  Ml 
DO  33  J  =  1 ,  K 
RDAT 1 ( J ) =RDAT1 ( J ) / XM1 
33  CONTINUE 
C 
C 

c 

C  ENSURES  ALL  ENDPOINTS  ARE  MARKED  &  STORED  IN  NO  EFT  (  I 

C 

C 

c 

1=1 

DO  36  J=1,KD 
I  1  =  2 

I PAR  =  DEFT( J  +  l ) - ( DEFTC  J ) +  LM I N ) +3 

I F ( I  PAR. I  E.O  )  GOTO  60 

GOTO  (37,37,37,37,37),  I PAn 

39  PART1  =  (DEFT(J  +  1)-DEFT(,!))*1.0/I1  +0.5 
PART2=PART1- LM I N+3 

I F ( PART2 . GT . 0  )  GOTO  61 
!  PT  =  J 
GOTO  62 

61  GOTOC  38, 38,38,38, 38),  PART 2 

11=11  +  1 

GOTO  39 

37  NDE  FT (  I  )=DEFT( J) 

GOTO  40 

38  NDEFTC I  )=DFFT( J) 

1=1+1 

41  NDEFTC I )=DEFT( J)+PART1 
DEFTC  J ) =NDFFT ( I  ) 

I F( I  1 . EO. 2  )  GOTO  40 
11=11-1 
1  =  1+1 
GOTO  41 

60  DEFTC J+1)=DEFT(J) 

GOTO  36 

40  1=1+1 

36  CONTINUE 

NDEFTC I )=DEFT(KK2) 


‘ 


' 


o  n  o  o  o  o  o 


1  =  1+1 

NDEFT ( I  )=EOS 


C 


CREATES  MATRICES  A1,AS  &  STORES  SCAN  DATA  IN  THEM 
M / N  ARE  THE  COORDINATES  OF  THE  VA! UE  ACM, N ) 


VmiTE(6,125) 

REA 0(5,115)  5 
1=1 
N  =  1 

NR =S* 10.0  +  0.5 
NT  =  NR 

52  1 A=NDEFT( I )-2 
I R  = I A-LENGTH  +  1 
K  =  1 

DO  50  J=l B, I  A 
I S  U  B  =  K  + 1  ENGTH* ( M- 1 ) 

I R  =  K+  LENGTH* ( NR- ]  ) 

A 1 ( I  SUB ) =RDAT1 ( 0 ) 

ARC  I R ) =RDAT1 ( J ) 

K  =  K  + 1 

50  CONTINUE 

I F ( I A*1 . 0 1 . GT . EOS )  GOTO  55 

1  =  1+1 

I C=NDEFT ( I )-l 
I D= I C- LENGTH+1 
N=N  + 1 
NR=NR- 1 
K  =  1 

DO  51  J=  I  D,  I  C 
J  J  = I C - K  + 1 

I  SUB  =  K  + 1  ENGTH* ( r f  —  1 ) 

I R  =  K+  LENGTH* ( NR- 1 ) 

Al(  I  SUB  )  =RDAT1  (  J-l ) 

ARC  I R ) =RDAT1 ( J  J ) 

K  =K+ 1 

51  CONTINUE 

I F ( I C *  1 . 0 1 . GT  .  EOS )  GOTO  55 
1=1+1 
N  =  N  +  1 
MR=NR- 1 
GOTO  52 
55  NZ=N+1 

DO  400  0=1,1  ENGTH 
DO  401  Jv!=NZ,NT 
I Z ERO=J+ LENGTH* ( 00-1 ) 

Al( IZERO)=0.0 
401  CONTINUE 
400  CONTINUE 

I FCNT.GE  .  2*N)  GOTO  71 

01  =  1 

►  ^St8  72 
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71  J1=0 
J  2  =  1 

72  DO  300  J=l, LENGTH 

N I NB  =  ( NT-N ) *01  +  M  * .  I  2 
DO  301  JJ=1,NINB 
IS=J+LENGTH*( JJ-1 ) 

AS  (  I  S ) = A 1 ( IS) 

301  CONTINUE 

N I NM 1  =  ( NT - N  + 1 ) *  J l  +  (M  +  1)*J2 
NINM2=N*J1  +  ( NT- N ) *U  2 
DO  303  J J=N I NM1 , N I MM2 
I  S=.l  + LENGTH*  (  UJ-1 ) 

AS(  I S  )  =  ( A 1 ( I S )  +  A  R ( IS)) *(11*1.0) 

303  CONTINUE 

N  I  N  =  ( N  +  l )  *01  +  (mt-M  +  1 )*J2 
DO  304  00=010,01 
I  S  =  0  +  LENGTH* ( 00-1 ) 

AS ( I S ) =AR ( IS) 

304  CONTINUE 
300  CONTINUE 

WP !TE(6,111) 

RE AD (5,115)  XSIZF,YSIZE 
C 
C 
C 

C  SETTING  PARAMETERS  FOR  PLOT 

C 
C 
C 

CALI  PLOTS 

CALL  PLOTCO. 0,1. 0,-3) 

M  =  LENGTH 


NC  =  7 

CVAL(1)=0.1 
CVALC  2 ) =0 . 25 
CVA L ( 3 ) =0 . 4 
CVAL ( 4  )  =  0 , 5 
CVAL(5)=0.6 
CVAL ( G ) =0 . 7  5 
CVAL ( 7 ) =0 . 9 
I  OPC 1 )  =  1 
I  0  P  (  2  )  =  1 
I  OP ( 3 ) =  0 
I  0  P  (  4  )  =  1 
I  OPC  5 ) =15 

CALL  MAPYCXS I ZE, YS I ZF , A1 , M , NT, CVA  I 
CALL  PLOTCC. 0,1. 0,-3) 

CALL  MAPY ( XS I Z  E , YS I Z  E , AS , M , NT, CVA  I 
009  CALL  PLOT (24.0,0.0,090) 

WR I TE ( 6 , 1 2 1  )  I  ENGTH , NT, LENGTH , N 
STOP 


NC,  I  OP) 
NC, I  OP) 


100 

200 

101 

103 


FORMAT ('  * , 1 0 A4 ) 

WR I TE ( 6 , 1 0 1 )  J 
FORMAT ( 1  THERE  HAS 
2D'  ) 

FORMAT (  '  ' , 10X, I  4, ' 


BEEN  AN  rRpON  IN 
prcupns  HAVr  BERN 


READING  THE 
p  r  Ap  |M») 


; H  RFCOR 
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104  FORMAT  (  '  ' , 1 0  X ,1 3  /  HAS  R  p  F  N  PETFnMlf'FP  AS  THF  MINIMUM  SCAN  LFNGT 
2!!'  ) 

105  FOpMA  T  (  '  '  ,10X,I3/  HAS  BE  Ep  P  F  T  F 0  M  I  p  E  0  AS  THF  SCAM  LENGTH  IN  THE 

2  X  D  I  R  E  C  T  I  0  N  9  ) 

106  FORMAT ( '  * /OX/THERE  ARF  ',14,'  RrCORPS  lM  THIS  SCAM '/////) 

110  FORMAT ( '  ’,14) 

111  FORMAT  (  '  '/THE  PIOTTIMG  IS  NOW  BEIPC  I”  I  T  |  A  1  I  ZFD  /, 

2'  ENTER  X  BY  Y  RIOT  P I  ME^'S  l  OT?S  <  INCHES/) 

112  FORMAT ( 1 0  * , *  MO .OF  CONTOURS  PFS!pFP?  < I  2  > ' ) 

113  FORMAT ( '  ' ,12) 

114  FORMAT  ('O',  '  VAIUF  OF  CONTOUR  NO.  '  ,12/  ?') 

115  FORMAT ( '  1 ,  3  G  8 . 3  ) 

116  FORMA  T(  1  0  1  ,  'MINI  MUM  NO.  OF  POIMTS  PEp  COr.'TOUR  PESInEP?  <  I  3  >  *  ) 

117  FORMAT/  '  ,  I  3  ) 

62  WRITEC6/18)  I  PT 

118  FORMAT (  1  0  '  /  STAN  ENDPOINT  ANOMAI.Y  AROUND  RECORD'/ 5) 

120  FORMATC  '  0'  /SCAM  LENGTH  CANNOT  BE  DFTFRM I N  FP  -  PLEASE  ENTER  IT:/ 

2  '  <  I  3  >  '  ) 

121  FORMATC  '  0' /THE  PI  OT  SIZE  IS  '/3/  BY  ',13/  -  X  BY  Y  *  / 

2 'THE  SCAN  SIZp  IS  '  ,13/  BY  '/3/  -  X  BY  Y') 

122  FORMATC 'O' /THE  FO!  I  OWING  ARE  THF  POOS!  RIF  EMP-OF-SCAN  PTS./) 

123  FORMATC ' 0' /  ENTER  COnRrCT  E^n-OF-SCAN  FROM'/ 

2' ONE  OF  THE  ABOVE  :  <  I  4  > 1  ) 

124  FORMATC 1 0 '  /  ENTER  PATA  INPUT  FORMAT ' ) 

125  FORMAT  (  1  Q  f  /  ENTER  DETECTOR  SEPARATION  IN  INCHES') 

126  FORMATC  '  0  '/  THE  MAXIMUM  COUNT  IS/F8.0) 

127  FORMATC ' 0 '/  XENON  SOURCE  -  ENTER  1.0  ,  POSITRON  SOURCE  -  0.0') 

STOP 

END 
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CYLINDRICAL  COLLIMATOR  MODELLING 
PROGRAM 


FTN 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


2 


PROGRAM  CLMOD 

DIMENSION  XIB0NC9)  ,XFB0N(9)  ,  Y0ONEC9)  ,ZB0NE(9)  ,RB0NE(9)  ,DR(b0) 
IN  I  EGER  LU (5) ,DATE (9) , TAPE, ISR (100) , IDR  (50) , IT1 (5) , IT2 (5) 

REAL  MUAIR,MIJT  1 S , MU8QN , MUCQL 
D A  T  A  ILU / 1 / , T  APE / 9/ , DR/ i 30*0 „ 0/ 

CALL  RMRaRCLU) 

IF CLU„GT,0)  ILUblU 


addresses 


meaning 


UNITS 


CRR  CRYSTAL  RADIUS 

H  COLLIMATOR  HEIGHT 

R1  cyl,  COLL,  crystal  End  RADIUS 

CYL,  coll.  FACE  end  RADIUS 
T  COLLIMATOR  thickness 

TSCES  TISSUE  SURFACE-CRYSTAL  FACE  SEP, 

S  V  1  ,  S  Y  1  ,  $  2  l 

SX2  ?  3Y2  #  SZ2  SOURCE  PLANE  CORNER 

SX3,SY3fSZ3  COORDINATES 

SX4,  SY4,S7.4 

NB  NUMBER  OF  ATTENUATING  BONES 

X I  BON  Cl)  INITIAL  X  COORDINATE  OF  BONE  I 
XFbON(I)  FINAL  X  COORDINATE  OF  BONE  I 
YROMt(I)  INITIAL  Y  COORDINATE  OF  BONE  I 

ZBONECI)  INITIAL  Z  COORDINATE  OF  BONE  I 

RBONECI)  RADIUS  OF  BONE  I 

MU A  I R  AIR  LINEAR  ATTENUATION  COEFF. 

MUTIS  TISSUE  LINEAR  ATTENUATION  COEFF. 

MUBON  BONE  LINEAR  ATTENUATION  COEFF, 

MUCQL  COLL.  MATERIAL  LINEAR  ATTEN,  COEFF, 
SIGMA  SOURCE  ACTIVITY  PER  UNIT  GRID 

GREN  GAMMA  RAY  ENERGY 


CM 

CM 

CM 

CM 

CM 

CM 


CM 


CM 

CM 

CM 

CM 

CM 


CM-1 
CM-  1 
CM-i 
CM-1 

PHOTONS/TIME/SO, CM 
KEV 


INPUT  OF  INITIAL  MODEL  PARAMETERS 

WRITE(ILU,90) 

READCILUf*)  ISCN 
WRITE  ULU,  91) 

READCILUf 92)  (DATE  Cl) # 1*1*9) 

WRITE (6, 93)  ISCN, CDATE(I) , I»li 9) 
WRITE  (ILU, 100) 

READCILU,*)  CRR 
WRITE  (6, 200)  CRR 
WRITE (6, 2011) 

WRITE (ILU, 103) 

READULU,*)  H,Kl,R2fR0 
WR I Tfc (6 , 203 1 )  h  »  K  1 , R2 
I  s  R  o  -  R  2 

wPITE (6, 2032)  r 
WRITE  ULU*  106) 


o  o  o 
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KEAOUUJ,*)  $X1,$Y1,SZ1 
WRXTEC6, 2US)  SX1  , $Yl,SZt 
WRITE (ILU,  10'/) 

READ (ILU, *)  $X2,SY2,SZ2 
WRIT £(6,207)  SX2,SY2,SZ2 
WRITE  CILU,  108) 

WEmOCILH,*)  SX3,SY3,SZ3 
WRITE  (  S  ,  2  0  8 )  SX3,SY3,SZ3 
WRITE (1LU, 109) 

W  F  A  D  C  1  LtJ ,  *  j  SX4,SY4,SZ4 
WRITE (6, 209)  SX4,SY4,SZ4 
wRi  TECILU, 10b) 

KF.ADUU!,*)  tscfs 
wRlTt(b,2v»5)  TSCES 
WPITECILU, 110) 

READ  (  ILU » *)  N 8 
IF(N6„EDt0)  GO  I U  3 
WRITE  (b, 210)  NB 

wRlTECau,lin 

DO  4  1  =  1 , N  B 

«PI  TECILU, 1 12)  I 

RE  ADC ILU,*)  XldONCI) , XFSON(I) ,  YBONE(I) ,ZBON£(I) , RHONE (I) 
WRITE (6, 212)  I, X  ISON  Cl) , XFBONCI) , Y 80NE  ( I ) , Z80NE C I ) , RBONE C I ) 

4  CONTINUE 

GOTO  5 

3  WRITE C 6, 211) 

5  WRITE  (ILU,  1040) 
wRITEUUU,  1041) 

READ (ILU,*)  hUAJR 
WRITE (ILU,  1 042  J 
READ C ILU r*)  NUT  IS 
wRITEULU,  1043) 

READCILU,*)  M  U  BON 
wRITECILU,  1044) 

READ C ILU,*)  MUCOL 
wRI  ItU>,2«M0) 

WRITE  (6,20  4  1)  HU AIR 
WRITE  (0,2042)  HUTIS 
WRITE  (6, 2043)  MUBON 
WRI  TE  (6 , 20  44)  HU COL 
WRI  TE  (ILU, 114) 

READCILU,*)  SIGMA 

W  R I f  t  ( 6 , 2  1  4 )  SIGMA 
WRITE  (ILU, lib) 

READCILU,*)  GREN 
WRITE CB, 21b)  GREN 
WRIT£(B,2b0) 


90 

FORMAT  Cl x , " 

91 

FORMAT  (IX," 

92 

FORMAT  C9A2) 

100 

FORMAT C 1  X ,  " 

1019 

FORMAT  (IX,  " 

103 

FORMAT  ( 1  X , " 

SCAN  NUMBER 
OATES  «-M) 


and  outside 

-_M ) 


L  RADIUS?  (CM) :  *-") 

COLLIMATOR  HEIGHT  (CM)  I  4-”) 

HEIGHT,  CRYSTAL  EnO  RADIUS  ,  PACE  END", 

RADIUS"/”  OF  COLLIMATOR  IN  THAT  ORDER", 


2"  RADIUS 
3"  (  C  M  )  :  t-  , 

IMu  FnRMAT(lx,»LINt*H  attenuation  coefficient  used 
1041  FORMAT  ( 1  X  ,  "AIR  :  <-") 


CCM-1) S") 


, 
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10  42  FORM  AT  c  IX,  "TISSUE:  ) 

1043  F  OKMaT  C  1  X  ,  "BONE  :  «-") 

1044  FORMAT  ClX,  "COLLIMATOR  MATERIAL!  «•’’) 

105  FORMAT  (1  X  , ’’TISSUE  SURFACE  -  CRYSTAL  FACE  SEPARATION?  (CM)!  4-") 

105  FORMAT  ( 1  x  ,  "X  ,  Y  ,  Z  COORDINATES  OF  SOURCE  PLANE  PT,  Pi?  (CM)!  *-”) 

107  FORMAT  (1X,"P2?S  «• "  ) 

108  FOkMAT  (1  X  ,  "P3?:  «-") 

109  FORMAT (I X , "P4?S  «.”) 

110  FORMAT ( 1 X , "NO,  OF  CYLINDRICAL  BOnES  DESIRED,  PARALLEL  TO", 

2"  X  AXIS  BETWEEN  SOURCE"/"  AND  CRYSTAL?:  «-") 

111  FORMATClx , "ENTER  INITIAL  AND  FINAL  X  VALUES?  INITIAL", 

2”  Y  AND  Z  VALUES?  AND  THE  RADIUS  FOR  EACH  BONE  (CM):  » ) 

112  FORMAT  (1  X,  "bONE  4  ”,ll,«S  4*") 

113  FORMAT (IX, "ENTER  linear  ATTENUATION  COEFFICIENT  OF  BONE", 

2"  FOR  DESIRED"/"  U A M M A  RAY  ENERGY  C CM-  1)1  <-») 

114  FORMAT (1 X , "SOURCE  ACTIVITY  PER  UNIT  GRID  (PHOTONS/UNIT  TIME)", 

2  "  :  ” ) 

115  F  OKMAT  ( 1  X  ,  "GAMMA  RAY  ENERGY  OF  SOURCE  PLANE?  (KEV)!  «•") 

C 

c 

93  FORMAT  (" i ", 10X , "SCAN  «  ” , 1 3 # 20X , 9 A2 , 5 ( / 1 X )  ) 

2 00  FORMAT ("0" , "THE  CRYSTAL  RADIUS  IS  ",F6,3,"  CM") 

2011  FORMAT ("0" , "T  HE  COLLIMATOR  SHAPE  IS  CYLINDRICAL") 

2031  FORMAT ("0", "COLLIMATOR  HEIGHT  !  "  ,  F  6 , 3  ,  "  CM?","  CRYSTAL", 

2"  END  RADIUS ?",F6. 3,"  CM?"/"  FACE  END  R AD  I  US : " , F6 , 3 , 

3"  CM") 

2032  FORMAT ("0" , "THE  FACE  END  THICKNESS  OF  THE  COLLIMATOR  IS  ", 
2F6.2,"  CM") 

2040  format  C'0" , "LINEAR  ATTENUATION  COEFFICIENTS  USED" 


20  4  1 

2042 

2043 

2044 
2u5 

206 

2  07 
208 

209 

210 

21  1 
212 
21  4 

215 

250 

251 


300 
30  1 

500 

501 
C 


2,"  CCM-n  5") 

FOKMAT  ( 1 X , 26X  ,  "AIR!  ", F9,6) 

FORMAT (1 X , 23 X , "TISSUE?  "  , F9.6) 

FORMAT  (.  1  X  ,  2  5  X  ,  "  8  0  N  E  !  "  ,  F  9 , 6  ) 

FORMAT UX, 10X, "COLLIMATOR  MATERIAL:  ",F9,6) 

FORMAT ("0" , "T ISSUE  SURFACE  -  CRYSTAL  FACE  SEPARATION:  ", 
2  P  6 , 3 ,  "  CM") 

FORMAT ( " 0 " , "X , Y  ,  Z  COUROINATFS  OF  SOURCE  PLANE  PT,  Pi:  ", 
2  3  F  8  „  3 , "  (CM) ") 

F  0  R  M  A  T  ( 1 X , 3  8  X , " P  2  •  ",3F8,3,"  (CM)") 

FORMAT (IX, 38X, "P3J  "  , 3  P 8 , 3 , "  (CM)") 

FORMAT (IX, 3ttX,  "P4?  ",3F8.3,"  (CM)") 

FORMAT ("W", "THERE  ARE",  12,"  BONES  BETWEEN  SOURCE  AND", 

2"  CRYSTAL  TO  ACT  AS  ABSORBERS  AND  SCATTERERS") 

FORMAT ("0", "THERE  ARE  NO  BONE  ABSORBERS  USED  IN  THIS", 


2"  SIMULATION") 

F 0  K  M A  T  ( 1  X  ,  "  BONE  4 " ,  1 1 , "  X  I  S " , F  7 • 3 , "  X  F : " , F 7 , 3 , 

2"  Y I ! ” , F  7 • 3 ,  "  Z I : " , F  7 , 3 , "  RADIUS? ",F7, 3,"  (CM)") 

FORMAT (" 0" , "SOURCE  ACTIVITY  PER  UNIT  GRID  AREA  USED!", 

2  F  1  0 , 1  ,  "  PHOTONS /UN  IT  TIME")  ti 

FORMAT ("0" , "SOURCE  GAMMA  RAY  ENERGY  USED : " , F 8 , 2 ,  KEV  ) 


FORM AT  ("1","  ") 

FORMAT  ("0", "EXECUTION  STARTED  AT 
213,"  SECS"/1 X  ,  10X , "STOPPED  AT: 
313,"  SECS") 

FORMAT  (/ J 
FOKMAT (IX , 16) 

FORMAT (15) 

FOKMAT  (1 X,6F8,3) 


" , 1 3 , "  HRS  " , 1 3 , "  MINS 
13,"  HRS  "  , 1 3 , "  MINS  ", 


ii 

, 


' 


n  n  n  n  n  n  n  n  o  n 
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INITIALIZING  DETECT OK  POINT  PARAMETERS 


CALL  EXEC  C  1  1  ,  TTn 
lMC*CRW/0,2  +  w  t  5 

+  i/J  ,  5 

X  INC  *  INC*  I  ,0/2,  Cl 
IT NC a  INC/ 2 
X  I NO* I NO*  1  ,0/2,0 
1 1  NO  =  I  NO /2 

IF (X  INC. L£, II  NCI  iNCalNCtl 
IFCXINU.LE.IINU)  INO«InOM 

ICNXl=lNC/2 
ICNYIsINU/2 
DZ=0,0 
XC=iU0,0 
Y  C  a  1  0  W  .  0 


INITIALIZING  SOURCE  POINT  PARAMETERS 


NTL=  CSQRT C  CSX2-SX \ ) **2+ ( SY2^SY  1)  **2  + 
2  USZ2-3Z1)  **2}  ) /0.5  +  0,5 
NT*I=  (SQR  T  (  CSX3-SX2)  **2+ (SY3-8Y2)  **2  + 
2CSZ3«SZ2)**2) )/l .0  +  0,5 
NKL=NTL+ 1 

mkwsnTw+I 

C 

C 

I )  0  9  K  m  a  1  ,  N  K  w 

TW  =  CRW-1.  )  *  1 , 0  / N T  to  *  J,  ,0 

X1=SX1+(3X4-SX1)*TW 
V 1 =S Y 1 +  CSY4-SY 1 ) *Tw 
ZiaSZl+CSZ4-SZl)*TW 
X2=SX2+  (SX3-SX2)  *Tw 
Y 2  =  S Y 2  +  (SY3-SY2)  *Tw 
Z2=SZ2+($Z3-SZ2')  *Tw 
I  S  »  0 
C 
C 

00  10  KL=1#NKL 
TL  =  (KL-1) *1 ,0/NTL*l .0 
3X  =  X1  +  CX2’*X1)*TL 
Sv  =  Y 1+  (Y2-T 1 ) *TL 
S7.  =  Z1+  CZ2-Z1 )  *TL 
C 

RHO  =  SQRT (  (SX-XC) **2+  (SY-YC) **2) 

CALL  SQL  CCRR,SZ,RHQ,3e>,SULAN) 

WRITE (5,501)  SOLAN 
I  S  »  I  S  +  1 
C 

c 

J  =- 1 
JJ  =  1 
SR =0.0 
JD  =  0 
W  T  »  0 , 0 


, 

, 


7 

6 

41 

A  Vi 

42 

16 

1  1 

13 

14 

16 

12 


MAIN  CONTROLLING  SEGMENT 
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00  50  lx  si, inc 

L/X=(-ICkxI  (IX-1))*0,4  +  1 14  0,0 
1F(JJ„GE.ING)  G  0  T  0  7 
vJ  a  J  ♦  1 

IJ=(XCkY!+J)*2  +  1 
IF  (II»LE»XnC)  GOTO  6 
J  J  =  J  J  i-  i 


JsJ“l 
11=11-2 
on  s  iy  =  i,  a 

UY=(-C1CRYI+J)  +  (IY-l))*0,4  +  1 0  0.0 

J  0  =  J  u  +  1 


AsuX-XC 
b  *  0  Y  -  V  c 
ALP=$X-OX 
BE  f  =vSY-OY 

GAMsSZ-DZ 

wsGAM/(ALP**2+BET**2+GAM**2)**1.5 
W  T  «  W  ♦  w  T 

U«W*StJLAN*SIGM  A 

IF  ( N  B  ,  E Q  •  0 )  GOTO  Av) 

BONAT>0.0 
XBUNT=0,0 
00  4  1  J  B  =  I | N  B 
bRnDY-YBONE  (JR) 

BC=DZ-ZBun£( JBJ 
KPsKBUNECU 

CALL  BTH(Ra,6B,BC,BET,GAM,8GNTH) 

BONA  r  =  EXP  C-MUBUNurHUMTH)  +  BONA T 

XBONTsXHON  I  +  6  0  N  T  H 

CONTINUE 

GOTO  42 

b  0  N  A  T  =  1 . 0 

XBUNTsa, u 

CALL  TTH  (TSCFS,,  DX,0Y,0Z,3X,SY,SZ,  ALP,BET,GAM,TISTH,XBONT) 

WRITE ( 6 ,  2042)  TISTM 
! ISA  I  =  tXP (-MUT 1S*TISTM) 

GOTO  11 

cola  t=exp  c-mucol*colth) 

WRITE (6, 2W44)  COL  TM 

CALL  ATHtALP,BET,GAM, AIR rH,TlSTH, BOOTH, CULTH) 
AIRAT=EXP(-MUAIR*AIRTH) 

WRITE  16, 2041)  AIHTh 
GOTO  12 

IF  CK2-R1)  13,14,15 

CALL  CCTHCOX, DY,DZ,SX,SY,SZ, A, B,C, ALP, BET, GAM, T,R1,R2,H,C0LTH) 
GO  It)  16 

CALL  CvSTH  (  OX  ,  i)Y  ,  OZ  ,  -SX  ,  $  Y  ,  SZ  ,  A  ,  B  ,  C  ,  ALP  ,  BET  ,  GAM  ,  T  ,  R1  ,  R2  ,  h,  COlTH) 

GOTO  16 

call  corn  cox , o  y , pz , sx , s y , sz , a , b , c , alp , be i , gah , i , k i , R2 , m , cqlthi 
GOTO  16 

A  F  I  N  =  Q  *  8  0  N  A  T  *  T  I  S  A  T  *  C  0  L  A  T  *  A  1 R  A  T 
DR(Ji))=AFlN  +  DR  (JO) 


. 


n  o  o 


S R  =  A F  I  N  +  SP 


147 


C 

*  CONTINUE 
bw  CONTINUE 

SUN  OVER  INDIVIDUAL  MATRIX  INDICES  AND  DUMP 
ON  PAPER  I  APE 
IFCJD,GTtbtd)  STOP 
DO  17  J  X  =  1  ,  J  0 
I  OR  C JX) =  QR  t JX) /wT 
17  CONTINUE 

ISk CIS! =SK/WT 
C 

I  F  ( I  $  R  C  I S  J  „  L  t  «  4 )  GOTO  HI 
WRITE  (6, 301)  I  SR (IS) 

C 

10  CONTINUE 

WRITE  (4, 50tf)  (ISRUin  , III=l,NKL) 

0  CONTINUE 

C 

WR 1 TE ( 4 , 500 )  (I DK  U 1 1) # II I  *  1 , JD) 

CALL  EXEC (3 , 1004B) 

CALL  EXEC  ( 1 1 , IT2) 

C 

WRITE  (6 ,25l)  TT1  (4)  ,  IT1  (3)  ,  IT  1  (2)  ,  IT2(4)  ,  IT2C3)  ,  IT2C2) 

WRITE (6,300) 

STUP 

EMI) 

C 

C 

C  SOLID  ANGLE  SUBROUTINE 

c 

c 

SUBROUTINE  SOL  (« , H, RMO, NC , SQL  AN) 

DIMENSION  EC90)  f  F  (951)  ,  G(90)  ,0(90)  ,  EM  (94)  , 

2 Y (90) , X (90) 

N»NC 

PI»3. 141S926 
N2*  (N-2) 72 
X  N  *■  N 
N  4  a  M  /  4 
IN  6  a  (  N  +  6 )  /  4 
N444s (N  +  4) / 4 
B  0  a  $  I N  l  P  I  /  X  N ) 

N?1 SN2+1 
N24» (N-2) 74 
N  4  4  =  (  N  -  4  )  /  4 
N 2  41 3  N 24+1 

SPRa  (PI /XN)*  (BO/COS  (PI /XN) ) 

IF  (SRR.LT. 0,0)  ST0P1 
X P  =  2  ,  *  R  *  S  Q  R  T  (SRR) 

C  DETERMINE  IF  n/2  IS  ODD  OR  EVEN 

XN2*N/4 

XNNaxN-4.0*XN2 
IF(XNN)20# 10,20 
C  000, EVEN, ODD 

C  CALCULATION  OF  OCIRJ  IE  N/2  EVEN 

10  1)0  13  1R*1,N44 

X  I  R  =  I  R 

13  DtIR)sXP/2)(0  +  XP*S  IN  (PI/4, 0-PI*XlR/XNj* 
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1  A 


c 


2  0 


2  1 


22 

30 


31 


35 

39 


5 14 

4  1 


51 

45 

42 

43 


C 

c 

c 

c 

c 


1 


2SIN  CRI*X  IR/XN+PI/4,0-PI/XN)  /BO 
0  f  N  4  J  =  X  P  /  2 , 0 

00  14  XR=N444,N21 
1RR  =  (N  +  2) /2-IR 
0  ( I R )  «  -  0  ( 1 R  R ) 
b  o  r  u  3o 

CALCULATION  of  O(IR)  IF  n/2  ODD 

DO  21  IR»l,N24 

XIK3IR 

U(1R)  »  X  P  a  S  I N  (PI/4.0-PI/  (2„0*XN)  +  PJ  *XIR/XN)  * 
2SlN(Pl/4.0-Pl/C2,f4*XNH'PI*XlR/XN)/bO 
0 CN241 ) =w,0 
DO  22  IH  =  Nt>,N21 
1RR  =  CN  +  2 ) /2-IR 
0  ( IR) s-O CIRR) 

CONT  In«JF 
DO  31  lRal,N2 
x  I R  s  i  R 

t  D  c  i  r  )  =  s  i  n  c  c  c  x  n«4 ,0txiR)*pn/(2#o*xNn/ 

2C0SC C  CXN-4„0*X1R) *RI) / C2.0*XN)) 

00  35  XR=i,N2 
X  1 R  s  I N 

X (I  R )  =  (XP/2,0) a  S  I  N  C  (2»0axlK-l ,0) a  P  X  /  X  N ) /BO 
00  39  IR=1,M21 
X  UR)  aRHO-O  (IK) 

DO  41  lRal,N2 

IOP*M*SwRTCxn«)**2  +  Y(IR)**2tHa*2) 

BOTaEM  CIR)  *  (X  UR)  **2  +  H**2)  -X  (IR)  *Y  (IR) 

FUR)  saTan  CTOP/BOT) 

IF  CF CIR) ) 50, 41, 41 
1  ( IR) *F  f IK) +PI 
CONTINUE 
1)0  4  5  I  R  a  1  ,  N  2 

SPR  =  X  CIK+1) **2  +  Y  t  IK) **2-2.0*EM(IR) a 
2x(lK)*YCIR)+CEM(IR)**2)*(XCIR)**2)+HaH+ 

3  t  F  M ( I R ) a  a£  j  a  t  X  ( I R  +  U  a  *2 ) +2 , afc"M ( I R ) a x ( I R+ 1 ) a 

4rUR)-2.0aCEMUR)**2)aX(IP  +  l)*XCIR) 

IF  CSRR.LT.  0.0)  ST0P2 

rnpaHaSQRT  CvSRR) 

60rsEM(IK)a(HaH-t,X(lR  +  l)*X(IR))-XllK,*,l)*Y(IR) 
LUR)  sATAN  CTOP/BOT) 

IFCECIR) )5l,45,45 
£  C  IR)  =H  (HO  +  PI 
CONT I NU& 

DO  42  I R  a  1 , N2 
GUR)  *  £  ( I R  )  -  F  UR) 

5  0  L  A  N  a  0  ,  'A 
DO  43  1  R  =  1  , N  2 
S  0  L  A  N  a  S  0  L  A  N  +•  2  »  0  a  G  (  1 R  J 
RF TURN 
£  N  0 


BONE  THICKNESS  subroutine 


SUBROUTINE  BTHtRBrBB,BC,BET,GAM,BONTH) 

TEST  a  CRB**2)  *  («ETaa2  +  fiAMaa2)-CbET  a  b  C  -  G  A  H  *  b  8 )  **2 

IFCTEST)  1,1,9. 
h  0  N  T  H  a  0  „  D 


i 


CJ  u  u  u  u  «“*  u  u  u  u  u 


R  F  1  U  R  N 

2  BONTrl  =  2.0*SQRT  (TEST ) /  (8ET**2  +  GAM**2) 
-ip  TURN 
END 
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r  ISSUE  THICKNESS  SUBROUTINE 


SUURUUT INE  TTH(T$CFS,OX,OY,OZ,SX,SY,SZ, ALP , BET , G AM , T I  STM , X60NT ) 
1FCTSCPS.UT.SZ)  GOTO  1 
TISTHafl.d 
PF  1  URN 

\  AXBDX+ALP*(TSCFS"DZ)/GAM 
AY*DY  +  6ET*  CT5CFS-OZ) /GAM 
3RRs((SX-AX)**2'f(SY-AY)**2+(SZ-TSCFS)**2) 

TISTH«SQRT (SRR) -X0ONT 
FORMAT  C1X,E10#6) 
return 

fc  N  0 


CYLINDRICAL  CONVERGING  collimator  subroutine 


SUbROLJ TINE  CCTH  (OX  , OY , DZ , SX , S Y , SZ , A , », C , ALP , BET , GAM, T , R 1 ,  R2 , H , 
2 COL TH) 

K  0  s  R  2  + 1 

'  Ks CCR1-R2) /H) **2 
D »  N  *  R 1 / t  R  1 «  R  2 ) 

C  e  "  0 

ENUMls-A*ALP’“8*BE f+R*C*GAM 

$PR  =  R*(.A*GAM-C*ALP)**2  +  R*(H*GAM~C*BET)**2~CA*BET-8*ALP)**2 
IF  CSKR.LT. 0.0)  ST0P4 
E  N U M  2  3  S  Q  R T (SRRJ 
lOEn=alP**2+8ET**2-R*(GAM**2) 

C 

1  FPL3  CENUMltENUM2) /EDEN 
TFMI s  CENUM J-LNUM2) /EDEN 
wPITE  (U,  100)  QX,UY,OZ,SX,$Y,SZ 
wP I  I  £ ( b , 1  MU )  A , B r C t ALP , BET , U AM 

W  P  I  T  E  ( 8  , 10  0)  R , 0 , ENUM1 , ENUM2  , EDEN , SRR •  f  EPL , TfcMI 
lvu  EORMATUX»10F8«3) 
t7=0Z+GAh*TEPL 
IF(EZ)  1,1,2 

1  EZ=DZ+GAM*TEMI 
7  IF(EZ-H)  3,6,6 

C 

c 

2  EZ2  =  0Z  +  GA.M*T£MI 
IFCEZ2)  7,7,3 

3  IFCEZ2-E7J4,7,7 

4  El-El  2 
GOTO  7 

6  COLTH»0,0 
RETURN 
C 
C 

5  rERARs CEZ-UZ) /GAM 
t.  XaOX  +  ALP*  TFP AK 


. 


iYst)Y«-bE  f  *  IFPAR 
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F  N  U  M  1  a  -  CA*ALP+D*bET  ) 

SPK=KQ**2* (ALP**2+bE  r  * *2 )  -  ( A *Bfc T-B*  ALP )  *  *2 
U  lSRRf LT ,  O,0l  STOPS 
f >*UM2  =  SQRT  (SRR) 

F0ENaALP**2+BE  T**2 

C 

J.F  (FOEN,  EU.0, 60  GOTO  13 
1  FPU=  (FNUM1+FNUM2)  /FQEN 
1  FMI*  (FNUM1-FNUM2)  / F {) E N 

fz=dz+gam*tfpl 

c 

l  F  (.  F  Z )  14,10,12 

10  FZ*DZ*GAM*TFMI 

12  IF(FZwH)  11,13,13 
C 

11  TFPAKs CFZ^OZ) /GAM 
C 

fx =ox+ alp*tfp aw 

F  y  a  D  Y  +  a  E  T  *  T  F  P  A  R 
GOTO  14 
C 

1 3  FX=OX+ALP* (H-OZ) /GaM 
F  YsQY+BE  I* (H-OZ) /GAM 
KZ  =  H 

C 

14  COLT rls SORT  ( (FX-EX)  **2  +  (FY-EY)  **2  +  CFZ-EZ)  **2) 

RETURN 

t  N  l) 

C 

c 

C  CYLINDRICAL  STRAIGHT  COLLIMATOR  SUBROUTINE 

C 

c 

SU8RUU TINE  CSTHCDX,DY,DZ,SX,SY,$Z,A,6,C,ALP,BET,GAM,T,R1,R2,H, 
2C0LTM) 

R  cr  1 

C  »  W  ,  0 

X  N  U  H 1 = - (A*aLP+B*BE1 ) 

D  E  N  a  A  L  P  *  *  2  ♦  8  E  T  *  *  2 
IFlD£N.Ed,0f0)  GOTO  6 

1  SRRsR*n2*(ALP**2fBtTa*2)'»(A*BET-B*ALP)**2 

IF (  S  R  R  »  L  I  •0*0)  STOPS 
XNUM2»SGRT  (SRR) 

C 

c 

TPL=(XNUMl+XNUM2)/UEN 
T  M i = (XNUM1-XNUM2) / 0 E N 

w P I T 1 1 6  # 1401  R , XNUMl , XNUM2 , DEN, SRR , TPl , TMl 
100  FORMAT (IX, 10F8.3) 

C 

IF(R-Rl)  2,2,3 
2  EZsDZ+GAM^TPL 
IF  {ED  4,4,5 

4  EZaOZ>GArt*TMI 

5  IFC£Z«M)  7, 6,6 

6  COLTH«0#0 
RETURN 


■ 


O  O  C3 
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7  T  F P  A  R  =  t  E  Z  -  0  Z  )  /  G  A  K 
EX  sDx«- AL.P*  TEPAK 
EYsDY  +  BEr*  TEPAK 

R  s  R  1  +  T 

G  0  10  1 
C 
C 

5  1 Z=0Z+Gam* FPL 

IF  C  F  Z  )  8,8,9 
ri  F7«0Z+GAB*fMl 
9  IFIFZ-H)  1-4,11,11 
C 

10  TFPARr (FZ-UZ) /GAO 
F X=DX*ALP*TFPAR 

P  Y  =  D  Y  +  6  E  T  *  T  F  P  A  K 
GOTO  12 
C 

11  FX=DX+ALP* ( H  -  D  Z ) /GAM 
F  Y  =  D  Y  +  BET  *  C H-DZ )  /  G  A M 
F7=H 

C 

12  CQLTHs8UKTCCFX-EX)**2+CFY'-EY)**2+CFZ-eZ)**2) 

RETURN 

CND 


CYLINDRICAL  diverging  collimator  SUBROUTINE 

c 
c 

SUBROUTINE  CDTM (DX,OY, DZ, SX, SY, SZ, A,B,C, ALP,eET , GAM, T,R! ,R2, H, 
2C01TH) 

K0=R2+T 

k=  t  ( K2-R 1 1 /H) **2 
D  e  m  *  R  2  /  ( R  2  -  R  1 J 
C  =  D-H 

C 

EMUMi=-(A*ALP  +  B*b£T«'R*C/GAM) 
$rR3K*(A*GAM-C*ALP)**2+R*C6*GA0-C*BET)**2 
2  -  (  A  *  6  E  1  -*  B  *  A  L  P  ) *  *  2 
IF tSRR,LT,O.0)  ST0P6 
ENUM2  =  SilRT  (SRR) 

EnEM  =  ALP**2  +  BET**2",R*GAM**2 

c 

T  F  P  L  =  (ENUM1+ENUM2) / E  0  E  N 
TEMIa  (ENtJMl-ENUM2)  /EDEN 

C 

EZsDZ+gAM*VEPC 
IFCEZ)  1,1,2 

1  E  2  s  U  Z  +  G  A  M  *  T  E  M  I 

IFCEZ)  0,6,2 

6  COLTh=0#0 
RETURN 

c 

2  IF(EZ-H)  4,6,6 
C 

4  I  FPARs CEZ-UZ) /GAB 
Ex«DX+ALP* TEPAK 
EY=0Y+6EI*rEPAK 


C 

C 
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F  HUMlc-  (  A*ALP  +  tf*BET) 

SRRsR0*rf2*  (ALP**2  +  BET**2)-CA*QET-H*ALP)**2 

1 F  (  S  R  R  ,  L  T  ,  H  ,  0 )  ST  0 P  7 
FNUM2»SQRT tSKK  J 
FHENb AUP**2  +  BE T *  * 2 
C 

T  F P L  =  C  F N U  M 1+FNUM21/FDE N 
T  F  M I  * (FNUM1 -FNJM21 /FQEN 
FZ=DZVgam* T FPL 
IFCFZ)  10,10,12 
10  FZ«DZ+GAM*TFMI 
C 

12  IF(FZ-H)  1.1,13,13 

13  TFPAKs  (FZ-OZ) /Gam 
C 

FxaDX+ALP*TFPAK 
F  Y  s  f)  Y  ♦  BE  f  *  T  F  P  A  R 
G 0  T U  14 
C 

13  FX  =  OX  +  ALP*  IH-DZ) /gam 
FY=DY*6E r * (H-OZ) /GAM 
F  Z  =  m 

C 

14  COLlHsSQKT  (  (F'X-EX)  **2  +  (FY*-£Y)  **£♦  (FZ-EZ)  **2) 

RE  TURN 

End 

c 

C 

C  AIR  THICKNESS  SUBROUTINE 

C 

c 

SUBROUT  I N £  ATH(ALP,8ET,GAM, A IRTH , T ISTH , BONTH , COLTH) 

0 1 S 3 3 (4 R T  (ALP**2+BET**2+GAM**2) 

AIK  riiauis-  r  ISTH-BONTH-COLTH 

RF 1  URN 

END 

EMU* 
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RECTANGULAR  COLLIMATOR  MODELLING 
PROGRAM 


PROGRAM  K C MOO 


DIMENSION  XIB0N(9)  ,  XK80NC9)  ,YB0NE(9)  ,ZB0NE(9)  ,RB0NE(9)  ,ORC50) 
IN  1  EGER  LU (5) ,UATE (9) ,TAPE, ISR  (100) , I0R  (50) , IT1 (5) , IT2(5) 

REAL  L 1 , L2 , KG  I , NG I ,  KG  I , MHI , NM I , KM  I , M J I , N J I , K J I , M 1 1 , N I  I , K 1 1 
REAL  MGO,NGO,KGQ# MHO,NHQ,KHQ,MIO,NIO,KlO#MJO,NJO,KJO 

REAL  MUAIR,MUTIS#  MU80N, MUCOl 

COMMON  DX,DY,OZ,SX,$Y,SZ, A,B,C, ALP,BET,GAM,T,RI,R2, H,RB,T$CFS 
COMMON  CULTH, AlRTH,TISTH#BONTHf MGI , NG1 r PGI#KGI # MGO , NGO , PGO , K GO 
COMMON  MH I, NHI, Ph I, KH. !, MHO, NHO,PHO,KHQ, MI  I, N IT, P II, K II, MI 0,NIQ 
COMMON  PIO,KIO,MJI ,  NJI, PJI , K J I r M JO , N JO , P JO , K J 0 # XBONT , BB , BC 
DATA  1LU/1/,T ARE/9/, PR/ 130*0,0/ 

CALL  NMPAW(LU) 

IF  (l'J,GT,0)  I  L  U  ®  L  U 


C 

C 

c 

L 


ADDRESSES 


MEANING 


UNITS 


c 

CNR 

CRYSTAL  RADIUS 

CM 

c 

n 

COLLIMATOR  HEIGHT 

CM 

c 

11 

RFC, 

CULL.  CRYSTAL  END  LENGTH 

CM 

c 

w  1 

RFC. 

COLL.  CRYSTAL  END  WIDTH 

CM 

c 

L2 

REC, 

CULL*  FACE  END  LENGTH 

CM 

c 

W2 

REC * 

coll,  face  end  width 

CM 

c 

T 

COLLIMATOR  thickness 

CM 

c 

TSCFS 

TISSUE  SURFACE-CRYSTAL  FACE  SEP, 

CM 

c 

SXI , SY 1 , SZl 

c 

SX2,SY2,SZ2 

SOURCE  PLANE  CORNER 

c 

S*3,SY3,vSZ3 

COORDINATES 

CM 

c 

SXR , SY4 , SZ4 

C  M  G  I ,  N  G I ,  P  G I »  K  G I 

L  MHI  ,  N  H  1  ,  P  h  I  »  K  M I 
C  M I  T  ,  N 1 1 ,  P  1 1 ,  *  1 1 

l  MJI,NJI,PJIrKJl 
C  MGO, NGU# PGO, KGO 
C.  MHO,  (N  H  0  ,  P  H  0  r  A  H  U 

C  M 1 0  ,  ,\|  1 0 ,  P  1 0  ,  K  T  U 

C  MJU, NJU, PJO, K JU 
C  Nb 

C  XT  BON  Cl) 

C  XFbON(I) 

C  Y RUNE  Cl) 

C  ZBONE(I) 

C  RBUNECI) 

L  MU A  IK 

C  MUTiS 

C  MUBUN 

C  MUCOL 

C  SIGMA 

C  GKEN 

C 


RECTANGULAR  COLLIMATOR 
PLANE  PARAMETERS 


NUMBER  OF  ATTENUATING  BOnES 
INITIAL  X  COORDINATE  OF  BONE  I 
FINAL  X  COORDINATE  OF  BONE  I 
INITIAL  Y  COORDINATE  OF  BOnE  I 
INITIAL  l  COORDINATE  OF  BONE  I 
RADIUS  OF  BONE  I 
AIR  LINEAR  ATTENUATION  COEEE, 
TISSUE  LINEAR  ATTENUATION  CUEFFr 
bone  LINEAR  ATTENUATION  COEEF, 

cull,  material  linear  atten,  cqeff 

SOURCE  ACTIVITY  PER  UNIT  GRID 
gamma  ray  ENERGY 


CM 


CM 

CM 

CM 

CM 

CM 

CM»  1 
CM-1 
CM*  1 
CM-1 

PHQTQNS/TIME/SQ.CM 

KEV 


t 


n  o  o 
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INPUT  or  INITIAL  MODEL  PARAMETERS 

WRITE (l'LU# 90) 

R  E  A  D  ( 1 L U #  * )  ISCN 
WE1TE(ILU,91 j 

ft  E  A  l.)  ( I L  U  #  9  2 )  ( U  A  7  E  ( I)  #  I  s  1  #  9 ) 

WRITE (6, 93)  ISCN,  (DATE  CD  ,  I  a  1 , 9) 

WRITE  (XL U,  im) 

READCILU#*)  CRR 
WRITE (6, 200)  CRR 
WRITE  C6, 2(412) 

WR1 T ECIL. U  t  1019) 

K  F  A  Q  (  I  L  U  #  *  )  H 
WRITE C  ILU ,  1020) 

READ  ClLU ,*)  MG  I, NG1 , PGI  , KG  I 
WRITE (ILU# 1021) 

RE  A  D  C I L  U  »  * ) 

WRITE  (ILU#  1022) 

READ  (ILU,*)  M  1 1 , M I , P 1 1 , k ii 
PPITECII Ur 1023) 

reaucilu#  *)  mjdnji,pji,kjj. 

WRITE C ILU, 1024) 

RE  AD  C  ILU #  * )  MGO , RGP , PGU, KGU 
WRITE  (ILU# 1021) 

READ (ILU# *)  MHO, NHU# PHO, KHU 
WRITECILU# 1022) 

H  E  A  D  U  L  U  ,  *  )  M  I  U  ,  N  I U  ,  P  1 0  #  K  I  U 
TRITE  (ILU, 1023) 

READCILU# *)  MJU,nJO,PJO,k JO 
WRITE  (6, 2021) 

WRITE (6, 2022)  MG  I  , NGI , PGI , *GI , MHI , NHI , PHI ,  KHl , 

2M I  I  , N I  I  , P  1 1  #  K  II  , M J I , N J I , P J I , K  J I 
WRITE (6,2023) 

WP I T E ( 6 , 2022 )  MGU  r  N60  #  PGO , KGO , MHO  #  NHO , PHO  #  KHO , M 1 0 , N 1 0  #  P 10  r  K 1 0 
2  #  MJO#  NJO, PJO,K JO 
L 1 =AHS (KHI-K J T ) 

L2=L1  +  AHS (H*  (PJI-PHI) ) 

W  1  =  ADS (Kll-KGI) 

w  2  =  w  1 

r=A6S (rgi-rgo) 

write (6, 2024)  H,L1 # Wl #L2# W2 
wPlTt (6, 2032)  T 
WRITECILU# 106) 

READCILU#*)  $  X  1  #  S  Y 1 #  SZ 1 
WR ITE (6, 206 )  SX1 , S  Y  1  #  SZ 1 
WRITE (ILU#  107) 

KFADULUi  *)  SX2#SY2#SZ2 
WRITE (6 , 20/)  SX2#  SY2#  SZ2 
WRITE ( 1 L  U  # 10b) 

READCILU#*)  SX3#SY3#SZ3 
wRITE(6#20H)  S  X  3  #  S Y  3  #  S 1 3 
WRITE ( ILU# 109) 

READCILU, *)  SX4,SY4,$Z4 

WP ITE (6 , 209)  SX4  #  SY  4 , SZ4 
WRITE (ILU# 105) 

READCILU#*)  TSLES 
WRITE (6,205)  TSCES 
WRITE  C ILU, 110) 

READCILU#*)  NB 
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1F(NBWE0„W)  GOTO  3 
w  R 1  J  t C  6 , 2 1 4 )  Mb 
w R 1 T t  Cl L U  ,111) 
on  4  1*1 #NH 
w  P  1  T  t  (  I  H 1 ,  1  1  2  )  x 

'  f  ?  ?  t  !LU:  *?  x  I  OOO  ( 1  )  ,  XFBON  ( I  )  ,  YHONE  (  [)  ,  ZB  ONE  ( I )  , RBONE  (1) 

A  y-]?)  1  '  *  1  HU',J  C  1  )  •  *F»ON  CD,  YbONE  II  )  ,  ZDONfc.  (I),  RBQNE  (  I ) 

on  to  o 

3  wPITfcCb,2il) 
o  *«  R  I  T  li  c  1 1 .  U  ,  1 1)  4  0  ) 

W  PITeULU#1041) 

RFAOULU,*)  MUAIK 

s  R  1  T  t  l  I  L  U  ,  1 0  4  2  ) 

RfADULU,*)  MU  TIS 
WRITE.  (ILU,  1043) 

R  F.  A  D  C 1 L  U  ,  * )  M  U  b  0  N 

>'■  P  1  T  t  (  1 1 U  ,  10  4  4,) 

RFADULU,*)  MUCOL 
wRiTfcl(?,20  40) 
w  R  1  T  E  ( 6 , 2  0  4  j  )  M  ll  AIR 
WRITE  (t>f2M42)  mu T  K, 
vi R  1  T  t  ( b ,  2 0 4 3 )  HUB 0 h 
wR  ITE Co , 2044)  MUCOL 
WRITE C ILU,  114) 

rfapulu,*)  SIGMa 
W  R  1 T  E  ( 6 , 2  1  4 )  SIGMa 
*RiTE  CJLUr 1 15) 

KFADCILU,*)  GREN 
■NplTE(e,21ij)  GWEN 
WPI Td  (6, 2S0) 

C 

C 

C 


90  FORMA)  ax  ,  "ENTER  SCAN  NUMBER  !  «-») 

91  f- o k n at  tix,  "enter  dates  «*") 

92  FORMAT (9A2) 

100  FORMAT  (IX,  "CRYSTAL  RADIUS?  (CM)  :  *") 

1019  FORMAT  (IX,  "ENTER  COLLIMATOR  HEIGHT  (CM):  <-") 

1020  FORMAT  U  X , ’'ENTER  COORDINATES  OF  RECTANGULAR  PLANES  -  ",C/1X), 
2"  INSIDE  PLANES:  M  n  P  K",C/lX),7X, 

3  "PLANE  G  S  «-") 

10  21  FORMAT  ( 1  X  ,  7X  ,  "PLANE  HI  4-") 

1022  FORMAT  (IX, /X,  "PLANE  IS  «-") 

10  23  FORMAT  U  X,7X,  "PLANE  JS  *■") 

1024  FORMAT  UX,"  OUTSIDE  PLANES!  M  N  P  K",(/1X),7X, 

2  "Plane  G!  «-") 

1040  FORMAT ( 1 X , "LINEAR  ATTENUATION  COEFFICIENTS  USED  CCM-1):") 


1041  FORMAT  (IX,  "AIR!  «-") 

1042  FORMAT  (IX,  "TISSUES  ♦•") 


10  43 
10  44 
105 
10b 
107 
1 0  8 
109 
1  l  0 

1  1  1 


FORMAT  ( 1  X  ,  "8UNE  :  <-") 

FORMAT  (IX,  "COLLIMATOR  MATERIALS  «•") 

FORMA T ( 1 X , "TISSUE  SURFACE  -  CRYSTAL  FACE  SEPARATION?  (CM): 
FORMAT (lx,  "X, Y,z  COORDINATES  OF  SOURCE  PLANE  PT.  PI?  (CM): 
format  ( l  x  ,  "P2? :  «-") 


F  OKMAT  (IX,  "P3?  S  ♦’") 


F  0  R  h  A  T  (  1  X  ,  "  P  4  ?  S  «* "  ) 

FORMAT (lx, "NO,  OF  CYLINDRICAL  bONES  DESIRED 
2"  X  AXIS  BETWEEN  SOURCE"/"  AND  CRYSTAL?: 
FORMAT ( 1 X , "ENTER  INITIAL  AND  FINAL  X  VALUES 


,  parallel 
<-") 

,*  INITIAL" 


T  0  »  , 


«•") 

«-") 


U  U  U  (J 
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?"  Y  and  l  VALUES;  AMO  THE  RADIUS  F  OR  EACH  BONE  (CM):  " ) 

11?  F  ORMaT  ( 1  X  ,  "NONE  «  ”,11,":  «.'') 

115  FORMAT (IX, "ENTER  LINEAR  ATTENUATION  COEFFICIENT  OF  BONE", 

2 "  PUR  DESIRED”/"  GAMMA  KAY  ENERGY  ( CM—  1  3  S  "  ) 

1H  FORMAT  (lx,  "SOURCE  ACTIVITY  RPR  UNIT  GRID  (PHOTONS/UNIT  TIME)”, 
2”:  <-") 

lib  FORMAT  ( I  X  ,  "GAMMA  RAY  ENERGY  OF  SOURCE  PLANE?  CKEV):  «-") 

L 

C 

93  FORMAT  ("1",  l^X,  "SCAN  *  " , I  A, 2«X , 9A2 , 5 ( / 1 X ) ) 

200  FORMAT  ("0" , "THt  CRYSTAL  RADIUS  IS  ",F6,3,"  CM") 

2012  FORMAT ("M"," THt  COLLIMATUR  SHAPE  IS  RFC T ANGUL AR " ) 

2021  FORMAT ("0" , "RECTANGULAR  COLLIMATOR  PLANE  COORDINATES:", 

22C/1X)  ,  1MX, "INSIDE  PLANES!  m  M  P  K") 

2022  FORMAT ("0" ,22X, "G *  " , 4  (F6 t 2 , 2X ) , / i X , 22X , " H i  "4 (F6.2, 2X)  , 
2/lX,22x,"i:  ",4(F6,2,2X) ,/lX,22X,"j:  "4(F6,2#2X)) 

2023  FORMAT  ("10  ",9X,  "OUTSIDE  PLANES:  M  N  P  K " ) 

2024  FORMAT ("0" , "COLLIMATOR  HEIGHT :", F6,3 , "  Cm;","  CRYSTAL", 

2 "  END  LENGTH: " ,F6, 3, »  CM;","  w IDTH X " , F6 , 3 , "  CM;",/, 

3"  FACE  END  LENGTH!  ",F6.3,»  CM;","  W  I OTH  :  " , F 6 • 3 , "  CM") 

2032  FORMAT ("a" , "THE  FACE  END  THICKNESS  OF  THE  COLLIMATOR  IS  ", 
?F(5,2,"  CM") 

204a  FORMAT ("0" , "LINEAR  ATTENUATION  COEFFICIENTS  USED" 


2  0  4J 

2042 

2043 

2044 
205 


206 

2D? 

2U8 

20  9 
2)  0 

21  1 

212 

214 

21b 
2b0 
2b  1 


3P0 
3  0  I 
50  0 


2 , "  (Cm- i)  :  "  ) 

F  OK M  AT  ( 1  X  ,  26X  ,  "  A  IK  :  "  ,  F  9 . 6  ) 

FORMAT ( IX , 23X , "TISSUE  :  "  , F9.fcO 
FORMAT  (1 X , 25X , "BONE !  ”  , F9#5) 

FORMAT (IX, IPX, "COLLIMATOR  MATERIAL:  ",F9.6) 

FORMAT ("0" , "TISSUE  SURFACE  -  CRYSTAL  FACE  SEPARATION:  ", 


2  F  6 , 3  , "  CM") 

FORMAT ("0" , "X , Y , Z  COORDINATES  OF  SOURCE  PLANE  PT.  Pi:  ", 
23F8.3,"  (CM)") 

FORMAT ( 1 X , 3  8  X , " P  2 :  ",3F8.3,"  (CM) ") 

FORMAT (1X,38X,"P3 I  ”,3F8,3,"  (CM)") 

FORMAT  ( 1 X , 38 X , "P4 S  ",3F8.3,"  (CM)") 

format ("0" , "THERE  ARE", 12,"  BONES  BETWEEN  SOURCE  AND", 

2”  CRYSTAL  TO  ACT  AS  ABSORBERS  AND  SCATTERERS") 

FORMAT ( "H" , " THERE  ARE  NO  BONE  ABSORBERS  USED  IN  THIS", 

2”  SIMULATION") 


FORMAT  (IX,"  BONE  #",I1,"  X  I  !  "  , F  7 , 3 , "  XF:",F7, 3, 

?"  Yi:», E/,3,"  ?  I J " , F  7 , 3 , "  RADIUS: " ,F7. 3, "  (CM)") 

FORMAT  ("Yl",  "SOURCE  ACTIVITY  PER  UNIT  GRID  ARtA  USEOS", 

2  F 1 H  f 1 ,  "  P HU TONS/UNIT  TIME") 

P ORM A f (" 0 "," SOURCE  GAMMA  RAY  ENERGY  USED : " , F b , 2 , ”  KPV") 

F  H  R  M  A  T  (  "  1  "  ,  "  "  ) 

FORMAT  ("0"  ,  "EXECUTION  STARTED  AT:  ”,I3,"  HRS  ”  ,  I  3  ,  MINS 
2X3,  "  SF.CS"/IX,10X, "STOPPED  A  I  :  "  ,  1 3 ,  "  HRS  '',13,”  OIOS  , 

313,”  SECS") 

FORMA T ( / ) 

FORMAT (IX, 16) 

FORMAT (IS) 


r 


c 


INITIALIZING  DETECTOR  POINT  PARAMETERS 


CALL  EXEC  C 1  1  ,  HD 
I MC -CRR / d  *  2  +  0*5 

I N  U  =  b . 3  *  i  N  c / 1 3  *  0  +  0*5 


r'r 


n  o  n  O  r1 
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XINCalNC*l ,0/2.0 
1 INC  a  INC/ 2 
XINU=INU*1 .0/2.0 
II  NO = I NO/ 2 

IF  UINC.Lfc.IINC)  I,NC  =  lNC+t 
IF  (XlNU#L£f»  1 1 N  U )  I  N  0  =  I  N  U  +  1 

c 

irwxx»iNC/2 
I  C  R  Y 1  =  1 NO/2 
i>2.  .  O 

Z  C  =  0 , 0 
X  C  a  1  d  0  „  0 
Y 0*100,0 
C 

c 

c  initializing  source  point  parameters 

c 

c 

N T o s  (SOW  rCCSX2"*Sxn**2+(SY2-SYl)**2  + 

2  CSZ2-SZ1)  **2n/0.5  +  0.5 
NTW=  tSNKl  C (SX3»SX2) *  *  2  +  (SY3-SY2) **2  + 
2(SZ3-SZ21**2))/1.0>  +  0.5 
NKLSNYL+  1 
N  K  W  *  M  T  W  ♦  1 
C 

c 

DO  9  Rw=l  rNKw 
lwsCrsW~lj*l,0/NTW*It0 
XlaSXl+(SX4»SXl)*Th 
Y  I  =  S  f  1  +  (  S  Y  4  -  S  Y  I  )  *  T  W 

ZJ aSZH-(SZ4-SZl)*Tw 
X2=SX2+ (8X3-SX2) *Tw 
Y?*SY2+  CSY3-SY2) *TW 
Z2SSZ2+  (SZ3-SZ2)  a'Tw 
1S®0 


i:>  O  10  K  L  *  1  »  N  K  l 

11.. s  ( k L -  1  .)  »i«0/NTL*l«0 

SX*X  1  +  (X2-X  1J  *TL 

SYSYi+tY2-YljaTL 

S7.  =  Z1  + tZP-ZU  *Tl 

WHO » SUNT C  C  S  X - X  C  5  *  *  2  + ( S  Y  »  Y  C ) *  *  2 ) 
CALL  SOL (CNR, SZ, WHO, 36, SOLAN) 
ISaIS+1 


J  e  -  1 
J  J  =  1 
$  P  =  0  ,  W 

jn  =  o 

w  T  s  0  •  0 

C 

C  main  CONTROLLING  segment 

c 

t )  n  5  o  i  x  =  1  ,  t  N  C 

oxs(~icwxi  +  ux-m*0.4  ♦ 

IFLJJ.GE.IN0)  GOTO  7 
J  5=  J  +  1 


i 1  =  (1CRYI  +  J) *2  +  1 
IF(Ii,LEtINC)  GOTO  6 
J  ,T  =  J  J  *  1 
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7 

fS 


c 


c 


c 


4  1 

41/: 

42 


C 

8 
b  W 
C 
C 
L 

c 


17 


C 


J  *  J  -  1 
11=11-2 
on  8  i y a i,n 

U  Y  =  (  -  1 1  C  R  v  .1  +  J  )  +  ClY-l))*Wf4  +  100,0 

j  n  =  j  o  i 

AsUX-XC 
ti  *  U  Y  -  Y  C 
C=UZ-2C 
At  P  =  S  X  -  0  x 
B  F I  =  6  Y  -  0  Y 
G A  h  =  S  2  ** u 2 

*  =  G  A  -1  /  tAuP**2+BFT  *  *  2  +  G  A  M  *  *  2  ?  **\  »  5 

A  T  s  «  +  W  T 

0  s  W  *  S  0  L  A  N  *  S  1  G  M  A 


IF  C  N  B ,  E  Q  »  0 j  GO  i  n  40 

8  H N A  I  so, 0 

XB ON  1=0,0 

1.)  n  4  1  J  B  S  1  ?  N  B 

HRs  I.)  r-YBUNFCJB ) 

»r=DZ-Z»UNfc(JB) 

R  B  =  W  8  0  N  F  u  J 
CALL  B  I  H 

b  0  in  A  T  a  t  X  P  C  —  M  U  B  U  N  a  B  U  N  1  H  )  +  BO  N  A  T 

a  RUN i sxBONT  +  BONTm 
cnu r i n  u  e 
Gniu  4  2 

d  H  N  A  f  a  1 , 0 
X R 0 Ml  a 0 , U 
CAUL  T1H 

1  T  S  A  I  =  h.  X  P  t  -  MU  T  J  S  *  T I  ST  H) 

CALL.  H  I  H 

cnu  A  i  =  LX  P  C-riuCUL*C tit  T H ) 

C  A  L  L  A  r  H 

A  I  K  A  r  =  fc  X  P  t  -  MU  A  1  R*  A  I  I  H  ) 

AF  1  N  a  Q  *  ft  U  N  A  T  *  T  ISAT*COLAT*AIRAT 
OR  UU)  =AFIN  +  0 R  (  J D I 
5  P  =  A  F  I N  +  5  K 

COinT  INUL 

c  n  n  t  i  n  u  fc 

Sun  OVER  INDIVIDUAL  MATRIX  INDICES  AND  DUMP 

UN  PAPER  tape 

IF  CJU.GT ,  b  0 )  SI  UP 
on  1 /  J  X  a  1 , J D 
inw C JX  J  a OR CJX) /A  1 
c n N TIN U £ 
iSK  t I S  J  a  S R  7  W T 


c 


if  usr (is)  .U:,y )  Guru  in 

m  R  ]  1  E  1 6  ,  A  o  1 3  I  b  K  C 1 3 ) 


159 


X  U  CONTINUE 

W  P  I  T  t  (  4  f  503)  (  1  SR  c  1  1 1)  ,  I  IT  *  1  ,  NKL  ) 

9  C  0  tv  T  l  N  U  E 

C 

1  T 1 1  4 , 5/1 .1  ]  a  OK  c  1  I  n  ,  II  I  =  1  ,  JO) 

CALL  E  X  £  C  ( 0  ,  1  0  0  4  B ) 

CALI.  E  X  E  C  U  1  ,  IT  2  ) 

C 

!PI  f  r.  (6, 25.n  I  T  1  (4)  ,  I  T  1  (S)  ,  H  H2)  ,  IT2(4)  ,  I  T2(3)  ,  1T2C2) 
w  R  1  r  £  ( b  ,  3  w  0 ) 

STOP 
L  N  0 

c 

c 

C  SOL  10  ANGLE  subroutine 

c 

c 

SUORU U TINE  SOL CR, H, RHQ, NC, SOLAN) 

DIMENSION  EC90)  ,  I -  (90) #  G  (  9  0 ) ,0(90) , EM (90) , 

2  Y (90)  , X (90  ) 

N  =  N  C 

P  T  =  0  .  1415  9  26 
N2=(N-2)/2 
X  N  *  N 
N  4  s  N  /  4 
N  6  =  (  N  4-  <>  )  /  4 
N  4  4  4  =  (  N  +  4  )  /  4 
b 0=SXNtPX/XN) 

•  J  2  i  =  N  2  +  1 
N?  4  = (N-2) / 4 
444= ( N« 4  J / 4 
N?4 1 =N24+ 1 

SRR  =  (Pl/XN) *  (BO/COS  CPX/XN) ) 

IF  (SRR„LT.0,O)  STOP  1 
Aps?,  *R*SURT (SRR) 

c  determine  IF  N/2  IS  onu  OR  EVEN 
XP2=n/4 

X  M  N  s  X  N  -  4  .  H  *  X  N  2 
IF (XNN) 20 #10, 20 
C  ODD r EVEN , OUU 

L  CALCULATION  of  D(IR)  IF  N/2  even 

10  00  IS  I R  a  1.  #  N  4  4 

X  T  K  u  X  R 

1  S  u  (  xR)  =XP/2„0+XP*SlN  (PI/4,  U-PixX]R/XN)  * 

2SIN(PI*XlR/XN+PI/4,0-Pl/XN)/en 
0  (  N  4  )  =  X  P  /  2  ,  kl 
DO  14  XR*N444rN21 
I P  K  s  (  N  +  2  )  /  2  -  1  R 
14  0 (I R ) = - 0 ( 1 R  R ) 

GOTO  00 

C  CALCULATION  OF  D(Ik)  IF  m/2  ODD 

23  DO  21  I R  = ! , N 2 4 
X  I  R  =  X  k 

21  i)  f  I  R  )  =  X  P  *  S  1  N  l  P  I  /  4 , 0  -  P  I  /  C  2  .  U  *  X  N )  ♦  P  1  *  <  I  R  /  x  5  * 

2 SIN  (PI/4„0'“PI  /  (2,0*XN)  +  PI*X.tR/XN)  /  BO 

0(0241) =0,0 
on  22  IRsN(),N21 
j.  P  R  =  (  N 1 2  )  /  2  -  I  R 

2  2  0  f  I  R )  s  - 1)  ( I  R  R  ) 

00  continue 


•  < 


r~.  o 


160 


c 

c 


3  i 


4  D 
3  0 


ini 
4  I 


5  1 
4  5 


4  2 


45 


!>n  31  1R=1,M2 

X  T  K  S  I R 


2 


t  )  -  SIN  «,  (  (  X  N  •  4  ,  W  *  X  I R  )  *  P  I  )  /  (2,0*XN)  1  / 

LnSUUN-4#w*XlW)*Pn/(2,0*XN)) 
on  35  IR=i,n2 


x  T  k  =  i  K 


i  Mk)-(XP/2,  \a  )  *  5  i  m  (  ( 2 , 0  *  X  I  k  -  1  .  tl )  a  p  1  /  x  N  )  /  b  f  i 
on  30  1R=i,n21 
X C IP) =RHO~D C Ik) 
on  4  1  I R  =  1 , N  2 

I  i 1 H  s  H  *  $  Q  R  T  t  X (Ik) *  *  2  +  Y (IR)**2+H**2) 

WPI  stMCIR)*(XUR)**2^H**2)-XCTR}*Y(I«) 

1  r Ik) = AT  AN  (TOP /BUT) 

1 1 ( F (lk) ) 50,41,41 
^IRjsFf  Ik)  +  p  j 

Cn.o  T  inuE 

011  4  5  IK  =  1,N2 

3kk  =  A(ik+l)^*2tY(,  Ih)**2«2.0*fcHUK)* 
?X(iKJ*y(IR)+(fM(IR)**2)*(x(lw)**2)+H*H+ 

■3 c e  n (ikj  **2j *cx  ctK+n **2)+2,m*[;kir) *x fik+n  * 
4  y (lk)-2,0*CkM(iR)**2)*x(  ( k  + 1 ) *  X  (1 R ) 

IF (Skk ,  LT*0.8)  SI  OP 2 
i OP=H*SQR T  (SRk ) 

6  n  T  =  E  M  1 1  H  )  *  (  H  ★  H  +  X  (  1 R  +  i  )  *  X  (  I K  )  )  -  x  (  X  K  +  n  *  Y  ( I R  ) 
t f lk) sATAmiTOP/BOT) 

IF (F (IR)I5l,45,4t> 
enk)s£UR)+pi 
i  n  .0  !  i  N  u  t 
on  4  2  I  k  c 1 , N 2 
b  HRJ  s£  UR)  -F  (  XR) 

S  n  L  A  N  s  14 . 0 
00  4  3  1  R  =  l  ,  N  2 
vSOL  A'MsSOL  AN  +  2  *  ( lk) 

k E turn 

t  N  U 


BONE  THICKNESS  SUBROUTINE 


c 

SUBROUTINE  b  T  m 

real  L  1  ,  L.2  ,  MG  T  ,  NG  T  ,  KG  !,  MH  T  ,  NH I ,  k  H I ,  M 1 1 ,  N I  I ,  K  II  ,  MJ I ,  NJ  I ,  K  J  I 

real  mgo,  ngq, kgo, mho, mho, kho, mio, nio, k  iq, mjjq, njq, k jo 

REAL  MU A  Ik , NUT  IS, MU BUN, MU COL 

COMMON  OX#()Y#DZ,SX,SY,  SZ,  A,b,C,  ALP,  BE T , G AM , T , R 1 , R2 , H , Rb , T SCFS 
COMMON  COL  I  El  ,  A  I R  X  H  ,  T  I  S T H  ,  BUNT H ,  MG  I ,  IMG  I ,  PG I  ,  K G I ,  MGO  ,  NGO  ,  PGO  ,  K GO 
COMMON  MH I , Nh I , P  H  T , K  H  I , MHO , NHO , PHO , KHO , M 1 1 , N II , P 1 1 , K 1 1 ,  M 1 0  ,  N 1 0 
COMMON  PIU, K 10, MJ  I  , NJ  T , PJ I , K J I , MJO, NJO, PJO, K JO, XbONT,  BB,  bC 
i FST = (Rb**2) * (dEI  **2  +  GAM**2  WBET*BC-GAM*BB)  **2 
IE  (TEST)  1,1,2 

1  ci  n  N  T  H  =  0,0 

RE  I  URN 

2  b  0  N  T  H  =  2  „  0  *  3  Q  R  T  (TEST ) /  (BE r**2+GAM**2) 

RE  1  URN 

E  N  0 
C 

c 

C  TISSUE  THICKNESS  SUBROUTINE 

C 

c 


, 


1 
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SUbKUU I  I Ufc  TTH 

'* r  A  L  1  • 1  »  L  2  i  M  G  I  ,N(»IfK6I#MHIfNHIfKMl,MII,NIT,KII,MJI,NJI,KJI 
i<F.  AL  MGU,NGQ,KGO#MhO,NHO,KHO,MIO,NIO,KIO,MJjO,NJO,KJO 

>1 E  A  L  M  U  A  .l  k  ,  r  i  i  j  I  1  S  ,  iv|  y  0  0  N  ,  M  U  C  0  L 

LiMil'Jf,  0X,uY,l)Z,SX,SY|SZ,A,B,C,AlP,BET,GAM,rrRl,k2,N,Rb,TSCFS 
C 0 r-i n U N  COL  I  H,  A  1RT  H,  TISTH,  BONTH,  MG  I ,  MG  I  ,  PGI ,  KG  I  ,  MGQ,  NGO,  PGO,  KGO 
COMMON  MH I , Nh I , PH  I , KH I , MHO , NHO , PhO , K HO , M XI ,  M II , P  X  I , K 1 1 , M 1 0 , N 1 0 
COMMON  PIO, K IQ, MJI ,NJI, PJ  I  #  KJI #  MJ0,NJO, PJO,K J0#X80NTfBB,bC 
IF tTSCFS.LTB$7J  GOTO  1 
f  1ST  M  a  v  i  .  i/ 
k  p.  I  U  '<  N 

1  AV*UX+ALP* CTSCFS-DZ) /GAM 
AVsOY  +  ttfc 1  *  CTSCFS-DZ) /GAM 

bPK= ( CSX-AXT  **2+ (SY-AY) **2  +  (SZ-TSCFSj **2) 

T I STrt=SQRT CSkW J -X80NT 
lev  FORMAT  (1X*E10,6) 

fv  p  I  U  R  N 

L  M  >,) 

c 

c 

C  RECTANGULAR  COLLIMATOR  SUBROUTINE 

C 

C 

SUBROUTINE  RTm 

K  FA  L  LUL2f  MGI,NGI,KGI#MrtJf  NHI,KHI,MII,NII#KII#MJI#NJI#KJI 
WEAL  HGO  ,  NGO  ,  K  GO  ,  MHO  ,  NHO  ,  K  Htl ,  M  I U  ,  N 10  ,  K  1 0 ,  M  J  0  ,  N  J  0  ,  K  JO 
COMMON  OX , UY , OZ , SX , SY , SZ , A , B , C , ALP , BET , GAM , T , R 1 , R2 , H 
COMMON  CUL  r  H , A  I R T  H , T  I  S  TH , BOOT  H , MG  I , NG I , PG I , K G I , MGO , NGO , PGO  ,  K  60 
COMMON  MH I , NH I , PH  I , KH I , MHO , NHO , PHO , KHO , H 1 1 , N I  I , P 1 1 , K 1 1 , M 1 0 , N I  0 
COM  MON  P 1 0 , K 1 0 , M J 1 , M J ! , P J I , K  J I , M JO , N J  0 , P JO , K  J 0 , TSCFS 
L)  F  N  s  M  H  I  *  A  L  P  +  N  H  I  *  B  E  T  T-  P  H I  *  6  A  M 
IF  tU£N,EUf  W.fti)  GOTO  1 
I  HI=  UhI«*MHl*PX-NMI*UY-PHI*OZ)  /DF.N 
C 

EZH=UZ+GAn*THI 
IF  CEZHT  1,1,2 

\  uENsMJ  J  *ALP  +  NJI*BET+PJI*GAM 
IF CDEN.EQ, 0f0)  GOTO  6 

rjI»tKJli-MJI<fOX**NJI*OY»PJI*OZ)/OEN 

tZJ=DZ+GAM*TJ I 
I F  U"  Z  J  )  6 , 6  ,  3 
C 

2  D  F.  o  =  i  .1  i  *  A  L  P  +  N  J  I  *  B  F_  1  +  P  J  I  *  G  A  ri 
IF (DEN. £0.0*0)  GOTO  4 

I  ,11s  ckji-mji*l)x-nji*uy*,pji*oz)  /o£n 

E  Z  J  s  0 1  +  G  A  M  *  T  J  I 
I F  t  a L J) 4 , 4 , 5 

H  L  F  C  t  Z  J  -  E  Z  H  J  3  f  3  /  4 

0  I FPAHa (EZ J-UZ) /GAM 
1  =  1 

GOTO  lv 

4  I FPARs (EZH-UZ) /GA < 

1=2 

GOTO  10 
C 

c 

10  EY =OY+BET*TEPAK 

lFtY2*LE*EY*ANU.tY.Lt.Y5j  GOTO  11 
GOTO  b 


C 
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1  1 


C 


6  0 


C 


12 


1.5 


C 


1  4 


1  6 
1  / 
2  1 


IS 


1  8 
22 


2  4 


C 


2  3 


2  b 


6 

34 


3  2 


C 


3  1 


tZ  »UZ  +  GAM*  TEPAK 
IF  t  E  L  -*  H 1  12,60,60 

COLT  Ms  Vi.*) 

K  F  T  U  P  N 

E  X  =  0  X  +  A  L  2  *  T  E  P  A  k 

lFU-i)  13,13,14 

UF.  N  =  0JU*ALP’H\iJU*BEl+PJO*GA0 

IFCOEN.EUtW.45)  GOTO  21 

1' J0  =  (K  JO-M  JO*DX*N  JO*pY-.PJO*DZ)  /  DEN 

F  Y  =  D  f  +  B  E  T  *  r  JO 

G  0  T  U  1  6 

j  f  ;  a  a  M  H  Q  *  A  L  P  +  N  M  U  *  B  E  1  +  P  H  0  *  G  A  M 

IF  lULN.EU,  W.45)  GOTH  21 

THUs  (KHO-MHO*nx-NHLl*OV-PHO*L>Z)  /  DEN 

F  V  =Oy+BET *  THO 

IF  i,Yi  .LE.FY.  ANU.FY.LE.  Y4)  GOTO  17 
G  0  T  LI  1  b 

r  7  ~  0  Z  ♦  U  A  H  *  C  F  Y  *L)  Y  )  /  b  £  ] 

IF  IF  Z-ri)  19,21,21 
r  V  =  L)X  +  AIP*  (H-iU)  /  GAM 
b  V  s  i)  T  +  U  E  r  *  C  H  -  0  L  )  /  G  A  4 
F  Z  =  M 
G  n  T  0  2  YJ 

F  v  sOx-FALP*  (FY-GY)  /BET 
GOTO  2  v) 

IFCOY-SY)  22,22,23 

i)EM  =  MXO*  ALP  +  NI0*BE  f +PIO*GAM 

IF CDEN.EG.0.W)  GOTO  21 

X  I  0  =  (  K  I  0  -  M  1 0  *  0  X  -  N  1 0  *  l>  Y  -  P  1  0  *  0  Z  )  /  D  l  N 

F  Z  =  DZ*GAim*T10 

I  F  C  F  Z - H )  2  4,21,21 

F  y  =  l)  X  +  A  L  P  *  T  T.  0 

p  y  s  j  )  y  +  Fj  E  T  *  T  1 0 

GOTO  20 

0 F N  =  1 G 0  a  A L 0 * N G U  * B E  T  +  P G U*G  A 3 

.1  F  l  i.i  t  h  .  P  0  •  YJ  ,  F  )  G  0  T  U  21 

f  G  0  s  t  K  G  0  »  M  G  0  *  0  X  «  N  G  u  *  D  v  -  P  G  0  *  0  Z )  /  0  fc  ^ 

F  7  a  D  l  +  G  A  M  *  T  G  0 
IFCFZ-H)  23,21,21 
F  X  =  0  X  +  A  L  P  *  T  G IJ 
F  Y  =  L>  Y  +  6  E  T  *  T  G  U 
b  0  1  U  2  O 


IF(OY-SY)  30,30,31 
I j  F  t  i  =  1 1 t  *  A  L  P  +  N I i *  B  E T  +  P  ll  * G  A H 
IF  COEN, EG. 0.0)  GOTO  00  % 

T  I  1  =  l  K  1 1  •*  O  1 1  *  0  X  »•  n  I  J  *  0  Y  •  P  II  *  0  Z  )  /  D  E  N 


t  7  5  U  Z  +  G  AM*  f  II 
IF CEZ-H)  32,60,60 
E  X  =  0  X  +  A  L  P  *  T  T  i 
t.v  =0  t  +BE  1  *  T  I  I 

IF  t  A L P  J  13,13,14 


U  f  ;j  =  n  G  1  *  A  L  P  +  N  G  I  *  B  £  I  +  J  G I  *  b  A 1 1 


0  0  0  0 
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IF  IDEn.EQ.M.M)  GOTO  60 
TGlstKGI-M6I*0X-NGI*DY-PGI*DZj/DfcN 
L  7  =  L)  Z  +  G  A  M  *  f  G 1 
IF  (FZ»N)  3 3 , 6 0 , 6 0 
33  t.y  SOX  +  ALP*  TGI 
fc.vnorfriET*TGl 
IF  (  ALP)  13,13,14 
C 

c 

2/'  C  0  L  T  H  =  S  Q  K  T  (  (  F  X  -  E  X  1  *  *2 +•  ( p  Y  -  F"  Y  )  *  *  2  +  CF  2"EZ)  *  *  2 ) 
k  F  T  U  K  N 

i  N  U 

c 


AIK  THICKNESS  SUBROUTINE 


SUbWUUTlNt  A 1 H 

K E  AL  L 1 , L 2 , MG  1 , NG I , K G I , MH I , NH I , K H I , M 1 1 , N II , K 1 1 , M J I , N J I , K J  I 
W  FA  L  MGO  ,  NGO  ,  K  GO  ,  OHO  ,  MHO  ,  K  HO  ,  HI  0  ,  N 1 0  ,  K  1 0 ,  M  J  JO  ,  N  JO  ,  K  JO 
KEAL  MU  AJK  ,  MU  I'lS,  MuBON,  MUCOL 

COMMON  OX  ,  l)Y  ,  PZ,  SX,  SY,  SZ,  A,B,C,  ALP , BET , G AM , T , HI , K2 , H , RB , TSCFS 
COMMON  COL  1  H,  AIHTH,  T I  ST H , 80N T H , MG  I  , NG I , PG  I  , K G T , MGO , N&O , PGO , K GO 

COMMON  MH 1 , NH l , PH  I f  KHI , MHO , NHQ , PHD , K  HO , MI  I , N 1 1 , P 1 1 , K 1 1 , M I  0  ,  N 1 0 

COM HON  P 1U, K 10, Mj I , NJI , PJ I , K J I , M J 0 , N J 0 , P J 0 , K J 0 , X 80NT , 88 , BC 
0IS=6QRT (AlP**2+BET**2+GAM**2) 

AlKTHsOIS-TISTH-BONTH-COLTH 
K  F  I  U  R  N 
END 
Emu  ,t> 


. 


